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Study on stationary solutions to the problem of phytoplankton dynamics considering 
transformation of phosphorus, nitrogen and silicon compounds 


Yu.V. Belova’, A.M. Atayan’, A.E. Chistyakov’ , A.V. Strazhko* 


'*34 Don State Technical University, Rostov-on-Don, Russian Federation 


McceqoBanne CTalMoOHapHelx pelieHHh 3a24a4n AHHAaMUKH PUTONIAHKTOHA C y4eTOM TpaHcwmopMallnn coeqn- 


7 OK 2k 
HeHHH Mocwopa, a30Ta HW KpeMHHA 


IO. B. Besopa’, A. M. Aras’, A. E. Wucraxos’, A. B. Crpaxko*” 


1,2,3,4 


Introduction. The solution to the problem of transformation of 
phosphorus, nitrogen and silicon forms is studied. This prob- 
lem arises under modeling phytoplankton dynamics in shal- 
low-water bodies including the Azov Sea. The phytoplankton 
dynamics model is formulated as a boundary value problem 
for the system of diffusion-convection-response equations and 
takes into account the absorption and release of nutrients by 
phytoplankton, as well as the transition of nutrients from one 
compound to another. To calculate the initial conditions and 
parameters of the equations under which the steady-state re- 
gime occurs, the software is developed, which is based on the 
model describing changes in phytoplankton concentrations 
without considering current effects. This model is represented 
by a system of inhomogeneous differential equations. Based 
on the developed software, the initial conditions and parame- 
ters of the phytoplankton dynamics model in the Azov Sea are 
calculated experimentally. 

Materials and Methods. A 3D model of phytoplankton dynam- 
ics 1S considered taking into account the transformation of 
phosphorus, nitrogen and silicon compounds based on the 
system of nutrient transport equations. The case of a spatially 
uniform distribution of substances is considered to specify the 
parameters of the model at which the stationary modes occur. 
Because of simplification, a system of ordinary differential 
equations solved through the Runge-Kutta method is obtained. 
Research Results. The software is developed to specify the 
initial conditions and parameters of the phytoplankton dynam- 
ics model considering the transformation of phosphorus, nitro- 
gen and silicon compounds. 

Several numerical experiments are performed under the as- 


sumption that the development of phytoplankton is limited by 


” The research is supported by the RSF (project no. 17-11-01286). 


JloHcKoH rocyapcTBeHHbIM TeXHuYeCKHH yHuBepcuterT, r. Poctos-Ha-Jjony, Poccuitckaa Dezepaina 


Beedenue. Padota TocBalleHa HUCCeOBAHHIO pellleHHsA 3ajla- 
4H TpaHchopmamuu dopm docdopa, a30Ta u KpemHua. Jan- 
Had WpoOsiemMa BO3HHKaeT MpH MOJeuHMpoBaHuu jWWHaMUKU 
(UTOMIAHKTOHA B MCJIKOBOJHBIX BOJOCMax, B TOM 4HCIIe B 
A30BCKOM Mope. Mojenb TMHaMuKH (UTONMAHKTOHA cop- 
MyJIMpoBaHa Kak KpaeBad 3ayjada JJId CHCTeMbI ypaBHeHuit 
muddy3uu-KOHBeKUMH-peakIIMH WU YAMTHIBAeT NOTMOMWeHve U 
BbIJ[eICHHe MMTATeCIbHbIX BeEL[CCTB (PUTOMAHKTOHOM, a TakKxKe 
TlepexOJ,| MWHTAaTeCJIbHBIX BeIIeCTB H3 OAHOTO COeMHeEHHA B 
mpyroe. Jia pacueTa HauasIbHbIX yCJIOBHN UM MapaMeTpoB 
ypaBHeHHi, Mp KOTOpbIxX HacTyMaeT CTalMOHapHbIN pexuM, 
pa3paOoTaHo TporpaMMHoe oOeciieyeHve, OCHOBOM KOTOPOTO 
NOcHyKUIAa MOJelIb, OMMChIBAIOMIad H3MCHeCHHA KOHIICHTpa- 
WH MUTOWIAHKTOHa Oe3 yueTa BIIMAHUA TeyeHn. JlaHHas 
MOJeIb IIpecTaBIeHa CHCTEMOM HeEOJHOPOAHBIX OOBIKHOBCH- 
HbBIX TuddepeHiwalbHbIX ypaBHeHui. Ha ocHoBe pa3pavo- 
TaHHOrO TporpaMMHoro oOOecre4eHHA IKCICPHMCHTAJIbHbIM 
00pa30M paccuHTaHbl HayaJIbHbIe YCIOBHA HU MapaMeTpbI MO- 
eu JMHAaMUKH (PUTOMIAHKTOHA B A30BCKOM Mope. 
Mamepuazoi u memoovi. PaccMaTpuBaeTcA TpexMepHad MO- 
eb TMHAaMUKH (UTOMIAHKTOHa C y4YeTOM TpaHcdopmalnu 
coeyMHeHHi docdopa, a30Ta HW KPeMHHA, OCHOBaHHad Ha CH- 
cTeMe ypaBHeHHui TpaHcnopTa OMOreHHbIX BelllecTB. Jia 
YTOUHeCHHA WapaMeTpOB MOJeIu, pw KOTOPbIX HacTyMaroT 
CTal|MOHapHble pexKUMBI, paccMaTpHBaeTcaA Clryyait IpocTpaH- 
CTBCHHO-OJHOpOHOTO paciipeszeneHua cyOcTtanHunui. B pe- 
3yIbTaTe YMpOleHuA MONyyeHa CuCTeMa OObIKHOBCHHBIX 
muddepeninalbHbIx ypaBHeHHH, KOTOpad pellieHa MeTOOM 
Pyure-KyTTbl. 

Pe3yibmamol ucciedoeanus. Pa3zpad0oTaHoO mporpaMMHoe 
oOecnieyeHve Jd YTOUHCHHA HavasIbHbIX yCJIOBMM U Tapa- 
MeTPOB MOJIeJIM JHHAaMUKU ()UTOMIAHKTOHa C yYeTOM TpaHc- 
(bopMallMu coesqMHeEHuM Mocdopa, a30Ta uv KpemHua. IIpose- 
JICHbI HECKOJIBKO YMCJICHHbIX IKCIICPUMCHTOB B TpeTO0%*Ke- 
HHH, YTO pa3sBHTHe MUTOMIAHKTOHa JIMMUTUpyeTcA eMH- 
CTBeCHHbIM OHOreHHBIM BelIeCTBOM. B pe3ysIbTaTe BbIYHCJIM- 
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a single biogenic substance. As a result of the computational 
experiment, it can be seen that with the obtained values of the 
initial concentrations and parameters of the equations, station- 
ary modes occur for the system of ordinary differential equa- 
tions describing the case of the spatially uniform distribution 
of substances. 

Discussion and Conclusions. The mathematical model of the 
transformation of phosphorus, nitrogen and silicon forms in 
the problem of phytoplankton dynamics is studied. Stationary 
modes for the system of ordinary differential equations are 
obtained, for which the values of the system parameters and 
initial conditions are determined. The results obtained can be 
used in further simulation of the phytoplankton dynamics con- 
sidering the transformation of phosphorus, nitrogen and silicon 
compounds with account for convection-diffusion, salinity, 


and temperature. 


Keywords: phytoplankton, phosphorus, nitrogen, silicon, bio- 
gen, chemical-biological source, convection-diffusion- 
response equation, Cauchy problem for system of ordinary 
differential equations, stationary mode. 


For citation: Yu.V. Belova, et al. Study of stationary solutions 
to the problem of phytoplankton dynamics considering trans- 
formation of phosphorus, nitrogen and silicon compounds. 
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TCJIbHOrO IKCIECPHMeHTAa BUHO, 4YTO IPH MOJYICHHBIX 3HAave- 
HHAX HaYaIbHbIX KOHI[CHTpallMu WU MapaMeTpOB ypaBHeHuii 
HacTyliatoT CTallMOHapHble Pe?KHMBI JIJId CHCTCMbI OObIKHO- 
BCHHBIX uddepeHuMaNbHbIX ypaBHeHH, OMMMCbIBAFOIIeH 
ciydai lpOCTPpaHCTBeHHO-paBHOMepHOro _ paciipewesIeHust 
cyOcTaHIUni. 

Oécyacodenue u 3aknrouenus. B padote uccieyqoBaHa MaTeMa- 
THYeCKad MOJIelIb TpaHcdopmMalnu dopm docdopa, a30Ta u 
KPeMHHAl B 3ayjade JMHaMMKH (UTOTAHKTOHA. 

IlonyueHbI cTal[wOHapHble pexKUMBI JJId CHCTeMbI OOBIKHO- 
BeCHHBIX JM@depeHHaIbHbIX yPaBHeHHH, JIA KOTOPbIX ompe- 
JesIeHbI 3HAYeCHHA TapaMeTPOB CHCTeMbI HW HadasIbHble yCJIO- 
Bua. [loryaeHHble pe3yIbTaTbI MOFyT ObITb HCIIOJIb30BaHbI B 
ipollecce MasbHeiiiero MOeIMpoBaHuA JMHaMHKH MuTO- 
IIIaHKTOHa C yueTOM TpaHcdopmallun coeyuHeHUuu docdopa, 
a30Ta HW KPpeMHHA C YYeTOM KOHBEKIMU-7uddy3uu, COeHO- 
CTH, TeMMepaTypbl. 


Ku1r0ueBble C10Ba: (UTOMIAHKTOH, Mocdop, a30T, KPeMHUH, 
OvOreH, XHMUKO-OMOJIOrMYeCKHH HCTOUHHK, ypaBHeHie KOH- 
BeKIMH-Tuddy3uu-peakiun, 3ayaya Komm ia CHCTeMBI 
OObIKHOBeCHHBIX TH@MepeHiwaIbHbIX ypaBHeHH, CTallMoHap- 
HBIM pexKUM. 


Oopazeu Ona uumupoeanua. benosa, tO. B. UccneqoBpanue 
CTaI[MOHapHBIX PelIeHHi 3aaqW TWMHAMUKU (PUTOTJIAHKTOHA C 
yueTOM TpaHcdopmMaliwu coeqwHeHuH dochopa, a30Ta u 
kpemuusa / IO. B. benosa [u gp.] // Becrauk JJonckoro roc. 
TexH. yH-Ta. — 2019. — T.19, Nel. — C. 412. 
https://do1.org/10.23947/1992-5980-2019-19-1-4-12 


Introduction. Because of the development of major cities on the coast of shallow water bodies and river systems 


that flow into these water bodies, eutrophication has become more frequent. The growth of algae in reservoirs is caused 
by an increase in the flow of nitrogen and phosphorus compounds from the adjacent land areas. Each water body is 
unique and requires a thorough study. Field investigations [1] and mathematical modeling are used to explore water 
bodies. Without downplaying the role of field experiments, we can say that mathematical modeling is less costly, and it 
allows us to predict the behavior of the ecosystem. 

To study the Sea of Azov, a three-dimensional model of hydrodynamics [2, 3] including the equations of motion 
in three spatial directions was developed. In [4], this model was made for the case of dynamic rebuilding of the compu- 
tational domain geometry due to the tidal effects. The investigation of this model accuracy is given in [5]. In [6—8], the 
reconstruction of an ecological catastrophe that occurred in 2001 caused by an excessive concentration of algae in the 
eastern part of the Sea of Azov is given. In [9], methods of controlling the suffocation phenomena arising in the Sea of 
Azov were proposed. The [10—12] papers are devoted to studying the dynamics of phyto- and zooplankton. 

The water condition in shallow water bodies is changing rapidly, and mathematical models need to be refined. 
The parameters determination of the three-dimensional model of the phyto- and zooplankton dynamics is laborious; 
therefore, it is proposed to use a simplified model to calculate these parameters. 

The work objective is to improve the parameters of the model of the phytoplankton dynamics considering the 
transformation of phosphorus, nitrogen and silicon compounds, under which stationary regimes occur with the assump- 
tion of a spatially uniform distribution of substances. 

Materials and Methods. The model is based on the system of equations for the transport of nutrients [15, 16], 


the form of which for each F model block is 


Mechanics 
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i aera Oh Oh = div(keradg,)+R., (1) 
Ot ax oy Oz : 


where g,is concentration of the i-th component, [mg/l]; ieM, M={F), FP, F'3, POz POP, DOP, NO3, NO2, NHy, 


Si}; {u,V, w} are components of the velocity vector of the water flow, [m/s]; & is turbulent exchange coefficient, [m/s]; 


R, 1s function-source of nutrients, [mg/I - s]. 
In equation (1), i index indicates the type of substance (Table 1). 
Table 1 
Biogenic substances in model of phytoplankton dynamics 


PO, 
Pop | speded organic phosphor 


por solubie organic phosphor 
NO. 
NO, 
NE, 





Chemical and biological reactions are described by the following equations R, =C,(1—Kyr)dp —Krpdp —KreWs » 


3 
Roop = LSrK 04 — B pn pop NI nop > 
3 
ae — LS pK rer, + K pV pop 7 KS h@nop ’ 
3 
Reo, = LSC, a -1)q, + Koy pop FR pop ’ 


Ga. 3 
Fr (4y04no4on,) +28, (ke + Ky )4 —Kodua, » 


3 
Ruy, = LSC (K,,-1) 
BY (Ivor Inoe2 Ion.) Ivo 


=K, 
f N (Gro, ’ Ao, ’ dv, ) Ano, - no, 


234 no, » 


q, + K, 


2 nit, 


3 
Ryo, = L Snr (K,.. -1) 


(1) 
five = vs.C. (Kee 1) £2 fron) 1 no, 


| 9g +K 
= i N Ce »WIno,» Inu, ) Ino, * Ino, 


23 dno, ’ 


R, =5,C, Ce -1)q, +8 K np , 

Here, K,, is specific breathing rate of phytoplankton; K,,,,is specific die-off rate of phytoplankton; K,,, is 
specific excretion rate of phytoplankton; K,,, is specific rate of POP autolysis; K,,, is coefficient of POP phosphatifica- 
tion; K,,,, is coefficient of DOP phosphatification; K,, is specific rate of ammonium oxidation to nitrites under nitrifica- 
tion; K,,1s specific rate of nitrite oxidation to nitrates under nitrification; s,, s,, s,,are normalization coefficients be- 


tween the content of N, P, Si in organic matter. 
Phytoplankton growth rate is determined by the following expressions: 


Ci. = Kup, min fe (Gro, ) a Cee Inn, )} 


C;, Sg min fe Ce ee (edatde, ne (ds ) 4 


where Ky is maximum specific phytoplankton growth rate. 
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Fig. 1. Model scheme of biogeochemical transformation of phosphorus, nitrogen and silica forms 
Functions describing nutrient content 


1 po, 


— for phosphorus: f,(g,,. )= 
(dr Ino, + Ko, 


where K,,, 1s phosphates half saturation constant; 
1 Is; 
— for silica: f, (45,)=—_ 
9a t Ks 


where K,. is silica half saturation constant; 


— for nitrogen: f, (Gro, Teste ) = fe Ce »WIno,» Inn, ) +90) (Cae ) ) 


(1) (Ino, ale Ayo, ) exp(—K dT.) (2) 
tn CR Reet ee See Ix CE =e 
Kyo, + (no, + Ino, ) Kon Cap, 
where K,,, is nitrates half saturation constant, A,,, 1s ammonium half saturation constant, A,,, 1s ammonium inhibi- 
tion ratio. 
For the system (1), it is necessary to specify the vector field of the water flow velocities, as well as gq. initial 


values of the concentration functions 
q, (x, y,z,0) =G (x,y,z), (x,y,z) eG , t=0,ieM. (2) 
Assume % boundary of G cylindrical region is piecewise smooth and XY =X, UX, Uo , where &,, is surface of 
the reservoir bottom, X, is still water surface, o is lateral (cylindrical) surface. Suppose uw, is vector component of the 
water flow velocity normal to 2 , and 7 is outward normal vector to X . Forg, concentrations, we assume: 


— on side boundary: 
g,=9, ono,if u, <0,ieM; (3) 


Oq. 
AG, og, 


=0, on o,if u, 2>0,ieEM; 
On On 


—on X, there is reservoir surface: 


OTe G. gen (5) 


IZ 


—on the bottom %,, : 


Og. ; 
ke = 6,,9¢,. t€{Fi, Fo, Fs}, 
Zz 
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k at =«,,g,, i€ {PO4, POP, DOP, NO3, NO>, NHy, Si}, (6) 
IZ 


where &, ,, €,,are sedimentation rates of algae and nutrients to the bottom. 





Stationary Mode. Consider the case of a spatially uniform distribution of substances (phytoplankton, forms of 
phosphorus, nitrogen and silica); then each of the equations (1) is simplified; and as a result, we get the following sys- 


tem of ordinary differential equations (ODE): 


dq, a 
dt = C, (l— Keds a Kode — Kaede 0H 195 (7) 





oe 
—_ => s pK ED qr, = Ky drop — By pop» 


At 
—_— => S»K BE q a K, pd pop ~ or pop 9 


7 
——. m= 9's,C (K, 2a, +K, vA pop + K, 1 pop » 

Ayn, -¥s.C =) fy (NH,) Ys, (Key + Kx )4y —Kud 
i a OO ENO .NO.NH) " a 2 Ue eee 
dq 3 (NO,, NO,, NH) q 
Se Di SG. (Kaeo a 1) fx (NO,NO,»NH,) Gr g + Ky yn, —K3.9yo0,> 
dt i=l I (NO,, NO, , NH,) Ano, T Ino, 

d ") (NO,, NO,, NH 

do, =>s,C, (Ky, 1) NO NO NH) doy qd, EK Gn : 

dt i=l fy (NO,, NO,,.NH, ) Avo, * Ino, 





“ = 5,C, Ca -1) 4, gore ae ae 


We solve the system of ordinary differential equations by the Runge — Kutta method [15-17]. We will conduct 
several numerical experiments, assuming that the development of phytoplankton depends on a single limiting substance. 


Research Results. For the ODE system (7), we calculate the initial conditions and parameters of the equations 


at which the stationary regimes occur. Let us take the initial concentration values: g, (0)=2.5 mg/l, q,, (0) = 2.6 mg/, 
d,,(0)=0.91 mg/l, dpop(0)= 9.07 mg/l, yop (0)= 9.07 mg/l, Apo, (0) = 0.005 mg/l, qy,,(0)= 0.11 mg/l, yo, (0) 
= 0.0178 mg/l, qyo, (0) = 0.304 mg/l, g,,(0)= 0.4 mg/l; coefficients: K,,. = 2.8 day", K,,=0.15 day”, K,,,= 0.05 day” 
. Kg = 0.15 day" , Kpp = 0.015 day" , Kpy = 0.02 day" , Kpn = 0.1 day" , Kap = 0.9 day” , Ko3 = 2.5 day“, Kysi = 1.46 


day, sp=0.01, sy = 0.016, ss; = 0.023, Kyo, = 0.024, Ky, = 3.0, Ky, =2.0, K, = 3.0. 


The obtained stationary modes of the ODE system (7) on the assumption that the development of phytoplank- 
ton is limited by a single nutrient (phosphorus, nitrogen or silica) are shown in Fig. 2-4, respectively. Fig. 2 describes 
the effect of phosphorus on the development of various phytoplankton species; Fig. 3 describes the effect of nitrogen on 


the development of various phytoplankton species; Fig. 4 describes the effect of nitrogen on diatom development. 
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chy, ms 
2.6 

















Fig. 2. Stationary mode of ODE system under assumption that phytoplankton development is limited by phosphorus: 
a) green algae (ChV), b) green-blue algae (AF-A), c) diatom (SC), d) suspended organic phosphorus (POP), e) soluble 
organic phosphorus (DOP), f) phosphates (PO,) 
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Fig. 3. Stationary mode of ODE system under assumption that phytoplankton development is limited by nitrogen: 
a) green algae (ChV), b) green-blue algae (AF-A), c) diatom (SC), d) ammonium (NH,), e) nitrites (VO>), f) nitrates (VO3). 


sc, mg/l 
10 


Fig.4. Stationary mode of ODE system under assumption that phytoplankton (diatoms) development is limited by silica: 
a) diatom (SC), b) silica (Si). 
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The result of the computational experiment shows that with the above values of the initial concentrations and 
parameters of the equations, stationary modes occur for the ODE system (7), which describes the case of a spatially 
uniform distribution of substances. The obtained values will be used in further simulation of the spatially inhomogene- 
ous distribution of substances, saltiness and temperature considering the movement of the aquatic environment [18]. 

Conclusion. A mathematical model of the transformation of forms of phosphorus, nitrogen and silica in the 
problem of phytoplankton dynamics is studied in the paper. The case of spatially uniform distribution of substances 
(phytoplankton, forms of phosphorus, nitrogen and silica) is considered. The system is divided into three systems of 
ordinary differential equations, each of which simulates the dependence of phytoplankton growth on a single nutrient. 
These systems are solved by the Runge-Kutta method (Fig. 2—4); stationary modes are obtained, for which the values of 
the system parameters and initial conditions are determined. 

The results obtained will be used for the further simulation of the phytoplankton dynamics considering the 
transformation of phosphorus, nitrogen and silica compounds, taking into account diffusion-convection, saltiness, and 
temperature. 


References 

1. Yakushev, E.V., Sukhinov, A.I., et al. Kompleksnye okeanologicheskie issledovaniya Azovskogo morya v 
28-m reyse nauchno-issledovatel'skogo sudna «Akvanavt». [Comprehensive oceanological studies of the Azov Sea on 
the 28th voyage of the ““Aquanaut” research vessel.] Oceanology, 2003, vol. 43 no. 1, pp. 44—53 (in Russian). 

2. Sukhinov, A.I., et al. Chislennoe modelirovanie ekologicheskogo sostoyaniya Azovskogo morya s prime- 
neniem skhem povyshennogo poryadka tochnosti na mnogoprotsessornoy vychislitel'noy sisteme. [Numerical modeling 
of ecologic situation of the Azov Sea with using schemes of increased order of accuracy on multiprocessor computer 
system.] Computer Research and Modeling, 2016, vol. 8, no. 1, pp. 151—168 (in Russian). 

3. Sukhinov, A.I., Sukhinov A.A. 3D model of diffusion-advection-aggregation suspensions in water basins 
and its parallel realization. Parallel Computational Fluid Dynamics 2004: Multidisciplinary Applications — 2005, pp. 
223-230. DOI: 10.1016/B978-044452024-1/50029-4. 

4. Sukhinov, A.I., Chistyakov, A.E., Shishenya, A.V., Timofeeva, E.F. Mathematical model for calculating 
coastal wave processes. Mathematical Models and Computer Simulations, 2013, vol. 5, no. 2, pp. 122-129. DOI: 
10.1134/S2070048213020087. 

5. Sukhinov, A.I., et al. Predskazatel'noe modelirovanie pribrezhnykh gidrofizicheskikh protsessov na mnogo- 
protsessornoy sisteme s ispol'zovaniem yavnykh skhem. [Predictive modeling of coastal hydrophysical processes in a 
multiprocession system based on explicit schemes.] Mathematical Models and Computer Simulations, 2018, vol. 30, 
no. 3, pp. 83—100 (in Russian). 

6. Sukhinov, A.I., et al. Matematicheskoe modelirovanie usloviy formirovaniya zamorov v melkovodnykh 
vodoemakh na mnogoprotsessornoy vychislitel'noy sisteme. [Mathematical modeling of the formation of suffocation 
conditions in shallow basins using multiprocessor computing systems.] Numerical Methods and Programming, 2013, 
vol. 14, no. 1, pp. 103—112 (in Russian). 

7. Sukhinov, A.I., Sukhinov A.A. Reconstruction of 2001 ecological disaster in the Azov Sea on the basis of 
precise hydrophysics models. Parallel Computational Fluid Dynamics 2004: Multidisciplinary Applications — 2005, 
pp. 231-238. DOI: 10.1016/B978-044452024-1/50030-0. 

8. Debolskaya, E.I., Yakushev, E.V., Sukhinov, A.I. Formation of fish kills and anaerobic conditions in the 
Sea of Azov. Water Resources, 2005, vol. 32, no. 2, pp. 151-162. DOI: 10.1007/s11268-005-0020-5. 

9. Nikitina, A.V., et al. Optimal'noe upravlenie ustoychivym razvitiem pri biologicheskoy reabilitatsi 
Azovskogo moray. [Optimal control of sustainable development in biological rehabilitation of the Azov Sea.| Mathe- 
matical Models and Computer Simulations, 2016, vol. 28, no. 7, pp. 96—106 (in Russian). 

10. Sukhinov, A.I., Nikitina, A.V., Chistyakov, A.E. WU. Matematicheskoe modelirovanie protsessov evtrofikatsii v 
melkovodnykh vodoemakh na mnogoprotsessornoy vychislitel'noy sisteme. [Mathematical modeling of eutrophication pro- 
cesses in shallow waters on multiprocessor computing system.] Bulletin of the South Ural State University: Series “Computer 
Technologies, Automatic Control & Radioelectronics”, 2016, vol. 5, no. 3, pp. 36 —53 (in Russian). 

11. Nikitina, A.V., Puchkin, M.V., Semenov, I.S. Differentsial'no-igrovaya model’ predotvrashcheniya zamo- 
rov v melkovodnykh vodoemakh. [Differential game of fish kill prevention in shallow water bodies.| Large-scale Sys- 
tems Control, 2015, iss. 55, pp. 343-361 (in Russian). 

12. Sukhinov, A.I., Belova, Yu.V. Matematicheskaya model’ transformatsii form fosfora, azota 1 kremniya v 
dvizhushcheysya turbulentnoy vodnoy srede v zadachakh dinamiki planktonnykh populyatsiy. [Mathematical transform 
model of phosphorus, nitrogen and silicon forms in moving turbulent aquatic environment in the problems of plankton 
populations dynamics.] Engineering Journal of Don, 2015, vol. 37, no. 3, pp. 50 (in Russian). 


Mechanics 


1] 


http://vestnik.donstu.ru 


Vestnik of Don State Technical University. 2019. Vol. 19, no. I, pp. 4-12. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk JTonckozo 2zocyOapcméennozo mexnuyeckozo ynueepcumema. 2019. T. 19, Ne 1. C. 4-12. ISSN 1992-5980 eISSN 1992-6006 





13. Degtyareva, E.E., Protsenko, E.A., Chistyakov, A.E. Programmnaya realizatsiya trekhmernoy ma- 
tematicheskoy modeli transporta vzvesi v melkovodnykh akvatoriyakh. [Software implementation of a 3D mathemati- 
cal model of suspension transport in shallow water areas.] Engineering Journal of Don, 2012, vol. 2, no. 4 — 2, 30 p. 
Available at: ivdon.ru/ru/magazine/archive/n4p2y2012/1283 (accessed: 12.12.2018) (in Russian). 

14. Samarskiy, A.A. Teoriya raznostnykh skhem. [Difference scheme theory.] Moscow: Nauka, 1989, 616 p. 
(in Russian). 

15. Sukhinov, A.I., Sidoryakina, V.V., Sukhinov, A.A. Dostatochnye usloviya skhodimosti polozhitel'nykh 
reshenly linearizovannoy dvumernoy zadachi transporta nanosov. [Sufficient conditions for convergence of positive 
solutions to linearized two-dimensional sediment transport problem.] Vestnik of DSTU, 2017, vol. 17, no. 1 (88), pp. 
5—17 (in Russian). 

16. Samarskiy, A.A., Nikolaev, E.S. Metody reshenitya setochnykh uravneniy. [Methods of finite-difference 
equation solution. | Moscow: Nauka, 1978, 532 p. (in Russian). 

17. Marchuk, G.I. Matematicheskoe modelirovanie v probleme okruzhayushchey sredy. [Mathematical model- 
ing in environmental problem.] Moscow: Nauka, 1982, 319 p. (in Russian). 

18. Belova, Yu.V., Chistyakov, A.E., Protsenko, U.A. O chetyrekhsloynoy iteratsionnoy skheme. [On four- 
layer iterative scheme.] Vestnik of DSTU, 2016, vol.16, no. 4 (87), pp. 146—149 (in Russian). 


Received 20.11 .2018 
Submitted 21.12.2018 
Scheduled in the issue 11.01.2019 


Authors: 


Belova, Yuliya V., 

Junior research scholar, Research Institute for Mathemat- 
ical Modeling and Forecasting of Complex Systems, Don 
State Technical University (1, Gagarin sq., Rostov-on- 
Don, 344000, RF), 


ORCID: https://orcid.org/0000-0002-2639-745 1 
yvbelova@yandex.ru 


Atayan, Asya M., 

teaching assistant of the Computer and Automated Sys- 
tems Software Department, Don State Technical Univer- 
sity (1, Gagarin sq., Rostov-on-Don, 344000, RF), 
ORCID: https://orcid.org/0000-0003-4629-1002 
atayan24@yandex.ru 


Chistyakov, Aleksandr E., 

professor of the Computer and Automated Systems 
Software Department, Don State Technical University 
(1, Gagarin sq., Rostov-on-Don, 344000, RF), Dr.Sci. 
(Phys.-Math.), professor, 

ORCID: https://orcid.org/0000-0002-8323-6005 
cheese _O5@mail.ru 


Strazhko, Aleksandr V., 

student of the Computer and Automated Systems Soft- 
ware Department, Don State Technical University 

(1, Gagarin sq., Rostov-on-Don, 344000, RF), 
ORCID: https://orcid.org/0000-0002-2449-853 1 
strajcko2@gmail.com 


Vestnik of Don State Technical University. 2019. Vol. 19, no. 1, pp. 13-23. ISSN 1992-5980 eISSN 1992-6006 
Becmunuk Jlonckozo zocyOapcméennozo mexnuueckozo ynueepcumema. 2019. T. 19, Ne 1. C. 13-23. ISSN 1992-5980 eISSN 1992-6006 





MATTIMHOCTPOEHHE HW MAHTHHOBEJIEHHE 
MACHINE BUILDING AND MACHINE SCIENCE 





UDC 621.382 https://doi.org/10.23947/1992-5980-2019-19-1-13-23 





Study on oil pilot circuit of adaptive hydraulic drive of tool advance in mobile drilling 
machine - 
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Uccaeqopanne ruspaBsMuecKkoro KOHTypa yipaBJIeHHA aanTHBHOTO THApONpHBOAa WO,a4uH HHCTpyMentTa 
MOOMJIbHO OypoBoi MalIMHbI 


B. C. Cugopenko’, B. H. [punrenko’, C. B. Paxystenko’, M. C. Mouemkun’, JI. J. bimouxnn ™ 


1,2,3,4,5 


Introduction. An adaptive hydraulic drive of the tool advance 
in a mobile drilling machine is studied on the example of the 
URB-2.5 installation. A typical technological cycle of the 
mobile drilling machine is considered; the performance 
criteria are defined. An original design of the adaptive 
hydraulic drive is proposed on the basis of the analysis. 
Adaptation of the hydraulic drive of the tool advance is 
carried out using an adjustable volumetric hydraulic motor 
with a hydraulic control circuit under discontinuous loads on 
the tool during the drilling process. 

and Methods. 


computational experiment in the Matlab Simulink program, 


Materials Through a __ preliminary 
the following parameters of the control loop devices were 
determined: a hydromechanical sensor and a hydraulically 
controlled valve, on the basis of which the experimental setup 
was implemented. The performed multifactor experiment 
allowed identifying the processes in the original hydraulic 
control circuit of the hydraulic motor under various modes of 
tool loading. 

Research Results. The kinematic and power characteristics of 
the hydromechanical system of a mobile drilling rig, the 
hydraulic control effect on the settings of the hydraulic 
control circuit devices were obtained and determined. The 
results enabled to specify the rational ranges of the 
hydromechanical system operation for a typical work cycle. 
Discussion and Conclusions. The results obtained can be 
used to create hydraulic systems of new drilling machines 
with various characteristics. The application of the developed 
techniques of research and processing of their results will 
reduce the time and costs involved in designing adaptive 
hydraulic systems for mobile technological machines, 
conducting 


creating prototypes and commissioning 


procedures. 


JloHcKoi rocyapcTBeHHbIM TeXHHU4eCKHH yHUBepcutet, r. Pocros-Ha-Jlony, Poccuitckaa Dezepaina 


Beedenue. CratTba NOcBaleHa UCCIeOBaHUIO ajlalTHBHOrO 
THAponpHBoya MWoOsawi WHCTpyMeHTa MOOWIbHOM OypoBol 
MalliMHbI Ha wWpwMepe ycraHoBKH YPb-2,5. Paccmotpeu 
THMOBOM TeXHOJOrM4YeCKHH IKI MOOMIbHOM OypoBOH 
MalUMHbI, OmpeyieyieHbl KpuTepuu dyHkuMOHMpoBaHua. Ilo 


pe3yIbTaTaM aHasiu3a TIpeWI07%KeHO OpPHTMHasIbHoe 
CX€MOTeEXHHYeckoe pellieHve aflanTHBHOrO TrHuyzpompHBora. 
Ayanrainu THIpPOMpuBoTa noyjaun HHCTpyMeHTa 


OCYIIecTBIAeTCA IIpH MOMOLIM perysupyemoro oOObemMHOrO 
THIpOABHraTesia C KOHTYpOM IpaBlMyeckoro  yrpaBsIeHHA 
TIpH W3MeHAFOMIeicaA Harpy3Ke Ha MHCTpyMeHTe B Iportecce 
OypeHuaA. 

Mamepuazbi u memoooi. IIpeqBapvitesIbHbIM BbIMMCIIMTeJIbHbIM 
9KCIepHMeHTOM B TporpamMe Matlab Simulink ompeyeseHsi 
llapaMeTpbl yCTpOvicTB KOHTypa ylipaBJIeHua: 
THIpOMexaHnyeckoro aTaKa WU THApOyNpaBsAMoro KlaraHa, 
Ha OCHOBe KOTOPbIX peayIM30BaHa 939KCIIePHMeCHTaIbHad 
ycTaHoBKa. BbIMOJIHCHHbIM MHOrOMaKTOpHbIM 3KCIIepHMeHT 
MO3BOIM UAeCHTHUIMpOBaTb MPOLeccbl B OpHIMHasIbHOM 
THIpaBsIMYeCKOM KOHType yiIpaBsIeHHA THIPOMOTOPOM IIpu 
pa3yIM4HBIX pexkKuMax Harpy2KeHuaA MHCTpyMeHTa. 

Pezyibmamel ucciedoeaxus. Tlomry4eHbl UU OlpeyiesieHbl 
KHHeMaTH4eckHe u CHJIOBbIC xapakTepHcTHKU 
THIPOMexaHvyecKOl CHCTeMbI MOOWJIBHOM OypoBOl yCTaHOBKH, 
BIMAHMe THapaBsIMyecKorO§ ylipaBsIieHHa Ha TapaMeTpbl 
HacTpOMKH yCTpOMCTB KOHTypa. Pe3yIbTaTbI MO3BOIIMIM 
OMpeeIHTb pallMOHasIbHble Mala30Hbl (yHKWMOHMpoOBaHHA 
THI pOMexaHv4ecKOH CHCTeMBI JIA THMOBOrO padouero IMKIAa. 
O6cyorcoenue u 3akioueHud. Ilosy4eHHble pe3syJIbTaTbl MOryT 
ObITb HCIOb30BaHbI TPH CoO34aHHH TUAPOCHCTeM HOBBIX 
OypOBbIX MallIMH C pa3IMYHbIMM XapaKTepHCTHKaMH. 
Vcnonp30BaHve pa3pa0oTaHHbIxX MeTOMK UccieyqOBaHuA 
OOpadoTKH UX pe3yIbTaTOB MO3BOUIMT COKpaTUTb 3aTparTbl 
BpeMeHH HU cCpesCTB pH TpoeKTMpoBaHvM = aslalTHBHBIX 


TUTpOCucTeM MOOMJIBHBIX TEXHOJIOTHYCCKHX MaluuvH, CO37aHvu 
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Introduction. Dynamic development of natural resources requires the improvement of existing and the 
creation of new automated complexes of processing equipment for drilling production with improved mechanical and 
energy characteristics. Analysis of the known circuit design solutions of hydromechanical systems (HMS) of mobile 
drilling machines (MDM) identified the key feature of their construction — a multi-engine system, in which it is 
important to consider the effect of the tool feed drive and the main motion drive during the technological cycle [1-4]. 
At this, the quality, productivity and safety of operation depend much on how the kinematic and force parameters (V, 
M, etc.) are matched under changing the tool loads [4, 5]. In view of the above, the work objective was to increase the 
efficiency of the hydromechanical system of MDM working motions through developing and studying its hydraulic 
control circuit. 

Performance Criteria. In the drilling process, an important criterion is the fulfillment of the basic production 
requirement — operation matching of the main motion and feed drive, which would ensure a stable tool advance 
revolution [6]. 

Using the basic laws of similarity of various process cycles, it is easy to apply the above reasoning to the 
technique of drilling various wells [7, 8]. 


The performance (II9) is determined as follows: 


Ho = 1/Tor (1) 

Hoey = Wo ky ky = Lex / Lon (2) 

Ton = Imam + Tec + To3 + Tosca + Trem. (3) 
Tam = Lon /Vai Va = So'Du (4) 

ny = 1000: Vgyp / 3,14:Dy. (5) 


To is reference cycle time; Ty3 is time of adjustment of the next reference operation cycle (increasing the tool 
length); Tam 1S machine time spent on the cutting operation (drilling) when moving the tool by the value of Loy with 
V, Speed. 

When the condition (Sp = const) is fulfilled, the tool life corresponds to Ty = [Tyy] standard with Dy diameter 
and ny frequency of tool rotation. In this case, the number of tool changes throughout Lox = Loyk, drilling depth is 
reduced. 

Tc 1S time spent on installation, commissioning operations, and tool change. It is regulated by the equipment 
type [7]; 

Toscy 18 scheduled time of change (replacement) of the tool after the development of its technological stability 
(Ty); 

Tpem 1S time for eliminating failures, 1t is reduced with the increased reliability; Vsgyp, S, are standards for 


drilling soils of various categories by o; or from the practice of drilling in each region [5]. 
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The tool life and capacity depend directly on the stability of S,(t) value (working process flow per revolution) 
[7, 9]. This is achieved in case, when S,,(t) minute feed rate performed by the drive of the MDM dependent feed, 
decreases synchronously with reducing the tool n(t) rotational speed. When the elastic deformations are not considered 


in the kinematic tool feed chain, then the feed amount is determined as follows: 


Sy(t) = an (6) 
Su(t) = (7) 
w(t) = 2m: n(t) (8) 


where v(t) 1s linear tool advance speed, m/s; w(t) = 27+ n(t) are angular velocities of tool rotation, rad/s; n(t) is tool 
rotation frequency, rev/s. 

As 1s well known, the possibilities of rigid stabilization of each of the working motions under the conditions of 
hydraulic drive application are limited. In the technological machines of this type, the task of adapting a volumetric 
hydraulic drive is not solved automatically, but 1s done manually by the operator [3, 10]. 

Circuit Solution. On the basis of the previously proposed structural scheme [11], as well as the analysis of 
typical MDM operation cycles, the MDM URB-2.5 schematic hybrid (Fig. 1) is proposed. It considers the behavior and 
composition of its mechanical subsystem. 

The machine power system consists of hydraulic fixed axial-piston pumps (H1, H4), installed on the chassis 
transfer gear through power take-offs (PTO). The PTO control is electro-pneumatic and it 1s performed through 
switches installed in the cab of the auto chassis [10, 12]. 

The hydraulic control circuit (HCC) of the installation receives hydraulic energy from H2 two-stage pump 
driven by an electric motor. The electric motor receives power from the auto chassis generator. The first section of H2 
directs the hydraulic energy to GMD1 sensor, the second one — to the hydraulically operated valve (HOV). 

The operation of both HCC circuits occurs off-load, so overheating of H2 two-stage pump motor is excluded. 
Each two-stage pump has a safety relief valve with electrical control (KP2, KP3). Pressure is controlled by MN3 and 
MNS gauges respectively. HCC has its own closed hydraulic tank (B2) [13]. The main hydraulic tank is B1 tank 
equipped with F1, F2, F3 drain filters, TO1 and TO2 heat exchangers. 

The MDM hydraulic system (HS) is divided into two large circuits. The auxiliary circuit includes an outrigger 
drive (TH 1, PH2), a mast lift-lower drive (TI 3), and a winch drive (IM 3). The main contour includes a drive of the 


tool main motion (rotary table) and a feed drive (Fig. 1). 
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Fig. 1. Hydraulic schematic diagram of MDM URB-?.5 


The tool feed drive entering the main circuit receives hydraulic energy from the hydraulic pump (H1). The 
pump is equipped with an electrically-controlled relief valve (KII1). Operating pressure is monitored by the pressure 
gauge (MH2). Reverse movement, stop are provided by a three-position electrically-controlled valve (P2) with a closed 
center. The hydromotor (M1) of the feed drive, adjustable through mechanical transmission, provides a longitudinal 
feed of the drilling tool. The heat exchanger (TO1) and filter (1) provide filtration and conditioning of the working 
fluid. 

The pump (H4) is rotated in the same way as the pump (H1) (Fig. 1). The valve (P3) initiates start, stop of the 
main drive or feed rate increase under the rapid traverse of the feed drive. In this case, the upper position of the valve is 
turned on, and the flow of the pump (H4) enters the drive supply circuit through the check valve (KO1), and is summed 
with the flow of the pump (H1) providing a quick lowering or lifting of the tool from the well. When the pump (H4) 
operates on the main drive (P3 in the lower position), the check valve (KO1) prevents the working fluid from being 
drained out of the supply circuit into the tank (B1) [13]. 
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A similar function is performed by the check valve (KO3) when the pump is operated on the feed drive. The 
valve (P4) provides reverse of the drive main motion. Using the choke (J[P2), the tool rotation speed is adjusted (at 
idle). 

The heat exchanger (TO2) and the drain filter (©3) provide filtration and conditioning of the working fluid of 
the rotational drive circuit. 

The HCC circuit includes two original devices — a hydraulically controlled valve (HOV) (4) and a hydraulic 
multiparameter sensor (HMS) (5) (Fig. 2) [14]. Under the rapid traverse of the feed drive described above, the pump 
(H2) (Fig. 1) 1s automatically shut off by the control system (AIIM). 
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Fig. 2. MDM general view: | is adjustable hydraulic motor; 2 is tool feed drive; 3 is uncontrolled hydraulic motor; 4 is 
HOV valve; 5 is HMS sensor; 6 is gearbox; 7 is driveshaft; 8 is coupling; 9 is auto chassis; 10 is system of outriggers; 11 is tool; 12 
is main motion drive; 13 is double drum winch; 14 is hydraulic cylinder of mast lifting (lowering); 15 is swivel; 16 is mast. 


The HCC principle of operation (Fig. 2) as part of the MDM hydraulic system is as follows: the shaft (HMS) 
(5) 1s mechanically connected to the shaft (TM2) through the gearbox (6). Pressure fluctuation at the HMS input (5) is 
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transmitted to the control input of the HOV (4) which is adjusted — by selecting springs — to the mode of operation 
whereby, at its output, the average operating pressure proportional to the oscillation amplitude and, accordingly, to the 
rotational speed of the hydromotor shaft (M1) is formed. 

With resistance moment increment on the hydromotor shaft (M2), the shaft rotation speed and the oscillation 
frequency at the HMD input decrease, and the pressure amplitude increases. As it increases, the average value of 
opening the HOV valve, which supplies a greater amount of fluid from the H2 to the gearbox (M1) of the feed drive, 
increases. The control pressure (Py) at the HOV input increases, which provides an increase in the hydraulic motor 
capacity and reduces the feed speed. The mechanism of the HCC contour as a component of the HMS with a dependent 
tool advance, as well as the operation principle of the (AIIM) unit, is considered in [15]. 

The accumulator (AK) smoothes the pressure control pulsation (Py) at the hydraulic motor (M1) input 
control. Using the setting of the chokes (JIP1, J[P2) and the selection of the pressure valve (KO3), the limiting values of 
the tool feed rate and rotation speed are adjusted, and the required functional relationship between the working motion 
speeds is formed. 

The original HOV differs in function of continuous regulation of the valve flow section. Therefore, it is 
possible to control the performance of the hydromotor at all speeds. 

Experimental Study. To identify the HMS parameters operating under the specified conditions, special bench 
equipment and accessories are developed [13]. The hydromotor flow rate was determined using a turbine sensor-flow 
meter connected through the converter board (IJAITI-ALIIT). The adjustment ranges of the parameters under the study are 
shown in Table 1. 

Table 1 
Parameter ranges under HMS study 


Unit of 
measurement 


0513 Phos oa 
Choke flange | dy ——*«? O82 | mm ~—~«YSS*~h ga 


Gap between nozzle and 
eee en cn i cc 


Fl 
gage tank 


Pressure in HWS MPa 
Characteristic design factor 
2-8 Slid 


The traverse speed and acceleration of the shaft (IM) (feed drive) was determined by sequential differentiation 


sensor 





of motion over time using the following formulas implemented by post-processing of data in the PowerGraph program 
[8, 11, 16]: 





ot) =f (9) 
c(t) = =. (10) 


where w is angular velocity, rad/s; ¢ is angular acceleration of (IM), rad/s. 


The results of the obtained experimental data processing through the known methods [17, 18] are presented in Fig. 
3-9, 


The graph in Fig. 3 explains the pressure fluctuation amplitude response of the HMS under a change in the 
hydromotor (M1) speed variation within the range of 45—125 rad/s. The data were obtained when testing nozzles with 
the diameters: d,,, = 2, 4, 6 mm. 


As a result, the dependence of the pressure change was almost linear in nature, and it fell as the speed 
increased, which was associated with the non-stationary mode of the working fluid outflow through the flapper-nozzle 
unit. 
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Fig.3. Dependence of pressure fluctuation value (P7;) on @ speed of rotation of hydromotor (M1), approximation: | is 
P=0.140°-3+24.8; 2 is Py=0.36m*-3.84@+20.4; 3 is Pyp=0.210°-2.80+16.6. 


The results presented in Fig. 4 show that with an increase in the characteristic structural parameter of the flow 
section (TMJ Xq) from 2 to 8 mm, the control pressure level changes from 1.9 to 1.4 MPa. This enables to make 
recommendations for further optimization of the HMS flow path [19], in particular, for some increase in (X/JI). 
However, its further increase 1s impractical because it increases the geometric dimensions of the modulator disk. 
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Fig. 4. Dependence control pressure (P,) variation on effect of characteristic structural parameter of cross section of 
flow part (TMJ Xx), approximation: Py=1.65X)+20.75 


Fig. 5 shows the effect of the design features of (X,) on the maximum pressure amplitude (Pymax) used for the 
hydraulic control circuit during advance to the HOV. 
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Fig.5. Dependence of change in maximum amplitude of control pressure (Py ax) on effect of characteristic design parameter of flow 
section of ((MJI Xq), approximation: Pyyax=21.46X °°” 
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During identification, an important step was the determination in the HCC circuit: HMS — HOV, the degree of 
impact of each control component setting. Thus, the HOV tuning property, under the control loop performance, 1s its 
spring constant (C,,), the selection of which significantly affects the sensitivity of the control subsystem [20]. 

The characteristic built in Fig. 6 and approximated by the results of the experiment, is linear, and explains the 


magnitude of the maximum control pressure drop from 1.7 to 1.4 MPa under changing the HOV spring force. 
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Fig. 6. Dependence of maximum control pressure (Py max) variation on magnitude of HOV spring force, approximation: 
Pymax=1 Cpt] 8.5 


The second tuning element in the circuit with the HOV is an adjustable choke in the shunt line (Fig. 1). Its 
setting parameter is the flow area regulated within the range from 2.5 * 10°” to 7.5 * 10° m’. The results of the 
experiment are shown in Fig. 7. 


P,, Max, 
MPa 























Oo rF NY WwW fF UI HD N 


25 3015 a) 6.25 


75 £107 m7? 
Fig. 7. Dependence of maximum control pressure (Py max) variation on flow area of choke (J[P1), approximation: 
PymMax=1 Sfyp°-27.6 : fypt 84.2 
Four operation modes of the tool-drive motor (hydraulic motor (IM 1)) under relay loading, and, accordingly, 
under changing @ angular velocity by: 22 rad/s; 49 rad/s; 71 rad/s and 85 rad/s, are studied. The dependence obtained 


from the experiment is shown in Fig. 8. 
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Fig. 8. t, dependence of HCC response time on @ angular rotation velocity of hydromotor (I'M1) shaft, approximation: 
ty=-0.0027-0.4@+0.83 
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Since HCC 1s a part of the MDM hydraulic system, its length has a special impact on the quality and time of 
transient processes [16]. In order to determine the degree of impact on the above parameters, the pressure variation time 
in the control line with varying of its volume was studied (Fig. 9). 
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Fig. 9. t, dependence of HCC response time on W volume in control flow line, approximation: 
ty=-0.02@°+0.1@7+0.12@+0.5 
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The results enable to determine the critical value of the flow line volume: W = 24-10% m’, over which the HCC 
response time does not exceed 1.5 s; that is due to wave processes taking place in the pipeline, as well as to the 
parameters of the lines themselves [20]. 


Conclusions. As a result of the research, the requirements were developed, and the HMS generalized structure 
of the MDM working motions was proposed, which made it possible to increase the efficiency of drilling production 
through introducing an internal kinematic hydromechanical connection between the main motion and the tool feed 
movement. 

A hydraulic control circuit is developed and implemented on the basis of a multi-parameter sensor and a 
hydraulically controlled valve, which ensures the coordination of working motions under the discontinuous process 
duty. Considers the characteristic features of the drilling process, the configuration of the reference operating cycle, 
which allowed for the development of an automated system for driving work motions, is proposed. 
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Introduction. The possibility of using monocarboxylic acids as 
a lubricant composition additive, and the effect of their 
concentration in lubricant on the evolution of the friction 
factor of a brass-steel couple, as well as the morphology of the 
film surface under friction is considered. The work objective is 
to study the effect of the concentration of carboxylic acids in 
the lubricant composition on the evolution of the friction 
factor of copper — steel alloy. 

Materials and Methods. Tribological studies of a brass-steel 
friction couple in aqueous solutions of monocarboxylic acids 
with the concentrations of 0.025; 0.05; 0.1; 0.2; 0.5 mol/l are 
carried out. Using scanning electron microscopy, we have 
studied the morphology of the servovite film surface that is 
formed on a steel disk after frictional interaction of a brass- 
steel couple in aqueous solutions of acids with the 
concentration of 0.1 mol/l. 

Research Results. Tribological characteristics of the brass- 
steel tribocoupling in aqueous solutions of carboxylic acids of 
various concentrations are studied. The optimum acid 
concentration in the lubricant composition is_ specified. 
Herewith, a selective transfer and a wearless friction regime 
are implemented under friction of the brass 59-steel 40X 
couple. A decrease in the friction ratio to 0.009 and 0.007 is 
found out under friction in aqueous solutions of valeric and 
caproic acids, respectively. The formation of an anti-friction 
film on the steel surface is identified through the scanning 
electron microscopy. It is established that the film formed in 


an aqueous solution of caproic acid has a denser structure in 
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Beedenue. CrTaTba TOCBAIIeHa UCCIeAOBaHHIO BO3MO2%KHOCTH 
HCHOJIb30BaHHA B KayecTBe pwcaaqKw K CMa304HOM 
KOMIO3HIIMH OJHOOCHOBHBIX KapOOHOBBIX KHCIOT 4H 
W3yYYeHHIO UX BIIMAHHA Ha IBOJIOMMIO KOIPMuUMeHTa TPeHHA 
llapbl JlaTyHb-cTalib, a TakoKe W3y4eHHIO MOpdouOrMUu 
(DOpMUPYFOIUXCcA IPH TPCHHM MOBeEpXHOCTH MJICHOK. 

Ilembro paOoTHI ABIIAIOCb W3y4CHHe BIMAHHA KOHUeCHTpaluu 
KapOOHOBbIX KHCJIOT B COCTaBe CMa304HOM KOMIMO3HUMM Ha 
S9BOJFOMMIO KOIPDUMUMeHTa TPeCHUA Mapbl CiMIaB MeJ{b-CTaJIb. 
Mamepuaibi u memodoi. Tpopefenbl tTpHOonormueckue 
WccIeOBaHHA Wapbl TPeHHA JlaTyHb-CTaJIb B  BOJIHBIX 
KapOOHOBBbIX 


pacTBopax KMCIIOT Cc 


KOHUeHTpauMaMu 0,025; 0,05; 0,1; 0,2; 0,5 Momnb/n. C 


OTHOOCHOBHBIX 


TOMOIIbIO PaCTPOBOM DICKTPOHHOM MUKpOcCKONIMU U3y4eHa 


Mopdouorna MOBepXHOCTH CepBOBHTHOM IIICHKH, 
(opMupyrollelica Ha CTasIbHOM JMCKe Tocue dpuKUMOHHOTO 
B3aMMOJeCHCTBUA Wapbl TpeHHA JIaTyHb-CTasIb B BOJIHBIX 
pacTBopax KHCJIOT C KOHUeHTpaunert 0,1 Mosb/I. 

Pe3yibmamoi ucciedoeanua. VU3y4eHbl tTpuOonormueckue 
XapaKTepHCTHKH TPHOOCONps.KeHHA JIATYHb-CTalIb B BOJIHbIX 
pacTBopax KapOOHOBbIX KHCJIOT pa3JIM4YHOM KOHIeHTpallunu. 
YctTaHOBJIeHa ONTHMaJIbHad KOHI[CHTpalluA KUCOTbI B 
COCTaBe CMa3KH, IPH KOTOPOM B pe3yiIbTaTe TpeHuA Mapbi 
maTyHb 59-ctanb 40X peamn3yertca u30upaTesIbHbIM MepeHoc 
WM JOCTHTaeTCA pexUM Oe3bI3HOCHOTO TpeHud. OOHapy2xxKeHO 
CHWKeHHe KOIddulUneHTa TpeHuA TO 0,009 u 0,007 npu 
TPpeHHH B BOJIHbIX pacTBOpax BaJIepHaHOBONM HU KallpOHOBON 
KMCJIOT COOTBETCTBCHHO. C  MOMOIIbIO pacTpoBor 
39JIEKTPOHHOM MMKPOCKOMMU BbIABIICHO (OpMHpoBaHHe Ha 
CTaJIbHOU aHTHPPUKIUMOHHOK ~=—s WICH KH. 
YctTaHOBJIeHO, 4TO TJIeHKa, (OpMUpyrolllaicx# B BOTHOM 


pacTBope KallpOHOBOM KHCIIOTbI, WMeeT OOsIee TIOTHYFO 


IIOBeEpXHOCTH 


CTpyKTypy, B CpaBHeHHu C MWIeHKOH, POopMupyrolleiica mpu 
TpeHHH B BOJHBIX pacTBOpax MacJIdHOM U KallpOHOBOU 
KMCIIOT. 
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comparison with the film formed under friction in aqueous 
solutions of butyric and caproic acids. 

Discussion and Conclusions. Thus, the tribological studies of 
a brass-steel friction couple in aqueous acid solutions show 
that the optimum molar acid concentration in the lubricant 
composition is 0.1 mol/l. At this acid concentration, the values 
of the friction factor characteristic of the wearless mode are 


attained. 


Keywords: friction factor, wear, selective transfer, servovite 
film, carboxylic acid, friction surface topography. 
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O6cyacdenue u 3ak0Y4eHUA. B pe3ylbTaTe ToHOouOrMyecKHXx 
UccIeqOBaHHH Wapbl TpeHHA JlaTyHb-CTaJib B  BOJ[HBIX 
pacTBOpax KUHCJIOT BbIABJICHO, UTO ONTHMasIbHOH MOJIAPHOH 
KOHIeCHTpalel KUCJIOTbI B COCTaBe CMa3KH ABJIAeCTCH 
KOHLeHTpauua 0,1 Mosb/n. [Ipu stow KOHWeHTpallwu KUCIOTHI 
WOCTHraroTcs 3HAaYeHHA KOOMDUUMeHTa TPeCHHA, XapaKTepHBIe 
JId pexuMa Oe3bI3HOCHOCTH. BbIABSIeHO, YTO U3MeHeHHe 
KOHICHTpalluH KUCIIOTbI MpwBOAMT HOO K yBeJINYeHHIO 
3HayeHHi KOIMUUMeHTa TpeHuIO, JIMO0 K He3HaYMTeJIBHOMY 
ero cHwxKeHHIO. IIpu 9TOM yBeM4eHHe KOHIeHTpallun 
KUMCJIOTbI COMPOBOXKTaeTcA KOPpO3HOHHbIMH Ipoljeccamu Ha 
MOBePpXHOCTH Tapbl TPeHHA. 


KoodulWeHT 
W30uUpaTeIbHbIN MepeHOC, CepBOBUTHAaA MWIeHKa, KapOOHOBad 
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Introduction. Modern high-developing machine-building industry advances new demands on lubricants. Most 
conventional methods of the friction and wear control are based on the use of solid and liquid lubricants [1—4]. The 
principal function of lubrication in friction units is to keep two contacting surfaces of machine parts from wear. To 
perform the necessary functions, basic fluids need modification by functional additives that change the antiwear 
properties of base oils through improving their lubricity among other. For this reason, additives are an integral part of 
the design of modern lubricants. The studies [2—7] show that to reduce friction and wear, both nanoparticles of metals 
and various organic components are used as additives to lubricant compositions. Some of them contribute to the 
formation of protective antifriction films on the tribocontact surfaces due to the presence of metal powders with particle 
sizes in the micro-range and nanoscale in the lubricant composition [7—9]; others — as a result of selective dissolution of 
tribo-conjugated surfaces (in the case of copper-steel alloy friction pair) under friction [10-17]. As is known, the 
lubricating medium plays a critical role in the selective transfer of copper to the steel surface under frictional 
interaction. For example, during the copper alloy — steel friction, various polar compounds, including carboxylic acids, 
are formed in the aqueous-glycerin medium. It was interesting to study the possibility of using them as additives to the 
lubricant composition for implementing a selective transfer and wearless friction. 

Materials and Methods. The evolution of the friction factor of the “brass 59-aqueous solution of carboxylic 
acid — 40X steel” system was investigated on the AE-5 type face friction machine. Aqueous solutions of saturated 
monobasic carboxylic acids with the general formula: R-COOH (R = C,Ao,+)), with the concentrations of 0.025—0.5 
mol/l were used as a lubricant composition. Before the tribological studies, samples of steel 40X and brass 59 were 
cleaned with abrasive paper, washed with distilled water, degreased with hexane, and air-dried. The friction assembly 
consisted of a steel rigidly fixed disk and three movable brass fingers arranged in a circle at the angle of 120° relative 
to each other. A friction couple of a ring steel specimen and brass fingers was placed in the lubricant composition into 
the working part of the friction machine made in the form of a textolite bath. A PHYWE Cobra force sensor was 
attached to the front face of the working part of the friction machine for continuous recording of the friction variation. 
Tribological studies were carried out under the following modes: sliding speed of the specimen was 0.45 m/s; axial 
load was 1.7 MPa; test time was 10 hours; temperature of the working environment was 37° C; slide path was 15,260 
m. 

The morphology of the servo film was studied using the scanning electron microscopy (SEM) on a Carl ZEISS 
microscope in the Electron and Optical Microscopy Laboratory and the Center for Collective Use of Scientific 
Equipment, REC “Materials” (http: //nano.donstu.ru). The studies were conducted under high vacuum conditions. 
Accelerating voltage under the scanning mode was 1-3 kV. 

Research Results. The dependence of the friction factor on the concentration of carboxylic acid for various 
volume acid densities under the brass-steel couple friction is given in Fig. 1 and in Table 1. From the analysis of the 
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data obtained, it follows that the dependence in Fig. 1 1s nonmonotonic, with a minimum concentration of acid in the 


solution of 0.1 mol/I. 
Table 1 


Dependence of friction factor on acid concentration in lubricant 


Lubricant Acid concentration, mol/I 
paecisnuacd Acid chemical formula: R- 
composition, acid | Cr where R | 0033 


solution friction factor, u 


-CH,CH; 0.195 0.129 0.087 0.239 0.335 





-CH,CH.CH; 0.112 0.083 0.037 0.184 0.307 
-CH»CH,CH,CH; 0.071 0.031 0.009 0.045 0.154 
-CH,CH>CH,CH>CH; 0.054 0.016 0.007 0.028 0.126 


Studying the effect of acid concentration on the friction factor variation in the solution shows that strengthening of 
formic and acetic acids in the lubricant composition from 0.025 mol/l to 0.05 mol/l leads to the increase in the acidity of 
the medium and, consequently, to an increase of the values of the friction factor from 0.3 to 0.35. 

As a result of frictional interaction, tribocorrosion occurs on the metal surface in the tribosystem. A further 
increase in the acid concentration to 0.1 mol/l causes a decrease in the friction factor, and then it is accompanied by its 
sharp increase and strong wear of friction couple materials as a result of the adhesion-mechanical interaction of the 








protrusions of the surface microrelief (Fig. 1, 2). 
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Fig. 1. Dependence of friction factor (1) on acid concentration (C) in lubricant composition under friction of brass-steel couple: | is 
formic acid, 2 is acetic acid, 3 is propionic acid, 4 is butyric acid, 5 is valeric acid, 6 is caproic acid 
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Fig. 2. Evolution of friction factor (uw) on concentration (C) in “brass-aqueous solution of acetic acid-steel” system 
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A similar dependence of the friction factor on the acid concentration in the lubricant composition is also 
observed under friction of a brass-steel couple in aqueous solutions of propionic and butyric acids (Fig. 3). 
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Fig. 3. Evolution of friction factor (uw) on concentration (C) in “brass-aqueous solution of butyric acid-steel” system 


At this, the friction factor decreases significantly at the acid concentration of 0.1 mol/l and does not exceed 0.1. 
A further increase in the acid concentration in the lubricant composition, as well as in the case of formic and acetic 
acids, is accompanied by its sharp increase (Fig. 1). As a result, wear products are formed in the lubricating fluid 
volume. In this case, the friction surface is subjected to corrosion-mechanical wear, and the friction factor has rather 
high values from 0.15 to 0.35 (Fig. 3). 

Analysis of the friction factor variation of a brass-steel couple in aqueous solutions of valeric and caproic acids 
with the concentrations of 0.025 and 0.05 mol/l determines rather low value to 0.07 (Fig. 1, 4). 

Application of valeric and caproic acids in a lubricant composition with the concentration of 0.1 mol/l enables 
to obtain the lowest values of the friction factor to 0.007; optimal conditions for self-organization on the steel surface; 
the formation of a visually detectable servovite film [13, 18]; and the establishment of a wear-free mode in the 
tribological system. At this, reduction of the friction factor in the copper-steel alloy couple is associated with the best 
damping of tribo-film stresses caused by friction. 
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Fig. 4. Evolution of friction factor (uw) on concentration (C) in “brass-aqueous solution of caproic acid-steel” system 


Triboelectrochemical reactions occurring in the frictional contact zone in aqueous solutions of valeric and 
caproic acids not only initiate the formation of a servovite film, which helps to reduce the friction factor, to heal surface 
microdefects, but intensify the chemosorption which enhances the ordering effect of the substrate on the orientation of 
molecules under the formation of adsorption acid molecule layer. 

Active —COOH polar groups in a carboxylic acid molecule lead to its interaction with metal surfaces with the 
formation of chemically adsorbed compounds of the bidentate-ligand type [12]. Aliphatic monocarboxylic acids (R — 
COOF) form layers in which the hydrocarbon radicals of the molecules form a closely-packed structure. The adsorbed 
carboxylic acid molecules on the metal surface not only reduce the friction factor, as compared to a clean friction 
surface, as it follows from the results obtained, but also increase it. At very low acid concentrations in the solution — to 
0.05 mol/l, the degree of filling the adsorption layer is very small; the molecules do not form a continuous film and 
move freely along the metal surface. Aside from the acid molecules, water molecules can adsorb on the metal surface; 
however, they do not have an effective shade action, therefore, the friction factor is almost the same as for clean 
surfaces. As the acid concentration in the solution rises to 0.1 mol/l, the occupancy of the adsorbed friction surface layer 
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increases, while the acid molecules run parallel to the surface reducing the friction factor as compared to the values for 
clean friction surfaces. The generated boundary layers at low acid concentrations reduce the friction factor. With further 
raise of the acid concentration, molecules are arranged perpendicular to the surface; the adsorbed layer increases the 
surface roughness; and the friction factor grows. The closedness of intermolecular bindings inside the adsorption layer 
leads to low surface energy and virtual absence of secondary adsorption of acids or other components [19-21]. 

Analysis of the results obtained through the scanning electron microscopy (SEM) indicates drastic structural 
changes in the friction surface in acid aqueous solutions when moving from formic to caproic acid (Fig. 4). As can be 
seen from the SEM visualization of the results obtained, with the relative movement of two surfaces in an aqueous 
solution of acetic and butyric acids, the surface of the tribocontact has a heterogeneous structure with a large number of 
pores and irregularities resulting from mechanochemical corrosion; as well as the occurrence of wear particles on the 
friction surface which causes abrasive wear of rubbing metals. At this, cracks and scratches appear on the friction 
surface, which significantly reduces wear resistance of the material. Micrographs of the sample surface after friction in 
aqueous solutions of acetic and butyric acids show major roughness of the friction track, which leads to an increase in 
the friction factor (Fig. 4a, b). The friction force growth also provokes an increase in temperature in the friction unit, 
and, as a result, thermal stress of the metal, which is also one of the causes of the crack formation, both on the material 


surface and in its volume. 





Fig. 5. SEM results of servovite film obtained under friction in “brass-aqueous solution of acid-steel” system: 
(a) is acetic acid, (b) is butyric acid, (c) is caproic acid 


When an aqueous solution of caproic acid is used as a lubricant composition under friction of a copper alloy 
over steel, the formation of a copper non-oxidizing layer on the surface of a steel disk is observed (Fig. 4c) [22]. This 
effect 1s explained by the catalytic effect of copper, which converts the lubricant monomers into polymers, ensuring that 
copper remains in an unoxidized state [12]. Copper film, which is formed through friction in an aqueous solution of 
caproic acid, has a dense structure with a minimum number of pores, ensuring wearless friction. 


Conclusion. As a result of the tribological studies of a brass-steel friction couple in acid aqueous solutions, the 
optimum acid concentration in the lubricant composition was found. It is established that a change in the acid 
concentration causes a change in its adsorption on the friction surface, leading to a change in the morphology of the 
tribo-conjugated surfaces. 
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Introduction. The electrodeposition of composite electrochem- 
ical coatings from electrolyte-colloid nickel plating containing 
ultradisperse zirconium boride powder is studied. The work 
objectives are as follows: to study mechanical-and-physical 
properties of the composites based on nickel and nanostruc- 
tured zirconium boride, and to determine optimal conditions 
for the application of such electrochemical coatings. 

Materials and Methods. Microhardness of composite electro- 
chemical coatings was measured using PMT-3 microhardness 
tester on samples with the layer thickness of 30 um under the 
indentation load of 100 g. A three-ball machine was used to 
determine wear resistance of the coatings. Sample tests were 
carried out under dry friction modes and with the use of 3% 
RV coolant. WSD values were measured by MIR-3 TU 3- 
3.1954-86 microscope. To determine the internal stresses in 
the coating, we used a flexible cathode method up to GOST 
9.302-88. 

Research Results. The electrolyte-colloid composition and 
modes of electrodeposition of composite nickel - nanostruc- 
tured zirconium boride coatings are developed. Mechanical- 
and-physical properties (microhardness, wear resistance and 
internal stresses) of the obtained composite electrochemical 
coatings are analyzed. Recommendations for use of the devel- 
oped electrolyte and the application of a composite coating on 
machine parts for their surface hardening are formulated. 


Discussion and Conclusions. Ni-ZrB, CEC (composite elec- 
trochemical coating) has high microhardness (10—11 hPa at the 
indentation load of 100 g), which exceeds the microhardness 
of pure nickel by 1.5—2 times. As the microhardness increases, 
the internal stresses of Ni-ZrB, CEC decrease. The proposed 
coatings were compared to chromium ones deposited from the 
environmentally hazardous electrolytes. The wear resistance of 


” The research is done within the frame of the independent R&D. 
“E-mail: degtiar@yandex.ru, i.ivanina96@mail.ru, iyuzh@mail.ru 
““Padota BBINOJIHeHa B paMkax HHuWMaTuBHOK HUMP. 


Beedenue. Uccneqopan mpotecc d1eKTpoocaxeHuA KOMMO- 
3HTHBIX 93JICEKTPOXHMHYeCCKHX TOKPbITHHM 3 3ICKTPOJIMTA- 
KOJWION a HUKeIMpOBaHHA, COJepxKalllero yIbTpayMciepcHbIM 
MOpoOwoK WHOopuya WMpKOHHA. 

Ilemu padoTsi: uccueqoBaHve (U3HKO-MexaHHdeCcKHX CBOMCTB 
KOMIIO3HTOB Ha OCHOBe HUKeJIA HM HAHOCTPyKTYpHoro WHOopu- 
fla IMpKOHUA, a TakKx%Ke OlIpeyjeseHve OMTHMAJIbHBIX YCIOBMU 
HaHeCeHHA TaKOFO IICKTPOXHMHYECKOLO MOKPbITHA. 
Mamepuanot u memodoi. MukpoTBepaocTb KOMIMO3HTHBIX 
3JIEKTPOXHMHYeCCKHX MOKPbITHM U3MepsJIM C MOMOIIbIO MHK- 
potBepzomepa IIMT-3 Ha oOpa3ilax c TONMMHOK cu04 30 MKM 
IIpH Harpy3kKe Ha uHeHTOp 100 r. Ja onpezemreHua U3HOCO- 
CTOMKOCTH TIOKPBITHM UCMOb30BaJIM TpexIWapHKOBYy!O MallIu- 
Hy TpeHua. Ucnbiranua oOpa3ill0B MpoBOAWIM B pexuMax 
cyxoro TpeHua UW Cc UpwMeHeHvemM 3% CMa304HO- 
oxaxKyaroljenh %#*uTKOcCTH PB. 3HayeHHaA MaMeTpa MATHA 
W3HOCa U3MepHIM To, MUKpocKoTIOM MUP-3 TY 3-3.1954-86. 
JIna ompeyeweHud BHYTpeHHUX HallpsxKeHHM B TOKpbITHU 
BOCIIOJIb30BaJIUCb MCTOJIOM TMOKOrO KaTOJa B COOTBETCTBHH C 
TOCT 9.302-88. 

Pe3yiemamoi ucciedoeanus. Pa3pa0oTaH COCTaB 93JIEKTPOIIN- 
Ta-KOJWIOWa WU PexXUMBI JICKTPOOCAaKCHHA KOMIMO3HTHbIX 
MOKPbITHM HUKesIb — HAHOCTPYKTypHbIM WHOOP WHpKOHHA. 
IIpopeyqeH aHamu3 Pu3HKO-MexaHH4ecKHx CBOMCTB (MUKpo- 
TBepOCTb, H3HOCOCTOMKOCTb HM BHYTpeHHve Halips2KeHus) 
MOJTYYCHHBIX KOMIMO3HTHBIX 3JICEKTPOXHMNYeCKHX TOKPbITHH. 
CcopMysIMpoBaHbI PeCKOMeHallMN MO UCIOIb30BaHHIO pa3pa- 
OOTAHHOrO IICKTPOJIMTa HW HAHECeCHUIO KOMIIO3HTHOrO MOKDpBI- 
THA Ha JeTaJIM MallIHH JIA UX MOBEPXHOCTHOTO yIIpoOuHeHH aA. 
O6cyotcdeHue u 3aKouenue. KOII Ni-ZrB, uMeerT BbICOKy!O 
MukpoTBepyoctTb (10-11 rila mpu Harpy3ke Ha WHAeHTOp 
100 r), 4TO MpeBbIaeT MUKPOTBeEpAOCTb YHCTOO HUKeIIA B 
1,5—2 pasa. [pu Bo3pacraHuu MUKpOTBepOCTH CHW KaIOTCA 
BHYTpeHHue HalpsKeHnA KOI Ni—ZrB,. IIpennaraempie mo- 
KPbITHA CpaBHUBaJIMCb C XPOMOBBbIMH, OCaKJaCMbIMU V3 IKO- 
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Ni-ZrB, CEC is 2—5 times higher than that of chromium coat- 
ings. Thus, instead of chromic coatings, it is recommended to 
use the proposed composition for surface hardening of parts of 
the specialty machinery and industrial equipment. 


Keywords: electrolyte, composite electrochemical coating, 
nickel, zircontum boride, microhardness, wear resistance, 
internal stresses. 
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JIOrH4YeCKH OMACHbIX BIeKTpoMTOB. M3HococTorwKoctTb KIIT 
Ni-ZrB, B 2—5 pa3 Oombille, 4eM Y XPOMOBBbIX NOKpbITHH. Ta- 
KMM 00pa30M, PeKOMeHyeTCA BMeCCTO XPOMOBBIX NOKPbITHH 
MCHOIb30BaTb TpeiaraeMbIM COCTaB JIA MOBepXHOCTHOrO 
yYIIpOUHeHHA WeTasei ciellMasIbHOW TeXHUKU VU TpOMBILLIJIeH- 
HOro OOopyOBaHHaA. 


Ks1104eBble CIOBA: SICKTPOJIMT, KOMIIO3HTHOE 3JIEKTPOXHMH- 
yecKOe NOKPbITHe, HUKesIb, WHOOpuA WHpKOHHA, MUKPOTBep- 
JJOCTb, H3HOCOCTOHMKOCTb, BHYTpeHHve Halips.KeHHA. 


Oodpa3zeu Ona YUmuposanua: OcodeHHOcTH PopmMupoBaHHA 
KOMIO3HTHBIX 9JICKTPOXHMHYCCKHX TOKPbITHH Ha OCHOBE 
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Introduction. Composite materials are widely used to strengthen the surface of typical nodes and parts of the 
petroleum-gas field, processing and other industrial equipment. As a result of a combination (composition) of dissimilar 
substances, a new material appears whose characteristics differ from the properties of its components [1]. 

To improve the performance properties of electroplating, particulate fillers are introduced into the solution. 
Thus, a composite electrochemical plating (CEP) with a matrix of the base metal is obtained. Precisely the creation of 
the CEP is one of the current areas of modern electroplating. To obtain the CEP, dispersed particles of various sizes and 
types from electrolytes codeposit simultaneously with metals. These particles are included in the coating; they improve 
drastically the operational properties (hardness, wear resistance, corrosion resistance) of the coatings, and give them 
new qualities (anti-friction, magnetic, catalytic). 

Various technologies of producing nickel-based CEP [2—5] are known. The selection of technique is deter- 
mined by the operating conditions of the composite coating (CC). At this, wealth of fillers for obtaining coatings with 
specific performance properties complicates the choice of dispersed particles [3, 5—9]. In [2, 3, 6—8, 10], the physicome- 
chanical properties of various nickel-based CC with particulate fillers (for example, nanocarbon material taunite, alumi- 
num oxide, silicon carbide, graphite, etc.) are considered. In the mentioned works, in particular, it was noted what de- 
termines the structural and morphological characteristics of the CEP [7—9], their resistance to corrosion [7, 11], the co- 
efficient of sliding friction [8], heat resistance, microhardness, and other properties [11] . Conditions and modes of CC 
deposition, composition of electrolytes, nature and degree of fineness of powder or other additive materials are men- 
tioned as determining factors. The use of micropowders with a grain size of more than 1 um and ultrafine powders or 
suspensions can change the nature of the codeposition of particles with the excreted metal, and, therefore, the character- 
istics of the CEP [9]. 

The use of micropowders with a grain size of more than 1 um and ultra-disperse powders or suspensions can 
change the nature of the codeposition of particles with the excreted metal, and, consequently, the characteristics of the 
CEP [9]. During the dispersion phase, the properties of electrolytes and coatings from them gain a number of ad- 
vantages [12, 13, 14] critical from the point of view of the applicability. 

Nickel CC obtained from electrolytes with additives based on zirconium dioxide to improve surface morpholo- 
gy, microhardness, heat resistance and other properties [5, 15] are of interest. Zirconium boride (ZrB,) as an additive to 
electrolytes for obtaining CC was not studied. However, it is known that zirconium and boron compounds improve 
drastically the mechanical and processing properties of electroplating [16]. . This work objective is to study the physi- 
comechanical properties of the CEP based on nickel and nanostructured ZrB,, as well as to determine the optimal condi- 
tions for its deposition. 

Materials and Methods. Reagents of Aldrich Company and distilled water were used for the fabrication of 
electrolytes. The content of ZrB, in the coating was determined by the weight (gravimetric) method. The CEP micro- 
hardness was measured with PMT-3 microhardness tester. For this purpose, samples with the layer thickness of 30 um 
were used under the indentation load of 100 g. 

The current output of the coatings was obtained using a copper coulometer. The wear resistance of the coatings 
was determined on a three-ball friction machine. Samples were tested under dry friction modes and with the use of 3% 
RV coolant [17]. The wear resistance of the CEP was investigated as follows. The balls of steel ShKh 15 with the area 
of 0.05 dm’ were covered with 30 um thickness coating. Washers of St 45 steel served as a counterbody. The wear scat 
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diameter was measured using MIR-3 TU3-3.1954-86 microscope. To determine the internal stresses in the coating, we 
used the flexible cathode method in accordance with GOST 9.302-88. 

Research Results. The electroplating process and some features of the CEC with ultradisperse powder (UDP) 
ZrB, (0.04—0.06 um) are studied. The creation of nickel-based CEC with certain physical and mechanical characteristics 
requires consideration of several factors. In this case, first and foremost, the challenge of the interaction of the dispersed 
phase, introduced earlier into the electrolyte, and colloid particles formed in the solution under the electrolysis or in the 
preparation of electrolyte is considered [12, 13]. 

For the studies, sulfate-chloride nickel plating electrolyte-colloid of the following composition was used, g/l: 
nickel sulfate — 250, nickel chloride — 60, a-aminoacetic acid — 20. ZrB, concentration varied from 1 to 60 g/l. 

The optimal electrolysis conditions in the presence of ZrB, UDP are specified. For this purpose, ranges of the 
cathode current density and the electrolyte pH, in which the CEC of good quality is formed, are determined. In the elec- 
trolyte composition under study, a-aminoacetic acid was present as an additive agent. Its maximum buffer properties 
were shown at pH 2-2.5: in this case, light matte coatings of good quality were formed. The presence of an effective 
additive agent in the electrolyte allows for electrolysis at the high cathode current density. 

If ZrB>, 1s absent in the electrolyte, then the limiting value of the cathode current density is 2 A/ dm*. When this 
value is exceeded, a dark contour is formed on the test samples. This indicates alkalization in the cathode layer, the 
formation of nickel hydroxide in a coarsely dispersed form, and its inclusion in the coating [12]. In this case, the physi- 
comechanical properties of the precipitation deteriorate, their fragility increases, and the coating cracks along the edges 
of the cathode. 

The introduction of ZrB, UDP in the electrolyte affects drastically the limiting value of the current density. 
Good precipitates are formed above the permissible current density characteristic of the electrolyte under study without 
adding ZrB>. In the experiments, coatings without a dark contour were obtained at 3; 4; 6 and 8 A/dm’. At the current 
density above 4 A/dm”, dendrites were formed at the corners of the cathode. Their appearance means that the delivery 
of nickel ions is the limiting stage of the process [12]. The concentration polarization is eliminated under increasing the 
concentration of nickel salts in the electrolyte. At the nickel sulfate concentration of 300 g/l and high current densities, 
dendrites were not formed at the cathode. 

The change in the current efficiency of the coating depending on the cathode current density is due to the ap- 
pearance of dendrites at the cathode. With increasing current density, the current efficiency increases, reaches the max- 
imum value at 3 A/dm’, and then decreases sharply (Fig. 1). 
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Fig. 1. Dependence of current efficiency of CEC Ni—ZrB, on cathode current density. Powder concentration in electrolyte is 40 g/l; 
electrolyte temperature is 50 °C; pH 2 


When washing the samples, the dendrites are partially lost; as a result, the current output values of the coating 
decrease at the current density above 4 A/dm’. The effect of increasing the upper value of the current density in the 
presence of ZrB, can be explained by the high degree of its dispersion. The particle size of the powder is comparable 
with the thickness of the electrical double layer, therefore, when particles move near the cathode surface, they can fall 
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into the double layer zone. In this case, colloidal nickel compounds are involved in the cathode layer. They are formed 
in the electrolyte under the preparation and electrodeposition, and they prevent sudden alkalization in the process of 
electrolysis. 

The buffer capacity of the electrolyte remains almost unchanged with the addition of ZrB,. pH in the cathode 
layer seems to decrease due to the addition of new portions of ZrB, from the electrolyte volume. The probability of this 
phenomenon is confirmed by the fact that the upper value of the operating current density does not increase in the pres- 
ence of ZrB, micropowders (5 um) in the electrolyte of similar composition. 

The optimal concentration of ZrB, in the electrolyte was determined from its content in the coating and from 
the microhardness values of the CEC. It was established that in the CEC under study, the content of ZrB, reached max- 
imum (1.1 wt.%) when its concentration in the electrolyte was 40-60 g/l. With an increase in the cathode current densi- 
ty at all the studied concentrations of ZrB, additive into the electrolyte, the content of ZrB2 in the CEC increases to the 
maximum limit (1.1 wt.%). CEC with stable microhardness characteristics is formed when ZrB, concentration in elec- 
trolyte is 30—40 g/l. 

For production and cost reasons, the powder concentration should not exceed 40 kg/m’. Otherwise: 

- mixing efficiency of the electrolyte decreases; 

- it is difficult to exclude the presence of particles in the workpiece surface areas which are hard-to-reach for 
the electrolyte circulation; 

- quality of coatings deteriorates. 

The dependence of the ZrB, content in the CEC on its concentration in the electrolyte is considered. It is shown 
that when the concentration of ZrB, particles is 10 g/l with the grain size of 2—9 um, the distance between the ZrB, par- 
ticles in the coating is approximately 1/3 of the distance between the particles in the electrolyte. As the ZrB, concentra- 
tion increases by 10—20 times, the distance between the particles in the electrolyte decreases, and the coating remains 
almost unchanged. The concentration of ZrB, in the electrolyte can be significantly reduced, and the content of particles 
in the coating will not decrease. Thus, the use of electrolytes with a high concentration of ZB, is impractical, since the 
maximum content of ZrB, in the coating is reached at its lower concentration in the electrolyte. 

The microhardness of the CEC with ZrB, rises with increasing the powder concentration in the electrolyte. It 
reaches its highest value (11 GPa) at the ZrB, UDP concentration of 50—60 g/l and the current density of 5—6 A/dm’. At 
low current densities (0.5 A/dm’ and 1 A/dm’), the microhardness gradually grows with increasing the concentration of 
the second phase. At the current densities above 2 A/dm’, the microhardness of the CEC increases sharply, and already 
at the ZrB, concentration of | g/l it exceeds the microhardness of pure nickel by 1.5—2 times. 

The increase in the CEC microhardness in the presence of the ultradisperse second phase is explained not only 
by the inclusion of solid and colloidal nickel particles in the coating, but also by a decrease in the crystalline grain at 
high current densities, as well as by the dispersion hardening of the nickel matrix. There is an opinion that a necessary 
condition for the dispersion hardening of the coating is the occurrence of ultrafine particles in it, which prevent recrys- 
tallization and the formation of coarse grains [18]. On the one hand, submicrometer-sized particles inhibit the growth of 
crystal grains; on the other hand, they form agglomerates in the electrolyte and cannot penetrate into in the coating as 
separate particles. Evidently, some part of the second phase is included in the sediment (CEC) in the form of individual 
particles and stimulates the formation of a finely-crystalline structure of a nickel matrix with high hardness. The proba- 
bility of this phenomenon cannot be denied if for no reason than because particles of the same nature, differing in the 
dispersion degree, exhibit different ability to increase the microhardness of the CEC. 

The possibility of obtaining CEC with a low content of the second phase and high hardness takes on particular 
importance for increasing the wear resistance of the surface, since, in this case, the matrix must have specific elasticity. 
A significant increase in the volumetric content of solid particles causes the fragility of materials. Thus, the use of UDP 
to obtain wear-resistant coatings is based on the practice. Equally important is the reduction of the consumption of the 
second phase, and, consequently, the process efficiency. 

It is known that the hardness of coatings allows for indirect assertions about some physical and mechanical 
properties of the CEC. In some cases, it correlates adequately with wear resistance [4]. However, maximum microhard- 
ness does not always correspond to high wear resistance. The latter depends on the interaction of various factors of the 
process, as well as on the specific friction conditions. Hence, the wear resistance of the CEC with ZrB, was compared to 
chrome galvanic coatings traditionally used as wear resistant [19]. At this, both lubricated and unlubricated friction 
modes were considered. According to the results obtained, the service lifetime of the parts with Ni — ZrB, coating is 2—5 
times longer than that of the parts with hard chrome plating. 
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Internal stresses of the Ni—ZrB, composite coatings decrease with an increase in the ZrB, concentration in the 
electrolyte and the coating thickness. Differences in the internal stress values of nickel CEC grow with increasing the 
sediment thickness (Fig. 2). 
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Fig. 2. Dependence of internal stresses on coating thickness and zirconium boride concentration in electrolyte, g/I: 
1 — 30; 2 — 5; 3 — without second phase. Cathodic current density is 2 A/dm”; electrolyte temperature is 50°C; pH 2 


For galvanic coatings, the following regularity is peculiar: the finer the crystal grain and the higher the hard- 
ness, the greater the value of the internal stresses. In the presence of ZrB, the microhardness of the coatings increases 
sharply, and the internal stresses decrease. To explain this common pattern mismatch, it is necessary to consider the 
effect of the second-phase particles (particles of ZrB, UDP and colloidal particles formed in the electrolyte) on the pre- 
cipitation hydrogenization ratio. The discharge of the hydrogen ions (and, consequently, a more complete removal of 
the molecular hydrogen) occurs faster on the second-phase particles than on the surface of nickel. So, the degree of pre- 
cipitation hydrogenization decreases, which leads to a reduction of the internal voltage of the coating. 

Discussion and Conclusions. In summary, the authors of the paper propose the following composition of 
nickel sulfate-chloride electrolyte-colloid for deposition of the nickel-based CEC and nanostructured ZrB,, g/l: nickel 
sulfate - 250; nickel chloride - 60; a-aminoacetic acid - 20; ZrB, UDP - 30-40. The electrolysis mode is as follows: 
cathode current density is 2-4 A/dm”; electrolyte pH is 2—2.6; temperature is 50° C; it is necessary to mix the electrolyte 
during the electrolysis process. 

The electrolyte composition is developed; and the process condition for the application of the nickel-based 
CEC with the addition of ZrB, is established. The physical and mechanical properties of the obtained Ni — ZrB, CEC are 
studied: microhardness, wear resistance, and internal stresses. 

CEC has high microhardness (10—11 GPa with the indentation load of 100 g), which is 1.5—2 times higher than 
the microhardness of pure nickel. The wear resistance is 2—5 times more than that of chrome coatings. As microhard- 
ness increases, the internal stresses of the Ni — ZrB, CEC decrease. 

Considering the listed parameters, the proposed method can be used for surface hardening of parts for specialty machin- 
ery and industrial equipment. 
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Dynamic damping under introduction of additional couplings and external actions 
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'?3 ipkyTCKHi TOCyapCTBeHHbIit YHHUBepcuTe;T IyTeli coobmenua, r. Upxytcx, Poccuiickaas Dexepanna 


Introduction. The dynamic interaction features in mechanical 
oscillating systems, whose structure includes additional 
couplings, are considered. In practice, such cases occur when 
using various optional mechanisms and motion translation 
devices under the formation of technical objects. The study 
objective is to develop a method for constructing mathematical 
models in the problems of dynamics of the mechanical 
oscillating systems with optional devices and features in the 
system of external disturbing factors. 

Materials and Methods. The techniques used to study 
properties of the systems and the dynamic effects are based on 
the ideas of structural mathematical modeling. It is believed 
that the mechanical oscillating system, considered as a design 
model of a technical object, can be compared to the 
dynamically equivalent automatic control system. The 
mathematical apparatus of the automatic control theory 1s 
used. 

Research Results. A method for constructing mathematical 
models is developed. The essential analytical relations for 
plotting oscillating systems are obtained, which enable to form 
a methodological basis for the integral estimation and 
comparative analysis of the initial system properties in various 
dynamic states. Dynamic properties of the two-degree-of- 
freedom systems within the framework of the computer 
simulation are investigated. The implementability of dynamic 
oscillation damping mode simultaneously in two coordinates 
with the joint action of two in-phase kinematic perturbations in 
the mechanical oscillating systems is shown. 

Discussion and Conclusions. The possibilities of new dynamic 


effects, which are associated with the change in the system 


“The research is done within the frame of independent R&D. 


Beedenue. PaccMaTpuBarotca OCOOeHHOCTH MHAaMM4eCKHX 
B3aMMOJIEHCTBHM B MCXAHM4CCKHX KOJICOaTeJIbHbIX CHCTEMAX, 
B CTPyKType KOTOPBIX VMMEeIOTCA JOMOJIHUTeIbHbIC CBA3H. 
IIpaxtuyecku Takue cuTyallMu BO3HHKaKoT ~——sdTIp 
HCHOIb30BaHHH B (OPMHMpOBaHHH TeXHHYeCKHX OOBEKTOB 


pa3sJIM4HbIx JOMOJIHHTCJIBHbBIX MCXaHW3MOB HU yCTPOHCTB AJA 


lIpeoOpa30BaHiA JBYKeHUA. Lemp HiccIeqOBaHHA 
3akIouaeTCHd B  pa3pa0oTKe MeTOZa  mOocTpoeHua 
MaTeMaTHMyecKHX MOseIei B 3ama4axX MHaMMKH 


MeXaHH4eCKHX KOJIeCOaTeJIbHbIX CHCTeEM C JOMOJHUTCIbHbIMU 
YCTpOMCTBaMH HM OCOOeCHHOCTAMH B CHCTeMe BHeIIHUX 
BO3MYIMarOllUx (akKTOPOB. 

Memooobl, UctiOsIb3yeMble JIA UCCHeOBaHHA CBOMCTB CHCTeM 
M W3y4eHHA WHHaMMyecKux IPDeKTOB, OCHOBAHbI Ha UeAx 
CTPyKTYpHOro 
IlomaraeTca, 4TO MexaHH4ecKOH KOIeOaTebHOM CHCTeMe, 


MaTeMaTHYeCKOLTO MOJICJIMpPOBaHHdA. 


paccmMaTpHBaeMow B- KkayecTBe pacueTHOM CXeMBI 
TEXHHYeCKOrO OOEKTA, MOXKHO COMOCTABUTb IKBUBAJICHTHY!O 
B J{MHaMW4eCKOM OTHOINICHHH CHCTeMy aBTOMAaTHYeCKOrO 
ynpapieHna. Ucnonmb3syetca MaTeMaTH4eCKHH allmapaT TeopHu 
aBTOMATHYeECKOLO yiIpaBJIeHHA. 

Pe3yibmamoi ucciedoeanusa. Pa3padOoTaH MeTO], TWocTpoeHua 
MaTeMaTHyecKHxX MoOjsieneHw. TlomyueHbl HeoOxoaUMBIe 
aHaJIMTHM4eCKHe COOTHOIMMCHHA AIA MOCTPOCHHA YACTOTHBIX 
yjMarpaMM KOJIeOaTeJIbHBIX 


CUCTCM, TIOSBOJIAFOWUIHE 


C(OpMUpPoOBaTh MeTOOJIOrH4eCK yO OCHOBY TA 
HMHTerpasIbHOM OICHKH MW CpaBHMTesIbHOrO aHasiv3a CBOMCTB 
HCXOJIHBIX CHCTCM B Pa3JIMUHbIX JJMHAMHYCCKHX COCTOAHHAX. 
IIpopeyeHbI UCCIeqOBaHHA JWMHAaMM4eCKHX CBOMCTB CHCTeM C 
BYMA CTeHeHAMM CBOOObI B paMKaXxX BbIYHMCMTeCIbHOrO 
MoyenMpoBaHHua. JloKka3aHbl BO3MO%KHOCTH peallM3alvu B 
MeXaHH4eCKHX  KOJIeOaTeJIbHbIX CHCTeCMaX  pe2xKHMOB 
JMHAMHYECKOLO FallleHHA KOUeOaHHM OJHOBPeMeHHO TIO JBYM 
KOOPJMHaTaM TIPH COBMeCTHOM JIeMCTBMM JIBYX CHHWa3HbIx 
KHHeEMaTH4eCKHX BO3MYIeHHH. 
O@6cyacdenue u 3akilro4eHuA. OTMeY4eHbI BO3MO%KHOCTH 
IIPOABJICHHA HOBbIX JMHAaMMyecKUuX IPdeKTOB, KOTOPBIe 
CBA3aHbI C W3MCHeCHH€M CTPyKTypbl CHCTeMbI IIpH 


OlpeyesIeHHbIX (MOpMax JIMHaMMYeCKHX B3aMMOeMCTBUM. 
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structure under certain forms of dynamic interactions, are PaOota lipesicTaBIaeT MHTepec JIA ClelMauHCTOB B OOsacTu 


noted. The study is of interest to experts in machine dynamics, AWHaMMKH MallIHH, POOOTOTeXHHKH, MeXaTPOHHKH, HaHO- 


Me30MeXaHHKH. 

robotics, mechatronics, nano and mesomechanics. 
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Introduction. Dynamic oscillation damping (DOD) is widely used in practice to change the state of technical 
objects and the local action on the forms of component interaction of mechanical oscillatory systems. In the papers [1-— 
3], the features of the approaches to the implementation of the dynamic oscillation damping modes were considered, 
options of the design and engineering solutions and calculation technique of the system parameters were proposed. 

A variety of dynamic tasks predetermines a wide variability of the proposed solutions, within the framework of 
which the dynamic features of the protected objects, the conditions of external disturbances of the original system and 
the design and engineering modes of the dynamic oscillation damping are considered [4—8]. 

DOD 1s used in the tasks of protecting instrumentation systems and tooling [9, 10]. However, some aspects of 
evaluating the dynamic properties of vibration protection systems have not obtained a proper level and detail of 
representations in the formulation of research tasks and parameter analysis for dynamic absorbers. This may be due to 
the introduction and use of additional links, considering the characteristics of external influences, as well as the effect of 
the simultaneous joint action of several external disturbing factors. 

In this paper, we develop a method of constructing mathematical models and the formation of dynamic 
oscillation damping effects in the chain mechanical oscillating two-degrees-of-freedom systems. 

I. Background. The original system 1s shown in Fig. 1. It represents two inertia members (m, and m2), which 
are interconnected by the elastic elements with rigidities (k,, k2, k3) and additional links in the form of motion 
translation devices (MTD) with reduced masses (L;, L2, L3). The system performs small oscillations under the action of 
external in-phase harmonic effects. Resistant forces are not considered. 











Fig. 1. Flowsheet of technical object in form of two-degrees-of-freedom lumped system 


To describe the motion, the coordinate system (yj), y2) 1s used in a fixed basis. Suppose that the kinetic and 
potential energies of the system are determined by the following expressions: 


Lua... &. l 8 1 ] ee 
T= 7 my +5 MVn ro —z) Fig a —y,)y ere: —z,)’, (1) 
ey ee, eee 2445 D 
9 iW, Z,) 5 OL, 5 es) : ( ) 


The set of differential equations of the system motion in the time domain is obtained on the basis of the 
formalism of the Lagrange equations of second kind, and it has the form: 

(™m, +L +L, )¥, + y,(k, +k,)—y,L, —k,y, = Lz, +K,Z, +O, ) (3) 

(m, Ls +L,)¥, + y,(k, +k,)—yL, —k,y, = 1,2, +k,Z, +Q, ° (4) 
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After the Laplace transformations under zero-initial conditions [11], the equation set (3)—(4) can be represented 
by a structural mathematical model in the form of a circuit that 1s dynamically equivalent to the automatic control 
system [12, 13] shown in Fig.2. 


| 1 





Fig. 2. Structural mathematical model of mechanical system shown in Fig. 1 


II. Mathematical Models Development. The type of transfer functions depends on the nature of external 


disturbances, that is, on whether the disturbances are force (O, and O, ) or kinematic (z, and Z, ). Further on, suppose 


that the force external actions (Q, and Q, ) have connectivity determined by the ratio 


O, =0-9,, (5) 
where a is coefficient of connectivity of external in-phase harmonic actions. 
For kinematic effects (Zz, and Z, ), it is assumed that 
Z,=B-Z, (6) 
where f is coefficient of connectivity of kinematic disturbances. 
The connectivity coefficients (a and B) can have positive, negative, and zero values. Special cases of the effect 
of external disturbances can be considered through zeroing o and £ properly. 
1. The case of joint force perturbation at B #0 (O, =0 and O, = ()) is considered. The transfer functions of the 


system in this case take the form: 
Yi (Lp +k, )[(m,+L,+L,)p° +k, +k,]+B(L,p° +4,)\(L,p +k) 


W.(p) = 7 
ae AD) 
y, L,p’ +k,)[(m,+L,+L,)p? +k, +k, ]+(Lp’ +k, (Lp +k, 
W,(p)= 2 =k sP tk )p Ca ie Ae a (8) 
QO;=Q,=0 Z, A(p) 
where 
A(p) =[(m, + £,+L,)p° +k, +k) [m+ £,+L,)p +k, +k]-(hp thy (9) 


is a system frequency standard equation. 

When considering the transfer functions (7), (8), it is assumed that the dynamic mode of oscillation damping is 
determined by the conditions for zeroing numerators (7), (8). Coordinate may cause two frequencies of dynamic 
oscillation damping. Along the coordinate (y,), the occurrence of two DOD frequencies is possible. Along the 


coordinate ( y, ), it is also possible to implement two DOD modes due to the virtual existence of roots of the biquadratic 


frequency equation. 
Assuming that the variable factor is B (coefficient of connectivity), a frequency diagram can be developed 
considering the following frequencies: 


1. partial frequencies: 


k +k 
n, =—__+_., (10) 
m +L, +L, 
k, +k 
ee (11) 
m, +L, +L, 
2. critical frequency of interpartial communication: 
k 
2 . 
Nowy = 7 (12) 


3. DOD frequencies determined from the solution of equations in y, coordinate: 
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P{L.(m, +L,4+L,)+BL,L,]+ p' [Lh (kh, +k,)+k,(m,+L,+L,)+P(A,L, +4,L,)]+ (13) 
+k (k, +k,)+Bk,k, =0: 


and in y, coordinate: 
Pp [BL,(m,+£,4+L,)+LL,]+ p [BL (4, +4,)+ BA, (m, + £,4+L,)+4,L,+kL,]+ (14) 
+Bk,(k, +k,)+k,k, =0. 

To construct a frequency diagram, the following parameters of the model problem are accepted: m, = 10 kg; mz 
= 10 kg; k,; = 5000 N/m; ky = 10000 N/m; 3 = 15000 N/m; L, = 5 kg; L, = 10 kg; L3 = 10 kg. 

The modes of dynamic oscillation damping in y, and y, coordinates are determined not only by the parameters 
of elastic inertia members, but also by the specificities of the external actions formation, in particular, by the 
connectivity value. 

From the equations (13), (14), corresponding DOD frequencies can be found showing that two dynamic 
oscillation damping frequencies can be found in each of y, and y, coordinates. The frequencies values, as it follows 


from (13) - (14), depend on the connectivity coefficient of kinematic disturbances (B). Fig. 3 shows the system 
frequency diagram. 
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Fig.3. System frequency diagram shown in Fig. 1 
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In the diagram, the solid line ( ) corresponds to (8) dependency graph. Since the frequency equations 


(13) - (14) are biquadratic, each of the equations has two roots. This is displayed in two graphs. For @,,,,,(8) graph of 


dependences, the solid line is marked with special symbols (xxx), and for the second root, respectively (°°°°); , sum (PB) 


2 
2 WMH 


dependency graphs are touching int. (£). Again, o;_ (8) dependency graphs are represented by dashed lines (— — —). 


©; ,,(B) dependency graph consists of two non-contiguous blocks. Mutual intersections of «,,,,(B) and @,,,,,, (B) 


2 1H 
dependency diagrams occur in tt. (1), (2), (3). Each of the considered points determines the amplitude-frequency 
characteristics associated with the DOD modes features. 

In the usual formulation of studying the dynamic oscillation damping, that is, under the action of a single 
perturbing factor correlated with a certain coordinate, one DOD frequency is determined in the two-degrees-of-freedom 
system. Such a frequency is determined by the partial frequency values of that system block, the movement of which 
demonstrates the dynamic oscillation damping (that is, “zeroing” the value of the corresponding coordinate). 

Under the action of several simultaneous disturbances, it becomes possible to implement two DOD modes in 
each of the coordinates for the system as a whole. When additional links are introduced into the system, in particular, on 
the basis of the MTD, specific properties occur when the dynamic oscillation damping becomes possible simultaneously 
in two coordinates. 

III. Comparative analysis of dynamic properties of the systems in the DOD modes. 

1. Fig. 4 shows frequency-response characteristics (FRC) of the system, which manifest themselves under the 


conditions corresponding to the intersection of  (B)and ow, (B) graphs in t. (1) in Fig. 3. The intersection 


1 quan 2 quH 
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) in Fig. 4 corresponds to ~*~ (@) dependency graph; the dotted 


1 





corresponds to the case when £6 = 0. The solid line ( 


line (+**" ) corresponds to ~* (@) graph. The parameters of the system as a whole are also shown in Fig. 4. 
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, (), 23 (@) 








Fig. 4. Frequency-response characteristics of the system for parameters determined by t. in Fig.3 


In t. (1) in Fig. 4, the frequency of the DOD mode is determined when y, and jy,coordinate values are 


simultaneously “zeroed out”. 


It is possible to implement another DOD mode in jy, coordinate corresponding to ~(@) graph in t. (2) (Fig. 4). At 


1 
this, the FRC reflect the properties of two-degrees-of-freedom systems. As follows from the FRC, it is possible to 
implement two DOD modes in tt. (1) and (2) in y, coordinate. For the FRC corresponding to y, coordinate, it is also 


possible to create two DOD modes at the double intersection of ~ () graph by the abscissa line after ,... 
Z 


1 


eigenfrequency. Thus, with non-degenerate FRC for each of y, and y, coordinates, two DOD modes can be 


implemented; while at one of the frequencies, there 1s simultaneous DOD in two coordinates. 


; I ae ; 
2. Fig. 5 shows the system FRC atB = 5? from which it follows that it becomes possible to restructure the 


system when one degree of freedom degrades at certain parameter ratios. In this case, in each of y,and y, coordinates 


(points (1) and (2) in Fig. 5), it is possible to implement the DOD modes, but this occurs when the system “degrades”; 


under the conditions when FRC  — © acquire limiting properties. 
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Fig. 5. Frequency-response characteristics of the system at B = = corresponding to t. (3) in Fig.3 


3. Fig. 6 shows the system FRC at B = a which corresponds to t. (2) in Fig. 3. For the given state of the system 


determined by the parameter values, a structural transformation of the system is also characteristic; the system 
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“degrades” to the status of a one-degree-of-freedom system. For the system on a whole, in each of the coordinates, the 
DOD implementation is possible, which corresponds to t. (1) and t. (2) in the graphs shown in Fig. 6. 
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Fig. 6. Frequency-response characteristics of the system at B = 34 corresponding to t. (2) in Fig. 3 
In the high frequency area (@ — ©), the system acquires limiting properties; at this, in comparison with the 
l ee . oe 1 
previous example, when B = ae the oscillation amplitudes ratio will have a different sign at B = on Consequently, 


changes in system parameters under the simultaneous action of two forces can change drastically the dynamic 
properties of the mechanical oscillatory systems. 


Conclusion. The simultaneous action of external disturbances in the presence of additional links in the system, 
implemented by the MTD under the kinematic disturbance, can have a great impact on the change in the dynamic 
properties of the mechanical oscillatory systems with several degrees of freedom. So, the authors in this paper have 
obtained the following research results: 

1. A technique of constructing mathematical models based on the use of methods of structural mathematical 
simulation, 1n which the mechanical oscillatory system is compared to the dynamically equivalent automatic control 
system, is proposed; 

2. A technique for constructing frequency diagrams that allows for the integral estimation of the 
interdependence of frequency characteristics is proposed for the case when the system parameters and their perturbation 
conditions for various power factors change; 

3. The analytical conditions for the implementation of the DOD modes simultaneously in two coordinates 
under the action of two interlinked disturbing factors are obtained; 

4. The possibilities to control the structural states when the original mechanical oscillatory system can change 
the number of degrees of freedom and the system of its dynamic properties are proposed. 
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Pa3padoTKa WM@poBsoro 7ABOMHHKA cTaHka c UITY Wa ocHOBe MeTOAOB MallIMHHOrO O0y4eHHA 


to. I. Kaéanaun’, JI. A. [lararun’, M. C. AHOCOB., A. M. Ky3bMnuinna’ 


| 2.3.4 HypxeropouckKHit rocyqapcTBeHHbIii TeEXxHHYeCKHM yHuBepcutert, r. HwKHuii Hosropog, Poccniicxaa Dezepanna 


Introduction. It is shown that the digital twin (electronic pass- 
port) of a CNC machine is developed as a cyber-physical sys- 
tem. The work objective is to create neural network models to 
determine the operation of a CNC machine, its performance 
and dynamic stability under cutting. 

Materials and Methods. The development of mathematical 
models of machining processes using a sensor system and the 
Industrial Internet of Things 1s considered. Machine learning 
methods valid for the implementation of the above tasks are 
evaluated. A neural network model of dynamic stability of the 
cutting process is proposed, which enables to optimize the 
machining process at the stage of work preparation. On the 
basis of nonlinear dynamics approaches, the attractors of the 
dynamic cutting system are reconstructed, and their fractal 
dimensions are determined. Optimal characteristics of the 
equipment are selected by input parameters and debugging of 
the planned process based on digital twins. 

Research Results. Using machine learning methods allowed us 
to create and explore neural network models of technological 
systems for cutting, and the software for their implementation. 
The possibility of applying decision trees for the problem of 
diagnosing and classifying malfunctions of CNC machines is 
shown. 

Discussion and Conclusions. In real production, the technolo- 
gy of digital twins enables to optimize processing conditions 
considering the technical and dynamic state of CNC machines. 
This provides a highly accurate assessment of the production 
capacity of the enterprise under the development of the pro- 
duction program. In addition, equipment failures can be identi- 
fied in real time on the basis of the intelligent analysis of the 


distributed sensor system data. 


"The research is done within the frame of independent R&D. 


Beedenue. B ctarTbe 10Ka3aHo, YTO IHpoBon ABOMHNHK (91eK- 
TPOHHbIM MaciopT) craHka c UITY pa3pabaTpiBaetca Kak Ku- 
Oepdu3sn4eckasd CHCTeMa. 

Ilemb padOoTbl — co3qaHve HeMpoceTeBbIx MOJIeJIeH, ompeye- 
JIAFOMIMX (PyHKIMOHUpOBaHne cTaHka c HITY, ero mpon3Boan- 
TeJIbHOCTb H JJMHAMMYECKYIO YCTOMYUBOCTb MPH pe3saHHu. 
Mamepuaabi u memodoi. PaccMaTpuBatOTCA BOTIPOCbI co31a- 
HUA MaTeMaTH4eCKHX MOJeeH MpoleccoB MexaHvyeckon 
OOpaOoTKH C UCHOIb30BaHHeEM CHCTeMbI CeHCOpoOB HU IIpo- 
MBIIIJICHHOrO HHTepHeta Belle. OleHeHbI MeTOJIbI MallIMH- 
Horo OOy4eHHA, NOAXOMANIMe JIA peasM3allMH Ha3BaHHBbIx 
3sanau. IIpemmoxena HelpoceTeBad MOJeIb WHHaMM4uecKou 
YCTOHYMBOCTH Tpollecca pesaHuA, MO3BOJAFOMIad ONTHMV3H- 
POBaTb Mpolecc MexaHuyeckKon OOpaooTKH Ha JTalle TEXHOJIO- 
THY4eCKONU NOATOTOBKH Tpou3BoycTBa. Ha ocHOBe NOAXOOB 
HeJIMHeCMHONM TMHAMMKH PeKOHCTPyHpoOBaHbl aTTPaKTOPbI IH- 
HaMH4ecKONM CHCTeMbI pe3aHHA U OlpeyeseHbI UX PpakTasib- 
Hble pasMepHocTH. BaiOpanbl ONTHMAaJIbHbIe XapakKTepHCTUKU 
oOopy0OBaHHA 10 BXOJJHbIM WapaMeTpaM HU OTIayke m1aHupy- 
eMOro TeXHOJIOrMYecKorO Mpolecca Ha OCHOBe LM@poBBIx 
JIBOMHUKOB. 

Pe3yibmamoi ucciedoeanuu. VUcriosmb30BaHve MeTOOB Ma- 
IWMHHOrO OOyYeHHA TMO3BONHIIO CO34aTb HW UCCIICOBaTb 
HelpoceTeBble MOJ[eJIM TEXHOJIOrMYeCKHX CHCTeM OOpadoTKH 
pe3aHveM HU IIporpaMMHoe obecreyeHue JIA UX peaM3allun. 
Iloka3aHa BO3MO2%KHOCTb IIPHMeCHeHHA JIepeBbeB PellleHHM NWA 
3aqa4uu jMarHOCcTHuKH UW KlaccHukallun HeucnpaBHocTer 
cTaHKkos c ITY. 

O6cyacdenue u 3akirouenua. B peasibHOM IIpOW3BOCTBe Tex- 


HOJOFHA UWAMPOBbIX JBOMHUKOB TMO3BOJIAeCT ONTHMHM3MPOBaTb 
Pe@xUMBI OOpaovoTKH Cc y4eTOM TeXHWYeCKOrO HW JMHAaMW4eCcKO- 
ro cocTtosHua cTaHKoB c UITY. Dro obecneuMBaeT BbICOKO- 
TOYHYIO OLCHKY MPOM3BOJICTBCHHBIX MOINHOCTeM lpewIpns- 
THA TIP COCTaBJICHHM IIpOH3BOJICTBeHHOM IIporpaMMBI. Kpome 
TOrO, Ha OCHOBE HHTEJJICKTYAJIbHOrO aHasIu3a JAHHBIX CHCTe- 
MBI paciipeJIesICHHBIX CCEHCOPOB MOXKHO BbIABUTb HeMCIIpaBHo- 


CTH OOOpyOBaHHA B PexKUMe peasIbHOrO BpeMeHH. 
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Introduction. The single core control platform proposed in [1] serves as the basis for the development of a 
new generation of the processing equipment control systems, and it also provides improvement of the numerical control 
systems (CNC) for digital productions. Open CNC systems of machines with long-lived computational resources and 
high speed processing of a large database (DB), embedded neural processor modules and communication modules with 
industrial Internet can be such a platform. In the latter case, we are talking about the applicability of cloud technologies 
for processing large amounts of data on enterprise servers (local networks) and providers. All this will create the basis 
of intelligent control for a wide range of the CNC process equipment. 

Since the nineties of the 20th century [2], the authors of the presented research have been developing software 
products for intelligent control of the processing equipment. In particular, the criteria were proposed for evaluating the 
dynamic stability of the cutting process, based on the methods of nonlinear dynamics and fractal analysis of the vi- 
broacoustic emission (VAE) signals. The authors have created a single platform for enhancing various technological 
equipment with CNC systems. It is implemented through embedding high-performance computing modules and deep 
learning of artificial neural networks into the CNC systems using CUDA nVidia technologies. In the framework of this 
work, the CNC process equipment accomplished with sensors is studied. It applies cloud technologies for collecting 
and processing information using the developed techniques. Such systems will be considered as cyber-physical ones [3]. 

It is assumed that the platform developed by the authors will become the basis for digitalization at all levels of 
the enterprise. To this end, it should not only analyze the data of equipment, systems, and devices, but also use the in- 
formation obtained in this way to reduce time on market launch of new products, to increase production flexibility, 
product quality and efficiency of production processes. 

Digital twin is a new word in the modeling of equipment, processes and production planning [3]. This is a set 
of mathematical models that describe reliably the processes and interrelationships on a single object and within an entire 
production enterprise using big data analysis and machine learning. 

The leader in the application of digital twins is Siemens [3]. According to its definition, a digital twin is an en- 
semble of mathematical models. They characterize various equipment conditions, production and business processes 
over time, in accordance with current production conditions. Neural network circuits hold a specific place among such 
mathematical models; i.e. a neural network model of a process or product is its digital twin [2]. 


The base unit of the digital production is a cyber-physical system (CPS) [1]. Its high adaptive and intellectual 





capabilities are due to the following features: 

- associated perception of information, 

- ongoing training 

- assessment of the current condition and forecasting of the future one. 

CPS is able to analyze multidimensional data considering even hidden factors of the real production. Based on 
this data, 1t can autonomously solve optimization problems and make right decisions. Therefore, CPS is the key element 
in creating a digital twin. In this case, the following is meant by the digital twin: 

- a set of mathematical models that characterize various equipment conditions, production processes over time, 
according to the current production conditions; 

- detailed 3D assembly object models, reflecting connections and interactions between nodes. 

From this point of view, a digital twin can be considered as a digital identity of the CPS, an electronic passport, 


which records all the data on the materials being processed, the performed technological operations, and testing. 
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Currently, digital twins are mainly created for commercial purposes. They operate successfully in the oil and 
gas industry. At the same time, there is no data on the twins of the equipment of machining industries based on CPS in 
the literature. 

Digital twins based on CPS can be obtained through the following: 

- traditional analytical approaches based on the mathematical description of physical processes; 

- modern statistical methods including machine learning. 

Materials and Methods. Machine learning methods used to build statistical models can be divided into three 
groups: regression analysis models, classification models, and outlier detection models [4—6] (Table 1). 

Table 1 


Basic Machine Learning Methods 


Regression analysis Classification Outlier detection 
Logistic regression Support vector method 
Decision Trees Forest Principal component analysis 


Decision hungle 
-_Neuralnetworks | _Bayespointmachine | 
One-vs-All 
/ Ordinal regression [| ___Newalneiworks S| SSS 


The selection of the machine learning method depends on the size, quality and nature of the data, as well as on 





the type of tasks to be solved. The existing methods require different computational performance, and they have varying 
degrees of accuracy. Mostly, they are evaluated by the possibility to achieve accurate approximation of the data and to 
identify the boundaries in the data space. The method of artificial neural networks (ANN) is the most universal and ac- 
curate one, it enables to operate a large amount of data and to build non-linear dependences. Neural networks use a 
large number of settings, which opens up possibilities for creating highly accurate models of processes operating in the 
regression analysis, classification, and outlier detection modes. 

The basic method for analyzing model quality is cross-validation. It allows for evaluating the statistical quality 
of the source data through constructing and comparing several models obtained from different training and verification 
samples. When building models of complex objects and systems, it is required to reduce the data dimensionality and 
eliminate the effect of multicollinearity of variables. The solution to such problems is possible due to the application of 
the principal component method, which represents multidimensional data in the form of a limited number of compo- 
nents. Such a generalized approach can be applied to eliminate retraining of models. The basic methods of machine 
learning are presented in Table 1. 

To improve the quality of the model, bagging and boosting algorithms are used. It is about building not one 
model, but a set of models solving the same task. The result of this work is a kind of integral assessment of the probabil- 
ity of some event. This assessment can be presented as a synergistic result of a set of models, each of which individually 
performs poorly. Thus, a digital twin can be a set of statistical models that use various combinations of the machine 
learning methods and have passed through various stages of verification and improvement. 

So, in general, a digital twin based on the CPS is a multifactor model of the equipment [1], which includes an 
ensemble of electronic (1.e., neural network) models. In this case, the following neural networks will be decisive: those 
of dynamic stability of the cutting process, of cutting forces, and of the machined surface roughness. Cutting forces 
cause elastic pressing in the “tool - workpiece” system, which specifies errors in shape and size. 

To solve this problem, a complex of statistical models was developed using machine learning methods. The re- 
sulting models are the basis for the digital twin of the CNC lathes. They provide solving the regression analysis prob- 
lems to predict the dynamics of the cutting process under various machining conditions, roughness of the machined 
surface, and cutting forces. In addition, these models help to solve the classification problems for evaluating the ma- 


chine current condition. A process diagram of creating a digital twin 1s shown in Fig. 1. 
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Development of equipment digital twin using machine 
learning methods and cloud technologies 


CNC lathes Industrial Internet of Things +» 
/ ‘Ne 
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Fig.1. Process diagram of creating a digital twin of CNC machine 

Consider in more detail the neural network model of the cutting process dynamics. The training set was ob- 
tained on the basis of telemetry data collected under machining through the distributed sensor system of the VAE sig- 
nal, dynamometer, and Industrial Internet of Things technology (IloT). The number of sensors, their type and spatial 
orientation were determined in accordance with the equipment layout. The presence of a large number of sensors with 
different spatial orientation in the system is explained by the heterogeneity of materials and structures, the signal propa- 
gation features, as well as the possibility of restructuring the oscillatory system under operation. Hence, the application 
of a heterogeneous sensor system could obtain the most complete dynamic picture of the processes in the n-dimensional 
“state — time” space. The standard TCP/IP protocol and JSON (JavaScript Object Notation) were used as the data trans- 
ceiving protocol within the IloT network. 

Analog and digital band pass filters, as well as wavelet filters, were used for the signal preprocessing. The ap- 
plication of wavelet filters provided: 

- elimination of the noise term impact in the VAE signals; 

- identification of periodic and chaotic components based on the entropy factors. 

Data pre-processing has been performed on an “edge” device (Edge), which calculates the values of the VAE 
signal parameters and forms data packets for pushing them to the virtual cloud storage. The development of mathemati- 
cal models based on the machine learning methods requires structuring and data separation. However, considering the 
IloT technology aspects, the construction of relational DB is not always possible; therefore, an approach based on 
NoSQL technologies was applied. For this, data storage and processing were implemented on a virtual server that simu- 
lates the operation of a computing cluster. A special distributed, scalable file system was deployed on this server. 

Freeware utilities, libraries and the Hadoop project frameworks (one of the most successful and common big 
data technologies) were used as a basis. In particular, the Hadoop Common software shell manages the HDFS distribut- 
ed file system and the HBase database [4]. To perform the distributed computing and processing of large volumes of 
data, MapReduce was used, which provides automatic parallelization and distribution of tasks on a cluster. nVidia 
CUDA graphics processors were used to speed up the training of neural network models. 

Research Results. A recurrent INS with the sigmoidal neuron activation function was selected to create a sta- 
tistical model of the cutting dynamics. The data obtained under the industrial equipment operation was used as a train- 
ing set. During training, the backpropagation algorithm was used. At the end of the learning process, the model obtained 
was verified on the basis of the mean square derivation values. Fig. 2 shows the monotonous decrease of error. Moreo- 


ver, its distribution is normal, and it is near zero value, which implies good quality of the model obtained. 
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Fig. 2. Assessment of model quality obtained: diagram of training error distribution (a); dependence of root-mean-square 


error on learning iteration (5) 


The ANN obtained (Fig. 3) consists of 17 input neurons that acquire information on processing conditions and 


on the previous dynamic condition. 
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Fig. 3. Neural network model of dynamic stability of cutting process 


The output layer of neurons gives information on the dynamic state of the cutting system. The hidden layer is 
formed by 38 neurons. This ANN is able to simulate the cutting system dynamics for various machining conditions at 
given instants. 

To assess the adequacy of the neural network model, a simulation of the machining process was carried out at a 
CNC turning center. A shaft of stainless 12X18HIOT steel fixed on both sides in the center, as shown in Fig. 4, a, was 
used as the workpiece. 
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Fig. 4. Neural network simulation results: 
treatment scheme (a), prototype cutting process (5), real-world cutting process (c), cutting process after optimization (d) 


The modes were defined in accordance with the production standards for machining parts of this type on a 
CNC lathe. The simulation of the tool motion relative to the workpiece was provided by a variable on the ANN input 
layer responsible for the machine support position along Z axis. This variable spacing was 0.01 mm at the distance from 
0 to 350 mm. Thus, the dynamic state of the cutting system was obtained for 35 thousand tool positions on Z axis. The 
recorded dynamic states are reflected in the corresponding spectrogram (see Fig. 4, b): there are two dominant self- 
oscillation frequencies (3.5 kHz and 6.5 kHz) and two phase transitions. 

To confirm the simulation results, a full-scale experiment was carried out with similar processing conditions 
and recording of the VAE signal. According to the experimental data obtained, a spectrogram of the machining process 
was constructed (see Fig. 4, c). It also shows two dominant frequencies in the 3.5 kHz and 6.5 kHz areas and two phase 
transitions. As a result, the simulation error in different machining areas was from 3% to 7%. As can be seen from Fig. 
4, the proposed production standards do not allow obtaining a product of a given quality. This may be the cause of de- 
fects, and, considering the workpiece runout, — of tool breakage. 

The use of the above-mentioned norms leads to an inadequate assessment of the plant capacity and, according- 
ly, to a poor-quality formation of the production program. Removal of this disadvantage at the design-engineering stage 
is fraught with additional time and financial expenditures. Therefore, the modes factored in the norms of production are 
adjusted. For this, the optimal values of the technological conditions vector are calculated using the ANN and the gradi- 
ent descent method. Objective optimization functions were the VAE signal amplitudes, fractal dimension and infor- 
mation entropy, which should tend to a minimum. The processing speed and the bevel value on the face surface were 
selected as optimization parameters. Upon completion of the optimization process, new treatment modes were obtained: 
V = 90 m/min; S = 0.1 mm/rev; fy = 0.2 mm. 

These modes were also evaluated according to the results of a full-scale experiment with recording and pro- 
cessing the VAE signal. As can be seen from Fig. 4, d, optimization of the processing conditions made it possible to 
increase the dynamic stability of the cutting process through reducing the amplitude of self-oscillations by half, and the 
fractal dimension (Df= 1.1). Hence, the dynamic machining quality was improved without loss of productivity. 

In a similar way, models were obtained for predicting the roughness values of the machined surface and cutting 
forces. Cutting forces determine the processing accuracy, causing an elastic pressing of the part and tool, and, conse- 
quently, the part shape error. These neural network models are multilayered artificial neural feedforward networks (Fig. 
5). 
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Fig. 5. Distribution of amplitude values of sensor signal: As_1 (a), As 2 (6), As 3 (c) 

Neural network models as part of the digital twin are applicable not only under real production conditions, but 
also when choosing and justifying the purchase of equipment. At this stage, it 1s possible to select the optimal character- 
istics of the equipment and debug the planned process on digital twins. 

To assess the current condition of the equipment, a classification model was developed using decision trees [5]. 
Values of the VAE signal amplitude at different points of time from three sensors located on the tool (As_1), at the 
front spindle support (As_3) and at the rear center (As _ 2) were selected as attributes for decision-making. Three condi- 
tions were used as classes: Normal, Tool wear, and Abnormal. Any condition, when the quality of the treated surface 
and the operation of the equipment were unsatisfactory, was implied as abnormal. The criterion for tool wear was con- 
sidered bevel on the rear surface (h = 0.15 mm). Upon receiving the learning sample, an exploratory data analysis was 
performed. 

Fig. 5 shows the correspondence of the distribution of the VAE signal amplitude values of the As_1 sensor to 
the three groups of conditions. The data obtained from this sensor differentiates well the normal condition and tool 
wear. However, the values corresponding to the abnormal condition have intersections with data from the other groups. 
Therefore, the information from the As_1 sensor cannot be used with high accuracy to classify conditions. 

Fig. 5, b, characterizes the distribution of signal amplitude values from the As 2 sensor. The data obtained 


have low information value due to two shortcomings: 
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- intersections of the values of the abnormal condition group with other classes; 

- weak differentiation of the normal condition and tool wear (this is due to their close location). 

Fig. 5, c, shows the distribution of signal amplitude values from the As 3 sensor. The information from this 
sensor identifies well the values of the abnormal condition, while there are no intersections of the values of all classes of 
conditions. Thus, the exploratory data analysis enables to draw a conclusion on the possibility of obtaining a statistical 
model for the classification of the CNC machine condition. 

The CART algorithm was used to train the model [3]. It divides the data into specific subsets at different lev- 
els. The objective minimization function of informational entropy is used as a separation criterion. Hence, each new 
level has a more specific content of attributes. A set of all levels obtained is represented as a hierarchical tree structure 
consisting of nodes and branches, and the branching conforms to the standard logical rules. 

From the entire set of attributes, a value of the division condition for the root node is chosen, which better min- 
imizes the informational entropy values. Then, branches are built from each node received, and a similar division into 
new nodes takes place. Such branching may continue until all values of the training sample are classified. In general, 
such a tree will represent a complex branched structure that ideally classifies the current sample. However, it may show 
very low accuracy when working with new data. Therefore, it is necessary to optimize the tree size. It can be carried out 
using a cross-validation, on the basis of which several candidate trees are built, and, in total, the tree that showed the 
best result is selected. Alongside with this, pruning is performed until this procedure increases dramatically the error. 

Fig. 6 shows the model obtained for classifying the equipment condition based on the readings of the VAE 


sensors using a decision tree. 





Fig. 6. Decision tree for classifying equipment condition 


The first condition chosen is the VAE signal amplitude readings from the As_ 3 sensor. If the VAE amplitude is 
greater than or equal to 2.3, then the condition is uniquely abnormal. This judgment is supported by the exploratory data 
analysis. In the contrary case, the subsequent splitting occurs, and the data from the As_1 sensor determines the classifi- 
cation. If the VAE amplitude of the As_1 sensor is less than 2.1, the condition is considered normal; if it is more, then 
we are talking about the tool wear. 

As can be seen from Fig. 6, the decision-making does not take into account readings of the As 2 sensor. Low 


information content of the signal from this sensor was also detected under the exploratory analysis. 
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High accuracy of the determination of the abnormal operation mode based on the analysis of the signal from 
the As_ 3 sensor may be related to the place of its installation (front spindle support). Most likely, the dynamic processes 
in the spindle assembly caused this condition, which also affected the machining quality. On considering the nonlinear 
processes of oscillations dispersion in the elastic machine system, this information reached other sensors in already 
heavily distorted form. 

At this, the As_1 sensor mounted on the tool specifies its wear with high accuracy. This can be explained by 
the dominant role of the cutter in the general oscillatory system. 

The As_ 2 sensor mounted on the rear center records information for classifying the condition with the least ef- 
ficiency. This can be explained by the distance of As 2 from the main dynamic processes. 

Hence, the described method can be applied not only to the task of classifying the equipment condition, but to 
the determination of the optimal number of sensors and their locations. However, it is necessary to consider the possi- 
bility of restructuring dynamic processes in the cutting system when changing processing conditions and natural tear 
wear of machine parts. 

Discussion and Conclusions. The use of digital equipment twins under manufacturing planning reveals bot- 
tlenecks in the technological operations, improves product quality, and reduces tool breakage and equipment failure 
risks. Digital twins can optimize the processing modes taking into account the technical and dynamic condition of each 
production unit. This approach provides a highly accurate assessment of the production facilities under the development 
of the production program. Besides that, real-time equipment faults are detected through the intelligent data analysis of 


the distributed sensor system. 
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PrucK-OpveHTHpoOBaHHbili NOAXO/] B CHCTCMEe «Me pCcOHaJI-MexaHhsMbI-IIpOnsBOICTBeCHHad Cpeta» Ha o0beKTAaX, 


oK ok ok 
IKCHIYAaATUpYIOWNx OalleHHbie KpaHbl 


E. B. Eresbekan., A. A. Koporknii’, 9. A. Tanusosa’, A. A. Kuwxu6a0B 


1,2,3,4 


Introduction. The paper discusses the applicability of a risk- 
based approach in the “personnel-machinery-production 
environment” system at the facilities running tower cranes 
through the introduction of IT-technologies that provide open 
communication of the staff, management team of the operating 
organizations, and the National Supervisory Authority, to 
reduce the accident rate. 

Materials and Methods. An example of a _ hazardous 
production facility running tower cranes is given. Materials on 
the analysis of operational status of tower cranes within the 
framework of the current legislation in the field of industrial 
safety are used. The necessity for innovations and 
transformations, one of whose methods is the risk-based 
approach allowing for the implementation of all required 
levels of control, is identified. 

Research Results. The stages of introduction of the risk- 
oriented approach for tower crane operators implemented 
through IT-technologies using Web-applications on safety 
management under the tower crane operation based on the 
hazard analysis and risk assessment in gamut, algorithmically 
associated with the electronic block key of its local security 
system, are determined. 

Discussion and Conclusions. The application of a risk-based 
approach in the ““personnel-machinery-production 
environment” system at the facilities running tower cranes 
through the introduction of IT-technologies will ensure the 
proper operation of all departments and each participant of the 
operation, as well as provide Supervisory bodies with an 
opportunity to access information on the operation for the 


implementation of control and supervision functions remotely. 


” The research is done within the frame of the independent R&D. 


JloHcKon rocyjapcTBeHHbIM TexXHuYecKHH yHUuBepcuteT, r. Pocros-Ha-Jlony, Poccuiickaa Denepaina 


Beedenue. B- ctTaTbe paccMaTpuBaeTCH BO3MO2%KHOCTb 
IIPHMeHCHHA PHCK-OPHeHTHPOBaHHOTO MOAXOAa B CHcCTeMe 
«TepCOHalI-MexaHH3MBI-IIpOu3BOJCTBeHHaA cpea» Ha 
OOBeKTAX, HCHOJIb3YIOUINX OallleHHble KpaHbl MOcpeCTBOM 
BHeapeHua IT-TexHonornv, oOecme4uBaIOlIMx OTKPbITOe 
PYKOBOAMIero 
opraHu3aluuw WM TocysapcTBeHHbIx 


B3aMMOJIeHCTBHe  mepcoHasia, cocTaBa 
9KCIUIYaTMp yrOWMx 
Haji30PHBIX OPraHoOB B I[esIA4X CHWMKCHHA YPOBHA aBaPHMHOCTH. 
Mamepuanot u memooovi. IlpupeyeH mpumMep omacHoro 
IIPOH3BOJICTBCHHOrO OObEKTA, SKCIIYaTUpylollero OallieHHbIe 
KpaHbl. Mcnob30BpaHbl MaTepiasibl aHasiv3a COCTOAHHA 
SKCIIyaTauMu OallleHHbIX KPpaHOB B YCJIOBHAX JTEHCTBYFOIIero 
3aKOHOAaTeJIBCTBa B OOACTH MpOMBILIICHHOM Oe30MacHocTH, 
BbIABJICHAa HeCOOXOJMMOCTb YHHOBAIMA UW Wpeodpa30BaHni, 
OHM H3 COCTABJIAFOUIMX Me€TOJOB KOTOPBIX ABJIACTCA PUCK- 
OPHeCHTHPOBaHHbI NOAXON, WO3BOJIAIOWIMM peasIM30BaTb BCe 
HeOOXOJMMbIe YPOBHH KOHTPOWA. 

Pezyibmamol ucciedoeanua. OnpefesIeHbl ITAabl BHeAPeHHuA 
noyxona Aa 
OalllieHHOrO KpaHa, peasIM30BaHHOrO ocpeycTBom == IT- 


PUCK-OpHeHTUpOBaHHOro MallIMHUCTOB 
TeXHOJOrMH C UCHONb30BaHHeM Web-IIpnioxeHuu 10 
ylipaBleHHtO Oe30MacHOCTbIO TIpH KCITyaTauMu OallleHHoro 
KpaHa Ha OCHOBe aHasIMv3a OMacHOCTeli HM OI|eHKH pucKa B 
CBA3AHHOW—C 


I[BeEeTOBOH TaMMe, aJITOPHTMUdeCKH 


IIICKTPOHHBIM KJIFOUWOM OIOKHPOBKU ero JIOKaJIbHOM CHCTeMBI 


Oe30acHOCTH. 

O6cyocoenue u  3akiiouenud. IlpumMeHeHue  pucK- 
OPpHeHTHpOBaHHOTO MOsxXoa B  CHCTeMe «ilepcoHal- 
MC€X@HH3MBI-lIpOM3BOCTBCHHad cpejla> Ha  OOBeKTAX, 
9KCIUIyaTHpyrollx $ OallieHHbIe KpaHbl $MOCcpeyCTBOM 


BHezpenusa IT-TexHoONOrMi, MO3BONMT OOeCHe4YHTb HajexKHOe 
(YHKUMOHMpoBaHve BCexX Hospa3szenieHui UW KaxOrO 
yuacTHukKa paOoT MO 9SKCIiIyaTalMu, a TakoKe lpesOcTaBHT 
BO3MOXKHOCTb Ha30pHbIM OpraHaM T1OJIyuaTb JOCTyN K 
WMHDOpMalwH TO 9KCIIyaTalMu WIA OCYIIeCTBIeHHA 
KOHTPOJIBHO-HaJ30PHBIX PYHKUMM JMCTAaHIMOHHO. 
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Introduction. The first decades of the XXI century were marked by serious achievements in science and 
technology, modernization and improvement of digital technologies, production processes and engineering capacities of 
the production potential of Russia, which serves to strengthen the economy of the state and to improve the living 
standards of the citizens. 

Against the background of the above achievements and prospects, issues of the operational safety and the 
professional competence of personnel at the production facilities, as well as supervision of the safe operation of 
hazardous production facilities (HPF), are challenge problems. 

The HPF condition as a whole is estimated by such indicators as accident rate and injury. The level of safety of 
the HPF is under the state control and is supported by Federal Law No. 116 [1]. Since its adoption in 1997, there has 
been a transition from the deterministic to probabilistic risk assessment. New techniques for analyzing technically 
challenging and hazardous objects are introduced [2]; it is supposed to evaluate them by the parameters of risk, 
survivability and safety (Fig. 1). 





Traditional 
Methods 


Fig. 1. History of development of methods for calculating HPF 


Research Results. Consider HPF running lifting devices, in particular, tower cranes. At present, such an object 
is ranged in the fourth hazard class (low hazard). In this regard, it is exempted from periodic inspections by the federal 
executive authority in the field of industrial safety, and the responsibility for complying with the industrial safety 
requirements rests in full with the owner or operator. The discussion of adequacy of the internal control over 
compliance with industrial safety requirements without the participation of supervisory authorities was held at different 
levels. 

The recent rise in accidents and injuries at the facilities running tower cranes, has led to the tightening of 
administrative measures, namely: 

- Order of Rostekhnadzor No. 146 of April 12, 2016, changes the Federal norms and rules for lifting devices 
under the operation of tower cranes. Article No. 141 assigned for the participation of an inspector of 
Rostechnadzor in the commission for decision-making on the possibility of bringing tower cranes in 
operation|[3 |; 

- Instruction of the Government of the Russian Federation No. AX-I19-682 of February 8, 2017, “On the 
organization and conduct of unscheduled inspections of the enterprises running tower cranes from 2017 to 
2019” regulates the procedure of random checks. 
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Additionally, the draft law proposes to classify enterprises using tower cranes to the third hazard class. This 
will allow for scheduled oversight activities against these facilities every three years. The operating organizations 
themselves also consider this measure appropriate. 

The need for transformation and innovation in the field of industrial safety has come to a head. Significant 
measures are as follows: the introduction of a risk-based approach during the organization and implementation of all 
types of state control (supervision) from 2018 (introduced by Federal Law No. 246-FZ of July 13, 2015) [4]; using the 
information and communication technologies to perform the duties in an electronic form during the implementation of 
state control (supervision). The regulating document is the Decree of the Government of the Russian Federation No. 
482 of April 21, 2018, “On the State Information System of “Typical Cloud Solution to Automation of Control 
(supervisory) Activities” [5]. 

The following tasks of the state policy in the field of industrial safety are formulated, among other documents, 
by the Decree of President of the Russian Federation “On the Fundamentals of the State Policy of the Russian 
Federation in the Field of Industrial Safety for the Period to 2025 and for the Future” of May 06, 2018: introduction of a 
risk-oriented approach under organizing the federal state control (supervision) in the field of industrial safety; 
annulment of obsolete, redundant and duplicate industrial safety requirements [6]. 

The operation of tower cranes is carried out, as a rule, at construction sites geographically distant from the 
workplaces of managers and experts, which allows the crane operator to make autonomous decision on the work 
commencement and performance. Periodic on-site visits of specialists responsible for the operable state do not allow for 
an objective situation assessment, which leads to numerous emergencies and injuries. 

The functions of a tower crane driver, indicated in the production instruction, presuppose a certain algorithm of 
his actions before work starting, related to the determination of the technical readiness of the crane for safe operation. 

The statistics of most accidents with tower cranes and the results of their investigations reflect the failed state 
of cranes (nonoperable condition of separate units), which is not recorded in the duty log during the inspection of the 
crane operator. Consequently, a specialist responsible for the operational condition of the lifting device has no 
information of it because of non-compliance with the instructions by a crane operator. The enterprise management 
learns about the incidents that have already happened as fait accompli — during the investigation of the causes of 
accidents (injuries). 

Multi-level control and interaction of all participants engaged in the industrial safety under the operation of 
tower cranes (Fig. 2) can be considered the only solution to safety problems. 
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Fig. 2. Personnel and managers interaction algorithm 
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The application of a risk-based approach provides the methods of organization and implementation of all levels 
of control, including measures of preventive control of the compliance with safety requirements concerning facilities 
Operating tower cranes. This approach can be implemented through IT technology using Web applications. These 
applications help to manage safety when operating a technical facility based on hazard analysis and risk assessment in 
color gamut. An algorithmic connection with the electronic block key of the local safety system is made. Thus, the 
compliance with the standards established by law in the field of HPF industrial safety is ensured through the remote 
technologies. 

One of the steps of introducing a risk-based approach will be the implementation of electronic checklists for 
meeting industrial safety requirements through individual mobile devices. At this, information on the fulfillment of the 
security requirements is sent in real-time via the Internet to individual mobile devices to the manager, specialists and 
service personnel (Fig. 3). 





Fig. 3. Electronic checklists in the form of mobile applications 


The next step involves monitoring and confirming the fulfillment of security requirements, determined by the 
duty and working instructions through individual mobile devices, where prompts in the form of tables on the controlled 
action and its scheduled time should popup (Fig. 4.) 


‘f 
+ 


Ti 
jftei 


fiyeaey 
\ 


"i rgi4s) ! 
ait Thi 
falegtate. 
Kk * ~ ~ 


fi i 





Fig. 4. Mobile screen with popup checklists 


The next step includes processing of the language data from managers, specialists and service personnel, as 
well as digital information obtained from the coordinate protection of the crane in the cloud space using a risk 
assessment algorithm based on the theory of fuzzy sets of L. A. Zade. 

The software to assess risk under the crane operation is displayed as a color gamut, for example, of three colors 
according to the scheme: 

e =red - “work prohibited”; 
e yellow - “work permissive, actions required”; 
e = green - “work permissive ”. 
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Fig. 5. Risk assessment under crane operation in color gamut 


Crucially, the access to information on existing violations should be implemented using an appropriate login 
and password, which is confidential information. This will enable to establish causes of the unfavorable cases (accident, 
incident), as well as a concrete employee who has failed to fulfill his duty or working instructions (Fig. 6). 





Fig. 6. Steps of displaying risk-based approach algorithm on mobile devices 


One of the stages of the implementation of the risk-based approach is to monitor the crane operation by the 
video surveillance signal through a remote CCTV camera — a “black-box” recorder. Its signal is transmitted via the 
Internet to the cloud space where the video information is stored (Fig. 7). 
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Fig. 7. Displaying panorama of object on mobile device using video surveillance 


Conclusion. The application of the risk-based approach in the “personnel-machinery-production environment” 
system at the facilities running tower cranes through the introduction of IT-technologies will ensure reliable operation 
of all departments and each participant of the operation. This will provide supervisors with the ability to remotely 
access to operating information for the implementation of the control and supervisory functions. Reducing the number 
of accidents and injuries during the tower cranes operation can be achieved only under the condition of open interaction 
of the personnel, management and supervisory authorities. 
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Deep convolution neural network model in problem of crack segmentation on asphalt images 
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Mojeab rtyOokon CBepToOuHOn HeHpoOHHO ceTH B 3aaue CerMeHTallMu TpelwH Ha H300paxKeHHAX acalbra 


Poa 


b. B. Cooonp’, A. H. CosoBbes’, Il. B. Bacuaes® , JI. A. ToaKkon3nna*” 


1, 2, 3,4 


Introduction. Early defect illumination (cracks, chips, etc.) in 
the high traffic load sections enables to reduce the risk under 
emergency conditions. Various photographic and video moni- 
toring techniques are used in the pavement managing system. 
Manual evaluation and analysis of the data obtained may take 
unacceptably long time. Thus, it is necessary to improve the 
conditional assessment schemes of the monitor objects through 
the autovision. 

Materials and Methods. The authors have proposed a model of 
a deep convolution neural network for identifying defects on 
the road pavement images. The model is implemented as an 
optimized version of the most popular, at this time, fully con- 
volution neural networks (FCNN). The teaching selection 
design and a two-stage network learning process considering 
the specifics of the problem being solved are shown. Keras 
and TensorFlow frameworks were used for the software im- 
plementation of the proposed architecture. 

Research Results. The application of the proposed architecture 
is effective even under the conditions of a limited amount of 
the source data. Fine precision is observed. The model can be 
used in various segmentation tasks. According to the metrics, 
FCNN shows the following defect identification results: IoU - 
0.3488, Dice - 0.7381. 

Discussion and Conclusions. The results can be used in the 
monitoring, modeling and forecasting process of the road 


pavement wear. 


Keywords: artificial neural networks, defect identification, 
segmentation, road pavement, cracks, IoU, Dice. 


JloHcKon rocyjapcTBeHHbIM TexXHuYecKHH yHuBepcuTeT, r. Pocros-Ha-Jlony, Poccuiickaa Denepana 


Beedenue. CBoeppeMeHHoe ycTpaHeHve jedeKTOB (TpelHH, 
CKOJIOB HW Ip.) Ha y4acTKax MOBBILNeCHHOM Harpy3KU OPO2xKHO- 
TO MOJOTHA MO3BOJIAeCT CHH3HTb PUCK BO3HHKHOBCHHA aBa- 
PHMHBIX cuTyauMi. B HacToallee BpeMA [JIA KOHTPOJIA COCTO- 
AHUA JOPOXHOTO MOKPbITHA MPHMCHAIOTCA pa3sIM4HbIe MeTO- 
bl POTO- U BUAeOHAOIFOAeHHA. OleHKa UW aHasM3 ToOuy4eH- 
HbIX JJAHHbIX B PYYHOM PexXKUMe MOFyT 3aHATb HeJOMYCTHMO 
MHOrO BpeMeHH. Takum o0pa30M, HeOOxOAHMO COBepIlicH- 
CTBOBATb IIpOlesypbl OCMOTpa HU OLCHKH COCTOAHUA OOBECKTOB 
KOHTPOJIA C TOMOIIbIO TEXHHYECKOLO 3PeHHA. 

Mamepuazoi u mMemoool. ABTOpaMH T1peAIO%*KeHa MOJeJIb Ply- 
OoKOH CBepTOYHOM HelMpOHHOM ceTH IA UAeHTUduKalNNn 
jedeKTOB Ha “300paxKeHHAX JOPOxXHOrO NoKpbITHA. Moyesb 
peamm30BaHa KaK OMTHMV3HpOBaHHbIN BapwaHT HanOosiee 
NOMYAPHbIX Ha JJaHHbIM MOMCHT TIOJIHOCTbIO CBePTOUHbIX 
HeMpOHHbIx ceTeli (FCNN). [loka3aHo nmocTpoeHue oOy4aro- 
Ie BEIOOPKH U TBYXSTAaHbIM Mpolecc OOyyeHHA CeTH C y4e- 
TOM ciielu@uku pemlaemon 3agaun. Ja WporpaMMHoln pea- 
TM3alWH pesOKXeCHHOM apXHTeKTypbI UCTIOJIb30BaIINCb 
(bpetimBopku Keras u TensorFlow. 

Pe3yibmamoi ucciedoeanua. IIpumMeHeHve lpeox*xeHHO 
aPXUTeKTypbI IPeKTHBHO axKe B YCJIOBHAX OFrpaHHyYeHHOro 
OOBeMa HCXOJHBIX JaHHbIx. OTMeueHa BbICOKAA CTeMeHb T0- 
BTOPACMOCTH Pe3yJIbTATOB. 

MojeIb MOxKeT ObITb HCHOUb30BaHa B pa3IMYHbIX 3ayjaqax 
cerMeHtTalMu. CorsacHo MeTpukam, FCNN noka3pIBaeT cJie- 
Myrolue pesyIbTaTbl UAeHTUduKaluU edextos: IoU — 
0,3488, Dice — 0,7381. 

O6cyorcoenue u 3aknroy4enua. Ilosy4eHHble pe3yJIbTaTbl MOLyT 


ObITbh HCIOJIb30BaHBbI B Ipouwecce MOHUMTOPHHTa, MOJICJINPOBa- 
HWA UW WpOrHOSHpOBaHHA TpOWeCCOB H3HOCaA JOPO?KHBIX IIO0- 


KpbITHH. 


Ku1104eBbIe CJIOBa: HCKYCCTBCHHbIC HeMPOHHble CeTH, UCH- 
THuKauUMA WedeKTOB, CerMeHTAallHA, TOPOxXHOe MOKpbITHE, 
TpemuHbl, loU, Dice. 





“The research is done on RFFI grants nos. 19-08-00074, 18-31-00024, 16-01-00390. 
“E-mail: b.sobol@mail.ru, solovievarc@Wgmail.com, lyftzigen@mail.ru, podkolzinalu@gmail.com 
“™ PaOota BbINOJHeHa Ip NossepxKe rpaHTos POOM 19-08-00074, 18-31-00024, 16-01-00390. 


Information technology, computer science, and management 


ON 
Oo 


http://vestnik.donstu.ru 


Vestnik of Don State Technical University. 2019. Vol. 19, no. 1, pp. 63—73. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk JTonckozo 2zocyOapcméennozo mexnuyeckozo ynuegepcumema. 2019. T. 19, Ne 1. C. 63—73. ISSN 1992-5980 eISSN 1992-6006 





For citation: B.V. Sobol, et al. Deep convolution neural net- Oopazeu Ona Yumupoeanua: Mojenb riryOoKou cBepTOuHOH 
work model in problem of crack segmentation on asphalt im- HeMpoOHHOHM ceTH B 3afaue cCerMeHTalluH TpellluH Ha 
ages. Vestnik of DSTU, 2019, vol. 19, no. 1, pp. 63-73 w300paxeHuax achanbta/ b. B. Codonp [u gp.] // BecrHuk 
https://do1.org/10.23947/1992-5980-2019-19-1-63-73 JloHcKoro roc. TexH. yH-Ta. — 2019. — T. 19, Ne 1. — C. 63- 


73. https://doi.org/10.23947/1992-5980-2019-19-1-63-73 


Introduction. The road pavement wear requires regular monitoring. Efficient monitoring strategies allow for 
the timely detection of problem areas. This approach improves significantly the efficiency of road maintenance, reduces 
maintenance costs and provides continuous operation. Technologies for identifying critical features of the condition of 
the road surface have developed from the manual photofixation methods to the high-speed digital technology [1]. 

Russia is one of the five countries with the longest automobile roads. To provide photo- and video monitoring 
of such a large-scale infrastructure special systems are required. They should be reliable and easy-to-use on-line sys- 
tems. It seems clear that in this case it is not about solutions that involve the data analysis in manual mode. Such an 
approach is unacceptable due to the considerable time spent on information processing and low quality of the analysis. 

The authors of the paper propose a new technological solution in the field of machine learning. Its implementa- 
tion enables to automate the process of road surface quality assessment. To this end, the convolution neural network is 
trained on the data marked up manually. The system learns to recognize and evaluate the basic failure modes of the 
monitor objects by this means. 

Research Digest. Many papers describe the improvement of algorithms for detecting defects in the structural 
components and infrastructure facilities. To solve this problem, computer vision capabilities are widely used. Their con- 
tinuous improvement is supported by the development of sensing technologies, hardware and software. However, it 
should be recognized that at present, computer vision is used sparingly. This is due to many factors, including: 

- heterogeneity of defects, 

- variety of types of surfaces, 

- complexity of the background, 

- junctions. 

The authors of a number of publications investigate automated methods for detecting cracks in images and 
propose their own solutions [2—9]. Some works consider the specifics of monitoring objects of road transportation infra- 
structure [10, 11], as well as bridges and structural systems [12, 13]. 

Until recently, mostly manual monitoring techniques were used to solve these problems, such as: 

- morphological operations [13], 

- analysis of geometric features [6], 

- application of Gabor filters [14], 

- wavelet transforms [15], 

- building of histograms of oriented gradients (HOG) [16], 

- texture analysis, 

- machine learning [4]. 

However, currently, the listed tool is used more and more seldom. It is replaced by a global spread of neural net- 
work technologies and machine learning supported by the computational capabilities of graphic processors. 

A convolution neural network (CNN) is a multi-layered artificial neural network architecture designed specifi- 
cally for imaging [17]. In this case, a subsampling layer (pooling layer) allows local capable fields to be realized 
through convolution layers and invariance with respect to small geometric strains. 

This architecture demonstrates outstanding results in solving the following recognition problems: 

- handwritten numbers [18], 

- house numbers based on the Google Street View house number (SVHN) dataset [19], 

- road signs [20]. 

The computation power rise of graphics processors allows for the use of deeper architectures of machine learn- 
ing models [21]. Now it is possible to avoid retraining [22]. This is facilitated by the development of such modern tech- 


niques as data augmentation, regularization, etc. 


Sobol B. V., et al. Deep convolution neural network model in problem of crack segmentation on asphalt images 


Co6oonb Bb. B. u op. Modes 21yookou ceepmounou NHeuponnou cemu 6 3a0aue CezemMeHnmauuU MpewWUH Ha uz0dpancenuAXx acdanoma 





Improving convolution neural networks opens up the possibility of more efficient study and generalization of 
the features of images (for example, image classification [23], object search [24], vehicle detection [25]). 

The flexibility and prospects of deep learning for the tasks of automatic detection of pavement cracks is shown 
in [26, 27]. 

In [28], the application of neural networks for the automatic detection and classification of cracks in asphalt is 
considered. The authors propose to use the universe mean and dispersion of the grayscale values. With consideration to 
these indicators, the image is divided into fragments, and then each cell is classified as a crack. Motivation for using the 
fall weight deflectometers (FWD) to assess asphalt cracks was shown. In 98% of cases, the system detects effectively a 
crack in the image. 

In [29], the use of a neural network for the detection of defects was investigated. The advantages of the method 
of clustering pixels as objects were clarified. It enables to increase identification accuracy and reduce noise. 

In [30], the authors used a deep learning architecture that includes the VGG-16 model. It was previously 
trained to identify features that allow distinguishing between classes of images. The model demonstrated excellent 
recognition quality even under working with images from the areas unknown to it.CNN VGG-16 was used as a deep 
feature generator of pavement images. The authors taught only the last layer of the classifier. They experimented with 
different models of machine learning, and showed their strengths and weaknesses. 

In [31], the use of CNN in an applied robotics problem is shown. It is referred to the autonomous detection and 
assessment of cracks and damages in the sewer pipe. CNN filters the data and localizes the cracks, which provides ob- 
taining a characteristic of their geometrics. 

The objective of [32] is to automate the sequential detection of chips and to give numeric representation of 
damage in metro networks. For that, an integrated model that implements a hybrid algorithm and an interactive 3D rep- 
resentation is created. Chipping depth prediction is supported by regression analysis. 

The paper [33] provides an overview and assessment of promising approaches that automatically detect cracks 
and corrosion in the civil infrastructure systems. 

In [34], an effective CNN-based architecture for detecting pavement cracks on a three-dimensional asphalt sur- 
face 1s described. The CrackNet architecture provides high precision of the data processing through an ingenious meth- 
od of representing the road surface geometry. CrackNet consists of five layers and includes more than a million of 
trainable parameters. The experiments using 200 test 3D-images have shown that the CrackNet accuracy can reach 
90.13%. 

Offered Method. To identify defects in pavement images, it is necessary to determine what a defect is and 
what is not. In other words, the image segmentation should be carried out, and the appropriate classes should be identi- 
fied. Recently, this type of problem has been effectively solved through purposely developed architectures of the convo- 
lution neural networks, such as SegNet [35] and U-Net [36]. 

The specifics of pavement images involve a small range of gray shades and a minor difference between the 
background and the target object. Also, the task is complicated by noise, defects and the occurrence of foreign objects 
in the image. 

Various data sets are used to learn neural networks [7, 37]. These sets include original pavement images and 
their corresponding mask images with or without defects. Images with defects on the road surface are specific, so the 
authors offer their own simplified model of a deep convolution neural network. For image segmentation, a fully convo- 
lution neural network (FCNN) [38] with an encoder-decoder structure is proposed. The image of the road surface is fed 
to the system input, and the output is a binary image. A segmented image showing the presence or absence of defects is 
resulted. 

Architecture of Deep Convolution Neural Network. Fig. 1 shows the architecture of the proposed deep con- 


volution neural network. 
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Fig. 1. Architecture of proposed network 


The neural network consists of two parts, convolution and deconvolution. The convolution part converts the 
input image into a multi-dimensional representation of features. In other words, it performs the function of extracting 
features. The deconvolution network serves as a generator, which creates a segmented image based on the characteris- 
tics obtained from the convolution network. The last convolution network layer with a sigmoid activation function gen- 
erates a segmented image — a map of probabilities of the defect of the same size as that in the input picture. 

The first part of the network consists of five convolution layers with filter sets (256, 128, 64, 64, 64). The 
“batch normalization” (BN) tool [39] is used. ReLU, rectified linear unit, is used as the activation functions. Then, of 
subsampling (pooling) layers with a 2x2 window follow. In passing through this layer, the image is decreased by half. 
The second part of the network is a mirror image of the first. The image size needs to be restored to the original one and 
to form a probability map based on the input image features. To this end, upsampling layers are used in combination 
with the convolution layers. The proposed neural network has 10 convolution layers and 929665 trainable parameters. 

Dataset Generation. The CrackForest data set [7] is used for training the constructed model. Its augmentation 
(artificial increase of the dataset) is carried out, since the training and operation of the neural network is based on the 
path-based approach, which involves the use of randomly clipping elements of the original images. 

So, the data set consists of 117 images. It is divided into training, test and validation samples. Fragments of 
64x64 are randomly selected for each image of the training and test samples. The studies have shown that gamma cor- 
rection of images improves the neural network quality within the framework of the task. Each image fragment under- 
goes rotation, reflection, and deformation. The optimal ratio of fragments with and without defect was established as 
95% to 5%. At this, defects that occupy at least 5% of the image area are considered. The sample size affects the learn- 
ing process and the network quality. The optimal ratio was established as follows: 15 200 fragments of the training 
sample and 3 968 — of the test one. Fig. 2 shows the images and the corresponding masks used for training the neural 


network. 





Fig. 2. Images and binary masks obtained from data augmentation 
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Neural Network Training. Intersection metrics between two detections (intersection over union, IoU, Jaccard 
index) and an equivalent binary similarity measure (SOrensen-Dice coefficient) are used to train and evaluate the neural 
network. The function 1 — J is used as cost: 


(A,B) = 


|ANB| 
|AUB| 


2|ANB| 
|A|+|B 





, S(A, B) = 


The initialization of weights in the layers of the neural network is carried out through the Glorot method [40]. 
To reduce the internal covariance shift, the batch is normalized through normalizing the input distributions of each lay- 
er. Adam-algorithm (stochastic optimization method) is used for learning [41]. 

At the first stage, the neural network is trained on a small amount of data (30% of the core set) during 5 
epochs. At the second stage, the network is trained on the full amount of data during the required number of epochs. 


The learning rate varies with each epoch according to the established dependence (Fig. 3). 
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Fig. 3. Variation of network learning rate depending on epoch 


As part of the task, it was established that the optimal number of training epochs is 25 (5 epochs at the first 
stage of training, and 20 — at the second one). With more epochs, the neural network accuracy did not change crucially. 

Keras and TensorFlow frameworks were used for the implementation of the developed architecture of deep 
CNN. 

Research Results. After training the neural network, it is validated on test data. Each image fragment is fed to 
the network input, and the output is a generated map of the defect probabilities. Fig. 4 shows the results of the trained 


network performance and their comparison with the true values from the test sample. 
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Fig. 4. Results of trained neural network performance 


In Fig. 4, the data is produced in four columns: 
1) image under study, 

2) neural network output, 

3) manually identified defect, 

4) difference between 2) and 3). 


Assume that the networks are compared with the true values. The values of IoU and Dice metrics are condi- 


tioned by the specific ratios of the following factors: 
- area of the defect and area of the whole image, 
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- binary (one-bit) mask and real (4-byte) generated image. 


It should be noted that when using IoU for defectless fragments, the metric values are 0 (Fig. 5). 


Network predict accuracy histogram 
Average accuracy (Dice): 73.81%, (loU): 34.88% 
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Fig. 5. Number of images identified with specified degree of accuracy 


The prepared data set quality affects drastically the learning and the neural network output. In some cases, the 
neural network indicates a defect, although the true image does not have it, or vice versa. This affects the overall as- 
sessment of the model quality. In general, the accuracy evaluation of the neural network according to the proposed met- 
rics can be subjective, so you should not take the data from Fig. 4 as absolute. 

Some models of FCN networks are evaluated in the framework of the presented paper. The results are shown 
in Table 1. 

Table 1 
Accuracy of some neural networks 


10 layers (256, 128, 64, 64, 64, ...), 929 665 parameters Dice: 73.81 %, IoU: 34.88 % 


16 layers (32, 32, 16, 16, 16, 8, 8, 8, ...), 43 441 parameters Dice: 70.40 %, IoU: 33.24 % 





12 layers (32, 32, 16, 16, 8, 8, ...), 37 537 parameters Dice: 67.57 %, IoU: 32.12 % 


Here, the number of filters on the first part of the network is indicated in brackets. The number of filters on the 
second part of the network is flipped (see Fig. 1). 


A sliding window method with a specified step that regulates the processing speed and detailing is used to pro- 
cess high-resolution images. This is how the resulting defect probability map for the whole image is developed. Several 
images from the validation set and the result of their processing by the neural network are shown in Fig. 6. 
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Fig. 6. Validation images processed by trained FCN 


Discussion and Conclusions. A model of deep convolution neural network is proposed for identifying defects 
in pavement images. The model is implemented as a simplified and optimized version of the most popular on-the-day 
FCN networks. The techniques for constructing a training set and a two-stage network learning process considering the 
specifics of the current problem are described. The work done shows that the use of such architectures is successful 
with a small amount of the input data. Fine precision is registered. The described model can be used in various segmen- 
tation tasks. According to the metrics, FCNN shows the following results: IoU - 0.3488, Dice - 0.7381. 
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Introduction. Digital data and analytics transform the role of 
the production equipment maintenance. Analytical information 
of sensors placed on the product allows continuous monitoring 
of the production machines operation and their timely servic- 
ing. Thus, defects in technical equipment are identified, the 
analysis of which enables to develop algorithms for monitor- 
ing and forecasting, and to prevent equipment from overshoot- 
ing the limits of the safe operation. 

Materials and Methods. Basic digitalization principles and the 
digital images structure are presented. A mathematical method 
is used to describe the digital image vector and the control 
system algorithm. 

Research Results. The achievements of the known systems of 
maintenance and digitalization of various machines are sum- 
marized. The application of a dynamic digital image made it 
possible to determine the desired levels of the production facil- 
ities maintenance. An optional version of monitoring the 
equipment state within the framework of the production digi- 
talization concept is shown. It is based on the proposed algo- 
rithm for an autonomous control of the process state. 
Discussion and Conclusions. The construction of machine 
digital images in accordance with the main stages of its life 
cycle is described. The task of automated maintenance of ma- 


chine tools based on digitalization is considered. 


Keywords: monitoring, digitalization, autonomous control, 
maintenance, digital image. 
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Beedenue. Uudposbie WaHHble UW aHasMTHKa WpeoOpa3yroT 
pOJIb TEXHH4eCKOrO OOCIY*XUBAHUA MPOW3BOACTBeEHHOTO 060- 
pyaHopanua. AHamuTHyeckad HHDOpMalluA WaTUMKOB, pa3Me- 
ICHHBIX Ha H3eIMH, MO3BOUIAeT HelmpepbIBHO HaOOTaTb 
(YHKIMMOHUpOBaHHe TeXHOJIOTM4YeCKHX MallIMH HW CBOeBpe- 
MCHHO OOCJIY2KHBaTb UX. Tak, BbIABJIAIOTCA JJe:eKTbI TeXHU- 
YeCKOFO OCHALCHHA, AHAIM3 KOTOPbIX NO3BOIIAeT paspadaTbI- 
BaTb aJITOPHTMbI MOHUTOPHHTA, MPpOTHO3HpOBaHHA U lpey- 
TIpexTaTb BBIXO OOOpyAOBaHUA 3a TipeeubI HayexKHOK 
padorTsI. 

Mamepuaazoi u memoool. IIpeactaBiieHbl OCHOBHBIe MPHHIMIIbI 
UMdpoBu3allMu MU CTpyKTypa MOcTpOeHHA WMpPOBBIX OOPa3zoB. 
Ucnosb3yeTca MaTeMaTH4eCKHH MeTOJ, ONMCaHuaA BeKTOpa 
IMpoBbix OOpa30B UH alITOPHTMU3allMu CHCTeMBI yiupaBJICHHA. 
Pe3yibmamol ucciedoeaHua. OOoOmleHbl JOCTHKeHHA 3- 
BeCCTHbIX CHCT€M TEXHHYeCKOTO OOCIIWKUBAHUA U WMdpoBu3a- 
IMM pa3sIM4HbIxX MalliMH. Ucnonb30Bnanve JMHAaMHYeCKOTO 
undpoporo oOpa3a MO3BOMIO OMpeeIMTb HeOOXOMMBIe 
YPOBHU JIA MOAAepxKaHUuA padoTOcMOCcOOHOCTU TexHOJOrM4e- 
ckux OObeKTOB. I]oKa3aH BO3MO2%KHbIM BapHaHT MOHUTOPUHTa 
COCTOAHHA OOOPyTOBaHHA B PaMKax KOHIenMMU WudpoBu3a- 
IMM TpOu3BOACTBa. OH OCHOBAH Ha IIpesIO%KeCHHOM asIrOpuT- 
Me aBTOHOMHOTO yiipaBJICHHA TEXHOJIOFHYCCKMM COCTOAHMeM. 
O@cyacoenue u 3akilio4venue. ONMCaHO MOCTpoeHHe 
IMpoBbix OOpa30B cTaHKa B COOTBeETCTBHH C OCHOBHbIMU 
CTaJIMAMU eTO %KM3HEHHOYO IHkKsIa. PaccmoTpeHa 

3aqaua aBTOMaTH3HpOBaHHOTO MOAepxaHuA padoTociocob- 
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Introduction. At the present stage of development of engineering and technology, issues of the production 
digitalization, in particular, the life-cycle management of production machines, the independent management of their 
Operation and maintenance, are particularly sensitive. According to K. Schwab, Head of the World Economic Forum, 
digital data and analytics transform the role of servicing [1]. We are talking about the analytical information of sensors 
placed on the product. It helps to conduct continuous monitoring, and to study equipment operation and its mainte- 
nance. Thus, defects in the equipment are identified, the analysis of which enables to develop monitoring algorithms, to 
make forecasts, and to prevent equipment from going beyond the limits of reliable operation. This approach generally 
improves the efficiency of production processes. Description of the product operation based on digitalization allows for 
the creation of digital twins of machines for various purposes. The concept of “digital twin” includes a temporal infor- 
mation copy of an object, artificial intelligence, information technologies, and software. All these components are in- 
volved in the creation and support of intelligent digital models (DM) of complex technical products. 

Publications devoted to this topic represent research in the field of intellectual monitoring and control of the 
technical condition of machines and machine tool systems [2-14]. 

Metal-cutting machines are provided with new properties that allow them to meet the growing demands for 
speed and accuracy of cutting, reliability and safety of operation under the conditions of fast-speed processes. As in 
mechatronics in general, in machine tools, a number of organizational and technology principles related to digitalization 
is implemented. Some of them are listed below. 

Openness and controllability. A machine is an open, self-controlled system connected with the outdoor envi- 
ronment. Information on the state of the equipment, the outdoor environment is used for control, and the object behavior 
is simulated. 

Originality of properties of the entire machine system. The characteristic of a system object extends further 
than a sum of features of its components. The system obtains properties that none of its components possesses. In other 
words, the functions of the components can synthesize new properties of the system. 

Self-sufficiency of the system. The system and its elements can function independently of external systems. 

Intelligent management. Procedures for controlling an object (machine) are based on an analysis of its states 
and manufacturing situations using information technologies and knowledge processing mechanisms. 

Information and intellectual support of the machine operation includes procedures for maintenance and repair. 
The development of digital automated diagnostic and control systems involves the creation of specific software for in- 
formation processing integrated with an external network communication environment. 

Main Part. In the presented paper, the digital embodiment of an object will be called its digital image (DI). 

The machine DI corresponds to various methods of constructing digital models (Digital model). The model se- 
lection is determined by the object life cycle stage — from design to disposal (Fig. 1). 


initial D, 


perfectioning D,, 
during succeeding operational period D,. 


Fig. 1. Stages of representation of machine digital images 





According to the notation of Fig. 1 and considering the sequence of transformations, an ordered set of digital 
models of the machine state during the life cycle is as follows: 
Dn, a (Di, Ds, Dn, D,+). 


Information technology, computer science, and management 


Nl 
er 


http://vestnik.donstu.ru 


Vestnik of Don State Technical University. 2019. Vol. 19, no. 1, pp. 74-80. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk JToncxkozo zocyOapcméennozo mexnuyeckozo ynueepcumema. 2019. T. 19, Ne 1. C. 74-80. ISSN 1992-5980 eISSN 1992-6006 





According to the standard rules and test methods, the work on creating D; includes the formation of an initial 
database for assessing the machine quality. These include the following: performance accuracy factors, dynamic charac- 
teristics, thermal deformations, reliability rating estimates. The database (DB) of the initial state of the machine forms 
the basis of the temporal machine DI built on the principles of the e-Mind Machine (e-MM) [10, 11]. 

If we are talking about the machine in use, its state 1s considered under processing typical test parts or on test 
modes subject to the technical requirements and a test program. The structure of the machine DI in the operation pro- 
cesses (D,), considering the e-MM, is shown in Fig. 2. 


MODULE D, ——————— 
e-Mind Machine 







“=a, 5. Machine state 
management 
4. Monitoring 
2. Observation and 
: Control KS 
1. Machine state 
monitoring l 


Fig. 2. Structure of machine dynamic DI at service stage (D,) 








BLOCKS 





The DI state specificity 1s determined by the alternation of the operation, maintenance and repair processes. 
The DI intelligence depends on the knowledge system (KS) located in the blocks of the machine temporal DI. 

Solutions based on the embedded expert systems (ES) and information obtained from the diagnostic and moni- 
toring units are synthesized in the modules for controlling the machine technical state. They identify the measured pa- 
rameters of the machine, the process and its result, which enables to diagnose and evaluate the state of processes and 
devices. In this case, a situational assessment of the system state, which is necessary for the management, adequate to 
the current situation, is performed, and the results are predicted. Based on the decisions made, control actions are gener- 
ated corresponding to the block assignments. Considering the huge amounts of stored data, optimizing analysis tech- 
niques can be used under the system operation — for example, online diagnostics with the aid of Mel-frequency 
cepstral coefficients [12]. 

At this, it is important to apply not quantitative, but qualitative assessments and concepts based on fuzzy pro- 
cedures of processing, accumulating and applying knowledge and meta-knowledge. It is advisable to use temporal ES. 
First, they work in real time. Secondly, training and adaptation algorithms are built into such ES, and they are able to 
improve the performance of the machine subsystems. 

Under the intelligent control with the KS units during the machine operation period, complete information on 
the DI state can be represented as an ordered set (vector) of possible states of the machine system and control: 

S=(S,,5,,S,,S,,5,). 
Here, S\,S,,S,,S_,S, are, respectively, sets of states of the machining process, the tool, the technical state of the 
machine components, the state of the product (workpiece - part) and control. The set of admissible vectors (S') is 
called the admissible set of DI states, and we denote it as S'. KS synthesizes sets of time-ordered regulated actions: 


U,¢ U . They provide transformation of the vector of the DI initial condition S" eS froma hypothetically identi- 


fied admissible set of states ( § ) to the specified goal state S,¢ S . This search is based on the analysis of knowledge 
about the functionality of a specific machine tool system. Actually, KS should realize 2, indication: 
dM, SxS UU, =A,(S',S, ], 


which, out of the set of admissible controls (U ), finds such control (U,¢ U ) that corresponds to the of the initial state 
vector (Ss < §) and to the specified goal state (S_¢ S ). In this case, Ui¢ U control enables the transition of the ma- 


chine from S’ €S'state to S,¢ S state. The management apparatus (U,¢ U ) is based on the formalisms of fuzzy alge- 
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bra and fuzzy sets. This approach allows for the formulation of a plausible hypothesis about the organization of the ap- 
propriate behavior of the machine tool system. In this case, a coherent transition from one state to another is observed, 
and at each transition, information about changes in the machine DI is generated [15]: 


S,,—(S,,1,]- 
Here, S,,is DI state at the beginning of j-th transition; S$, is DI state obtained as a result of j-th transition; /, are 
changes to be made in the descriptions of DI state obtained as a result of j-th transition, U,, is control ensuring the tran- 
sition from S, , state to S’, state. 
The current state (5, ) obtained as a result of the j-th transition is determined by the state of: 
- product (workpiece - part) (S, ), 
- tool subsystems (S, ), 
- object (S;, ). 
That is 
S, =(S,,5,55,)- 


a ij 

Models of technical condition management, built on this KS base, are the framework for building a digital sys- 
tem for autonomous control of machine tools. 

Developers and researchers pay serious attention to the automated maintenance of machine operability, espe- 
cially, under the conditions of computerization of production and improvement of information support [1 ]. 

The digitalization systems developed for the technological machines maintenance require the application of in- 
tegrated and/or remote software and hardware support components. These systems are created to fulfill the following 
tasks: 

- to obtain information on the equipment state in real time, 

- to forecast the development of the state of various devices/units (D/U), 

- to give the personnel warning of emergency and other dangerous equipment conditions, 

- to carry out self-maintenance and troubleshooting, 

- to conduct maintenance or updating of the control programs in the process of equipment operation. 


Within the framework of the e-MM concept, an automated autonomous control of the technical state of ma- 
chine tools (AUTS) was created [16]. AUTS means independent automated control. Special tools and informational 
links allow it to assess the state and role of the machine. Based on these estimates, the AUTS provides certain signals or 
acts on a specific machine device to maintain or restore its operability. 

The AUTS features are listed below. 

First, the system should be equipped with a sufficient number of sensors that transmit reliable information on 
the status of units and devices. 

Secondly, when recognizing the obtained data, it is necessary to exclude noise and determine the parameters 
informatively describing the state of the equipment. 

Thirdly, the forecasting systems should flexibly adjust the indicators in real time, focusing on the information 
obtained from the diagnostic unit. 

Fourthly, the decision should include the possibility of reducing the risk from the machinery breakdown or 
wear reduction of the equipment components, and it should comply with the selected treatment program with minimal 
departure from the optimal process. 

In view of the above, the structure of the system and the autonomous control of the machines technical state 1s 
developed [15]. The algorithm of the state monitoring system is presented in Fig. 3. 
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Fig. 3. Flowchart of state monitoring system (AUTS) 


The algorithm provides the following actions. 
l. 1. Search for information on the D/U inheriting states (1.e., at the preceding stage — n - 1) 1n the 
knowledge system, and basing the n stage states into the monitoring unit DB. 


2. Obtaining current operational data of diagnostics of the D/U state at n stage. 

3. Creation (change) of the current image of the D/U state in the working memory of the knowledge system of 
n—KS, stage of the machine condition monitoring unit. 

4. Determination of the ES belonging to the KS of the monitoring unit, and tendencies of changes in the D/U 
State. 

5. Comparison of the obtained state estimates with the parameters of the fuzzy boundaries of the D/U operabil- 
ity state extracted from the database. 

6. Predictable assessment of the preservation of permissible state parameters within the limits of the D/I opera- 
bility under the condition of work on the input control program of the CNC device. 
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7. Estimation of the accessibility of fail-safety under the condition of correction of the processing program 
components. 

8. Decision-making on the accessibility of fail-safety. 

9. Decision-making on technical condition management. 

10. Correction of the D/I state without processing interrupt. 

11. Part process shutdown for the correction (restoration) of the state. 

12. Selection of the offline means to prevent (eliminate) malfunction of the D/U or external services. 

13. Conducting off-line servicing. 

14. Support call to the external services for troubleshooting of the technical condition. 

15. Diagnosis and assessment of the D/I status after maintenance. 

16. Logging and input of the list of the operations performed and the resulting assessments of the D/I state to 

the KS. Accumulation of experience in assessing the technical condition and autonomous state control in the 

KS. 

In this case, the decision-making and autonomous control of actions to maintain working capacity is made on 
the basis of the presented approach to the digitalization of machine tools, and this is the key feature of the unit for moni- 
toring the object state and of the algorithm for monitoring the AUTS state. 

Conclusion. The construction of machine digital images in accordance with the main stages of its life cycle is 
proposed. The automated maintenance of machine tools based on digitalization through the creation of autonomous 
technical condition control systems is considered. These systems transmit information on the equipment state in real 
time, signalize purposefully about the equipment state, and maintain (restore) its performance. 
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Introduction. The estimation of efficiency of methods and 
algorithms for solving optimization problems with a vector 
criterion and a set of nonlinear constraints 1s considered. The 
approach that allows proceeding to an optimization problem 
with a single objective function (i.e., an unconditional optimi- 
zation problem) after equivalent transformations is described. 
However, the objective function obtained in this way has 
properties (nonlinearity, multimodality, ravine, high dimen- 
sion) that do not allow classical methods to be used to solve it. 
The presented work objective is to develop hybrid methods, 
based on combinations of the algorithms inspired by wildlife 
with other approaches (gravitational and gradient) for the solu- 
tion to this problem. 

Materials and Methods. New methods to solve the specified 
problem are developed. A computer experiment was conduct- 
ed on a number of test functions; its analysis was performed, 
showing the efficiency of various combinations on various 
functions. 

Research Results. The efficiency of hybrid algorithms that 
combine the following approaches is evaluated: genetic and 
immune; methods of swarm intelligence and genetic and im- 
mune; immune and swarm with gravity and gradient. 
Discussion and Conclusions. The hybrid algorithms in optimi- 
zation problems are studied. In particular, decisions can be 
made on their basis under the management of compound ob- 
jects in the military and industrial sectors, in the creation of 
innovative projects related to the digital economy. It is estab- 
lished that the type of the objective function affects the result 


much more than the combination of algorithms. 


Beedenue. CTaTba NOcBALeHa Ol|eHKe 9(~PeCKTHBHOCTH MeTO- 
JOB UW aJIFOPHTMOB pellleHHA ONTHMM3al[MOHHBIX 3ajla4d C BeK- 
TOPHbIM KpHTepHeM MW CHCTeEMOM HeJIMHeCMHBIX OrpaHvyeHun. 
OnvcaH MOAXO, MO3BOJIAIOMIMM MOce NpOBeeHHuA IKBUBA- 
JICHTHBIX TWpeoOpa3s0BaHH MepeiTH K ONTHMH3al[MOHHOU 
3aqjaue C OJHOM WeueBou PyHkKiHen (T. e. K 3aage OesycJIOB- 
HOM ONTHMM3aIlMn). OHakoO TOsyaeHHad TaKHM CIIOCOOOM 
lesepadt (dyHKiMa oOOsaqaeT CBOMCTBaMH (HeJIMHeMHOCTS, 
MYJIBTHMOJ@JIbHOCTb, OBPaxKHOCTb, OOJIbINAat Pa3MepHOCTS), 
He NO3BOJIAFONIMMU MCHOJIb3OBaTb JIA ee peleHHA KIaccu4e- 
CKHe MeTOJIBI. 

[lens lpeyzcTaBileHHoro UccieqOBaHHA — pa3pa0oTaTb JA 
pellleHua TJaHHOU 3aya4du THOpUAHbIe MeTObI, OCHOBaHHbIe Ha 
KOMOMHAalMAX aJITOPHTMOB, MHCIMPHPOBaHHbIX %KUBOM Tpu- 
powou, c ApyruMu nNoAXosaMu (rpaBuTallMOHHbIM HU rpajin- 
CHTHBIM). 

Mamepuazoi u memoooi. Co3aHbl HOBbIC MeTOJIbI JIA pellie- 
HHA yKa3aHHon 3ayjauu. IIpopeyqeH KOMMBbFOTepHBIN 9KCIIepu- 
MeHT Ha psye TeCTOBBIX (PYHKIMM, BbIMOJHeEH ero aHasin3, 
NOKa3bIBalOMMM IPeKTHBHOCTh pa3sIM4HbIX KOMOMHALM Ha 
pa3sIM4HBIX PYHKIUNAX. 

Pe3yibmamel ucciedoeanua. OweHeHa 3()eKTHBHOCTb TH- 
OpHAHbIX aJIFOPHTMOB, KOTOPble KOMOMHUPyIOT ciiewyroulue 
MOAXObI: TeHeTHYeCCKHM C MMMYHHBIM; MeTOJIbI POeBOrO HH- 
TeJIIeKTa C reHeTHYeCKHMH WM WMMYHHbIMM; MMMYHHbIe U 
poeBble C TpaBuTallMOHHbIM HU rpayMeHTHBIM. 

O6ocyorcoenue u 3akoueHue. U3y4eHbl BO3MOXKHOCTH THOpu- 
HbIX aJITOPHTMOB B OMTHMH3allMOHHBIX 3ayja4ax. B 4acTHocTn, 
Ha UX OCHOBe MOFyT IIPHHUMATbCA PellleHHA pu yupaBsIeHHHn 
CJIOXKHBIMM OObCKTAMH B BOCHHOM UH IIPOMBILIJICHHOH cdepax, 
IIpH CO327aHHN MHHOBALMOHHBIX MpOeKTOB, CBA3AHHBIX C IHd)- 
POBOH 9KOHOMHKOM. YcTaHOBJIeHO, YTO BU WeseBon MyHK- 
I[MH BIIMAeCT Ha pe3yIbTAaT ropa3q0 OoOee CyIeCTBeHHO, YeM 
KOMOHHAIMA aITOPUTMOB. 
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Introduction. A general optimization problem is multicriteria with a number of restrictions in the form of 
equations and inequalities. To solve a multicriteria task, the following types of algorithms are used: 

- a posteriori (sequential assignment algorithms); 

- aprior; 

- sensing; 

- adaptive; 

- approximation algorithms based on estimated accuracy of Pareto frontier. 

The listed approaches enable to narrow down the initial set of feasible solutions to one or several points [1, 2]. 

Let us identify the features of the methods mentioned above. The use of a posteriori, aprior, and adaptive algo- 
rithms reduces the initial problem with a vector criterion to the problem with a single objective function, which can be 
solved through bioinspired algorithms or their combinations. The sensing method requires a lot of computational re- 
sources. The approximation approach is characterized by high efficiency, which is due to the possibility of parallelizing 
the computational process [2]. 

The solution to the unconditional optimization problem is based on one of two algorithms: penalty functions and 
(or) sliding tolerance. Both of these approaches are well represented in [2] and [3]. The task is to find the global opti- 
mum of the unconstrained function. 

Note features of the objective function which helps to describe a large class of technical and economic problems. 
The objective function is often non-linear, in most cases it is not differentiable, not unimodal, and it has a complex to- 
pology of the tolerance region. Therefore, to find the global optimum of the presented problem, the authors of this paper 
have created hybrid algorithms. At this, pairs can combine: 

- two bioinspired algorithms; 
- a bioinspired and a classical algorithms; 
- a bioinspired algorithm and one based on physical laws. 

Thus, a hybrid based on Fourier series and a firefly algorithm was proposed in [4]. In [5], gradient and immune al- 
gorithms are combined, and in [6], swarm and gravity algorithms are paired. The combination of genetic and swarm 
approaches [7, 8, 9] was also considered for solving various optimization problems with an objective function (includ- 
ing the problems of learning neural networks of various topology). 

The efficiency criterion of the bioinspired algorithm should be considered the number of iterations (steps) at 
which: 

- the algorithm finds a result close enough to the optimum, 

- the number of steps (time) is acceptable, 

- the required accuracy of the algorithm is provided. 

Materials and Methods. The combinations of bioinspired algorithms developed by the authors [4—10] were 
tested on Rosenbrock, Rastrigin, Griewank, and Schwefel [2] functions. It should be noted that the combinations rein- 
forced the advantages of each algorithm of the pair and leveled their shortcomings. For each function, the hybrid 
showed the best performance and accuracy compared to single algorithms. However, the hybrids that produced the most 
effective result on one test function work much more poorly with other test functions. This observation leads to the fol- 
lowing conclusion: the efficiency of bioinspired algorithms and their combinations with other algorithms depends much 
more on the objective function than on the features of the combination of algorithms. The same conclusion can be 
drawn from the NFL theorem (short for “no free lunch’) [11]. 

Consider that it is necessary to find the global optimum of some function that has the properties described 
above. There is a bank of bioinspired algorithms and (or) their combinations. And there is a bank with an objective 
function that does not have the unimodality property, with a large number of variables and parameters. What algorithm 
can solve the optimization problem for this function as smoothly as possible? For what types of functions will this algo- 
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rithm show the best result, and for which of them will it be inappropriate? It is obvious that it is impossible to develop 
such an algorithm. Algorithms may perform well on one function and be completely unacceptable for other functions. 

Concerning the structure of bioinspired algorithms, it should be noted that a large number of empirical parame- 
ters used by them does not enable to evaluate the efficiency of such algorithms and their combinations in advance. (The 
opposite example is deterministic algorithms that successfully solve problems with linear, quadratic, strictly convex, 
unimodal functions). 

For alternative estimates of bioinspired algorithms, public libraries with test problems are used, which allow 
comparing the known and new algorithms and their combinations [11, 12]. Sobolev Institute of Mathematics is one of 
the largest collections of test problems [12]. It presents various approaches to solving complex NP-complete optimiza- 
tion problems considering estimates of computational complexity. Among other libraries, narrow-focus ones can be 
noted: an open-source library for GAlib genetic algorithms [13]; libraries for building EAlib [14], Perl [15] evolutionary 
algorithms, and Java frameworks [16]. 

Research Results. The study on the features of genetic and population algorithms, performed by the authors, 
makes it possible to define the following advantages and disadvantages. 

Genetic algorithms (GA) are based on the identity of the behavior principles of biological and technical sys- 
tems. The GA use the principle of choosing the best decisions from the population available, which allows us to find the 
optimal solution to the problem. These algorithms give good diversity, as far as the information on solution point sets is 
processed in parallel. In these points, the optimum is based on the application of the objective function, and not on its 
various increments. 

Ant colony optimizations simulate the principles of vital activity of the ant colony. There, the principle of au- 
tonomous functioning of each agent is combined with the activities of the colony as a whole, which allows solving ra- 
ther complex optimization problems. The combination of ant colony algorithms with local search algorithms enables to 
quickly find the starting points of the optimum search. Such algorithms give the best results for large-scale problems. In 
terms of efficiency, they are very close to the problem-oriented and metaheuristic algorithms. This approach has shown 
good results in solving various practical tasks, for example, problems of traveling salesman, of optimal design of elec- 
tronics, etc. Due to the possibility of using adjustable parameters, ant colony algorithms are applied in solving distribu- 
tion and transport problems. 

The following positive features of ant colony algorithms should be noted: 

- for some tasks, they provide a more efficient solution than genetic ones or algorithms based on neural net- 
works; 

- the genetic algorithm stores information only about the previous generation, whereas ant algorithms store in- 
formation about the entire colony, which is more efficient; 

- random routing in the ant colony algorithm enables to exclude non-optimal initial solutions; 

- the selection of parameters responsible for changing the optimization step enables to successfully use this al- 
gorithm in dynamic applications. 

The disadvantages of ant colony algorithms include the following: 

- theoretical analysis is difficult due to a sequence of random decisions, which is caused by changes in proba- 
bility distributions during iterations; 

- the algorithm convergence time cannot be predetermined, and to solve this problem, the ant algorithm is sup- 
plemented by local search methods; 

- free parameters for work adjustment under solving a specific task are determined only experimentally. 

When identifying extremes for composite multidimensional nonmonotonic functions, the swarm algorithms 
with equal probability determine the optimal element (an element with given properties) at any iteration. These algo- 
rithms are effectively used to optimize nonmonotonic functions in NP-complete problems, including distribution and 
transport ones. They are searching the only optimal element that determines the extremum (or sets of such elements). 
The swarm algorithm implements a search in the neighborhood of the best and selected sites in parallel. At this, it does 
not have any of the drawbacks of evolutionary methods, for example, it does not require a significant amount of 
memory to store a population of solutions. 

Discussion and Conclusion. Analysis of the known hybrid algorithms has shown the following. One combina- 
tion can give a good result on some test function (for example, Rosenbrock), but this hybrid will be significantly less 
efficient than other combinations on the other function. Hence, it is fair to say that the type of the objective function has 
much more significant impact on the result than a combination of algorithms. Therefore, the study of the form and type 
of the objective function enables to select the best combination. 
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Pocculickaa Deyepanna 


Introduction. The fast program algorithm of arithmetic coding 
proposed in the paper is for the compression of digital images. It is 
shown how the complexity of the arithmetic coder algorithm de- 
pends on the complexity measures (the input size is not consid- 
ered). In the course of work, the most computationally complex 
parts of the arithmetic coder algorithm are determined. Perfor- 
mance optimization of their software implementation is carried 
out. Codecs with the new algorithm compress photo and video 
records obtained through the remote probing of water bodies 
without frame-to-frame difference. 

Materials and Methods. In the presented paper, a selection of 
satellite images of the Azov Sea area was used. At this, the soft- 
ware algorithm of the arithmetic coder was optimized; a theoreti- 
cal study was conducted; and a computational experiment was 
performed. 

Research Results. The performance of the software implementa- 
tion of the arithmetic coder is increased by the example of the 
VP9 video codec. Numerous launches of reference and modified 
codecs were made to measure the runtime. Comparison of the 
average time of their execution showed that the modified codec 
performance is 5.21% higher. The overall performance improve- 
ment for arithmetic decoding was 7.33%. 

Discussion and Conclusions. Increase in the speed of the latest 
digital photo and video image compression algorithms allows 
them to be used on mobile computing platforms, also as part of 
the onboard electronics of unmanned aerial vehicles. The theoreti- 
cal results of this work extend tools of the average-case complexi- 
ty analysis of the algorithm. They can be used in case where the 
number of algorithm steps depends not only on the input size, but 
also on non-measurable criteria (for example, on the common 


RAM access scheme from parallel processors). 


” The research is done within the frame of RSF project no. 17-11-01286. 


™ E-mail: roman.arzum@gmail.com 
“ PaOota BBINOHeHa B paMKax mpoekta PH® 17-11-01286. 


Beedenue. IIpeqnoxKeHHbii B cTaTbe ObICTpbIM MporpaMMHBIii 
aIrOopHTM apw@MerTHyecKoro KOJMpOBaHHA WpeyHasHadeH JIA 
OKaTHA WM®poBEix Uu300paxeHui. IloKa3aHo, KaKUM OOpa30M 
CJIOXKHOCTb aIFOpUTMa apuMeTHYeCKOLO KOJepa 3ABUCHT OT KpH- 
TepHeB CJIOXKHOCTH (Mp 3TOM pa3Mep BxOya He yauTHIBaeTcs). B 
Tpoljecce paOoTs! ompeyesIeHbI HavOosIee BEIYHCIIMTeIbHO CJIOXK- 
Hble YacTH aropuTMa apwdMeruyeckoro Kojepa. BpmouHeHa 
ONTHMU3AlIMA TPOW3BOMTeIbHOCTH WX TIpOrpaMMHOW peasi3a- 
IMM. 

Koylekvl C HOBbIM asIrOPHTMOM CKUMAarIOT Oe3 yaeTa Me?xKKaJIPOBOH 
pa3HUIIbI (OTO- MH BUCOMAaTepHasIbI, NOJYYeCHHbIe pH TWcTaHun- 
OHHOM 30HJIMPOBaHHH BOJIHbIX OOBEKTOB. 

Mamepuazl u memoooi. B tpeyictaBsieHHow HaydHon padote uc- 
TIOJIb30BaHa TOOOPKa CIYTHHKOBbIX CHHMKOB akBaTOopHH A30B- 
ckoro Mops. IIpv sToM onTHMH3HpOBaH TIporpaMMHBbIM asIropHTM 
apH@MeTHYECKOLO KOepa, IpOBeJIeHO TeOpeTHYeCKOe UCcIIeIOBa- 
He, BbIMOJIHCH BEIYHCJIMTeJIbHbIM IKCIEPHMEHT. 

Pezyiemamol ucciedoeanud. YBerM4eha TPOW3BOJMTeJIbHOCTb 
TIporpaMMHOU peasM3alun apudMeTH4eckKoro KOyepa Ha TpuMepe 
BuyeoKoyeKa VP9. Jia u3MepeHia BPeMeHH BBITIOJIHCHHA TIPOH3- 
BeJICHbI MHOTOUNCICHHbIe 3allyCKH 3TAIOHHOrO HW MOMUIMpO- 
BaHHOrO KOIeKOB. CpaBHeHie cpeyHero BpeMeHH HX MCHOsIHeHHA 
MOKa3asO, YTO MPOW3BOAUTeIbHOCTh MOJMMUIUMpOBaHHOrO KoeKa 
Ha 5,21 % Benue. IIpupocr oOnjei mpon3BoquTebHOCTH JV1A 
apH@MeTHYECKOLO JeKOAMpOBaHuA cocTaBu 7,33 %. 

O6cyorcoeHue u 3aKiuoueHud. YBelIM4eHHe CKOPOCTH paOoTHI HO- 
BeHIUMX arOPHTMOB KATA WAdpoBbIx (POTO- U BUACOM300paxKe- 
HWM MO3BOJIACT TPHMeHATb UX Ha MOOWJIbHbIX BbIYMCIIMTCJIbHBIX 
TlaT@opMax, B TOM YHCIIe B COCTaBe OOPTOBOM 3IEKTPOHUKU Oec- 
TIMJIOTHBIX JleTaTeIbHbIX alinaparoB. Teopetmueckve pe3ysIbTaTbl 
JaHHOW padoTbI paciiMpsAKOT MeTOJIbI aHasIN3a CJIOXKHOCTH asIro- 
PUTMa B cpeyHeM cyry4ae. OHH MOryT HCHOJIb80BaTbCA B CHTyallMHu, 
KOry{a KOJIMYeCTBO WlarOB aJIFOPHTMa 3aBHCHT He TOJIbKO OT pa3Me- 
POB BxOjJla, HO MU OT HEH3MepUMBIX KpuTepHves (HalipuMep, OT cxe- 
MBI OOpallleHua K OOMeH OlepaTHBHOM MaMATH CO CTOPOHBI! Ma- 
PAJUICIIBHBIX MPOLeCcopos). 
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Introduction. Monitoring of the condition of the water body is often carried out through unmanned aerial vehi- 
cles (UAV), conducting aerial photography in the visible and infrared bands. The performance capabilities of the mo- 
bile camera, which the UAV has, impose a number of restrictions on the equipment that processes and stores the foot- 
age until the UAV returns. Specifically, the following factors should be considered. 

1. Energy efficiency of the equipment that encodes the footage since it depends on the off-line operation time of 
the UAV. 

2. Coding gain of photo and video data during the flight. High resolution images occupy a significant amount of 
the software memory, and this limits the amount of information that a UAV can accumulate. 

Virtually all equipment for photo and video hardware supports the most common codec for compressing JPEG 
images. However, it is inferior to the most modern HEVC and VP9 codecs, which not only support video sequences, but 
also allow for better compression of single images. Thus, GoogleVP9 shows similar JPEG visual quality by the SSIM 
metric (structural similarity) and at the same time compresses images by 25—34% stronger [1]. 

Compared to the same JPEG, HEVC codec enables to improve the compression rate by 10-44% in PSNR metric 
(peak signal-to-noise ratio) [2]. However, a high compression ratio with a smaller bitstream size results in greater com- 


putational complexity of the HEVC and VP9 codecs [3, 4]. The architecture of the JPEG codec is presented in general 
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in Fig. 1. 
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Fig. 1. JPEG codec flow chart 


Blocks 
8x8 


Here, VLC means compression through variable-length codes (variable length coding). The input frame is divid- 
ed into blocks of fixed size (8x8). Each of them is subjected to direct discrete cosine transform (DCT), quantization of 
transform coefficients and subsequent entropic compression using the Huffman algorithm [4]. Discrete cosine transform 
is performed in the integer form [5]. Since the adoption of the JPEG standard in 1992, a number of fast algorithms have 
been developed. They allow the conversion to be carried out entirely in the CPU registers. Huffman’s entropy compres- 
sion 1s also not a complicated problem, so, even mobile processors in program mode can compress and decode JPEG 
images [6]. 

HEVC [6] and VP9 are fundamentally similar, and they are hybrid block codecs with splitting a frame into 
blocks of indeterminate length, intraframe prediction, discrete transform, and subsequent filtering to eliminate blocking 


artifacts. The blockdiagram of the HEVC codec is shown in Fig. 2. 
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Fig. 2. HEVC codec flowchart 


Here, the abbreviation SAO denotes the sample adaptive offset. In addition to the above image reconstruction al- 
gorithms, both codecs use context-adaptive binary arithmetic entropy coding, which is much more complicated than 
Huffman compression. With a high level of visual quality, it 1s arithmetic coding that occupies a significant part of the 
total decoder operation time. In general, the flowchart of the arithmetic coder is shown in Fig. 3. 
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Fig. 3. Arithmetic coder flowchart 


Codecs of a new generation are supposed to be used, among other things, for compressing images obtained 
through remote sensing of natural objects. In this case, it is necessary to optimize the operation of hybrid block codecs 
to ensure the processing of photo and video UAV data in real time. It is necessary to put more focus on the optimization 
of binary context-adaptive arithmetic coding, due to its essential complexity. 

The objective of this study is to speed up the operation of the arithmetic decoder on mobile processors of the 
ARM architecture. This will improve the performance of Google VP9 video codecs and will allow for their application 
in the UAV on-board electronics for remote sensing of water bodies. The use of enhanced tools for compressing photos 
and videos obtained from aerial photography enables to increase the amount of stored data, improve their visual quality 
and resolution. 

Main Part. Core components of a modern video codec include the binarization of syntactic elements of the 
bitstream and the adaptive binary encoding of these elements into the bitstream. This stage is not safe for the 
vectorization and parallelization, but it can be optimized through the statistical analysis of the input data. The main way 
to predict the running time of a program is to analyze the complexity of the corresponding algorithm. Distinction is 
provided between complexity at the best, worst and average case. = x is input data of A algorithm which is used to 
calculate the output of y algorithm. The time function of the algorithm is denoted by C} (x), and the memory cost 
function — by C7? (x). In the worst case, we will call T,(n) = maxy,j-n C4 (x) functions of numeric argument the time 


and space complexity A 
Sa(n) = maxyx =n C3 (x). 
Consider a finite set of n size inputs: 


Xp = tx: [lll = nh. 
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Vx € X, corresponds to the probability: 
P(x) € 10,1): dixeXn P(x) = 1. 

The following expectation is called the average-case complexity: 

I, = yixeXy P, (x) CA (x), 
Sa = yixeXy P, (x) C3 (x). 

The described approach is classical for analyzing average-case complexity of the algorithm, and it is described in 
detail in [7, 8]. Note the reasons why the application of this method in practice may be difficult or impractical. 

1. Difference between the number of steps of the algorithm in theory and the number of processor cycles re- 
quired to perform the step in practice. Thus, the majority of modern central processors produce addition and multiplica- 
tion per cycle, while the remainder of the division is calculated for dozens of cycles. 

2. Hardware features of the memory system. Modern computers use a multi-level memory hierarchy. Its compo- 
nents operate at different speeds. Memory calls take considerably longer than register operations. 

3. Optimizing compilers and hardware planners. When building executable files, optimizing compilers transform 
dramatically the code without changing the program state engine. CPU hardware schedulers change the procedure of 
executing instructions for greater performance and predict conditional transitions, while cache controllers read from 
memory in blocks. 

4. In case of the software algorithms implementation on the general-purpose processors, software components 
interfere. For example, the task scheduler shares processor time, and parallel processes that have multiple threads can 
run on a variable number of processor cores. 

The proposed modification of the algorithm complexity serves as a theoretical addition to practical tools for 
measuring performance, such as, for example, profiling and instrumentation of the program code. For the first time, the 
method of splitting the algorithm inputs into complexity classes was presented in [9]. Consider A algorithm and a set of 
all possible inputs: 

G: {91 G2 + }, 
and also all possible samples from G, different in size and composition: 
GiGi» Gi 3. 
A set of criteria for the complexity measure of the algorithm implementation (for example, the number of processor 
cycles, run-time, etc.): 
X;:g; 2 R. 
A set of complexity measures: 
ACG Oy. cc} 
This set has the following properties. 
1. VWa,,«,€ A: «x,4#«, — all A elements are different. 
2. W«,€ A splits G into the set of equivalence complexity classes 
G(a;) = {gi N gf ..}- 
3. All samples from G(a;) are equally complex: 
a; € A, gi € G(a;),a;: 9% > 7k € [0,[1GCa,)II],7% € R. 

Thus, all elements of A are different, and they can be reordered so that the complexity function of the criterion is 
nondecreasing throughout the whole set of criteria. The expected complexity is similar to the estimates of the average- 
case complexity of the algorithm for the discrete and continuous probability of complexity. For the discrete case: 

R(A) = diajeA TY; Pi> 
for the continuous case: 
R(A) = J, rdF(r). 

The method in question is applicable for analyzing the complexity of an arithmetic codec. The process of entro- 
py compression [10, 11] can be divided into component parts. 

1. Binarization or conversion of the coded character (syntactic element of the compressed bitstream) into a string 
made up of zeros and ones (bit string). 

2. Context modeling for compressing syntactic elements in the normal mode. This step is not performed for syn- 
tactic elements whose statistical distribution is close to normal, and they are coded in the bypass mode. 

3. Arithmetic coding of a bit string. 

Consider in more detail the arithmetic decoding scheme of the Google VP9 codec, namely the part that is associ- 
ated with the subexponential coding of syntactic elements. 
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Imagine an algorithm of the subexponential coding in general [12]. The first step includes the calculation of the 
variables: 
_ ena 2" 
~ (log, n|:n = 2*, 
Fe O:n < 2* 
b—k+1:n > 2*, 
where k is parametrical value, it is 4 for the Google VP9 codec. 

In the second step, the unary code u(u + 1) bit is complemented by n low-order bits. The code length 1s equal 
to: 

k+1:n< 2k 
2|log,n| —k+2:n > 2*. 

Hence, literal decoding is reduced to decoding the bits that compose it in the loop. To optimize the performance 
of this algorithm, it is important to know the probability distribution of literal lengths. Literals occupying the largest 
number of bits in a compressed bitstream (such as inverse transform coefficients and motion vectors) are coded in se- 
ries, so there is a high probability that the distribution of literal lengths in a compressed bitstream will be constant with 
multiple repeats of elements of the same value. To test this hypothesis, experimental data on the distribution of literal 
lengths in a set of satellite images of the Sea of Azov is acquired (Table 1). 


ut+l+n=| 


Table 1 


SO 





/Probabiiy.% «oss fo | 7x8 | amas [0 | nae 


The literals of 3, 4, and 6 bits in length are most probable. The maximum possible literal length for this sequence 
is just 6 bits. This fact is essential for software optimization of the function of subexponential decoding of a literal. As 
part of optimizing a real codec, the run-time criterion is of prime interest. To obtain the set of difficulties: {7,...,7%}, 
we will profile the program performance. 

A set of unique elements (RX) will make up a set of criteria of the run-time complexity (A ). 

The following approaches that are based on the data obtained are applied for the optimization. 

1. Memoization of calculating the literal length to decode a series of literals of equal length. 

2. Unwinding of a cycle of subexponential decoding of a literal. 

3. More efficient algorithm for calculating the number of bits in a literal. 

4. More efficient use of the processor registers immediately within the arithmetic decoding function. 

The implementation of points 1 and 4 is fairly obvious; therefore, we consider in more detail points 2 and 3. 
Compiling of a decoded literal by bits decoded from a compressed bitstream occurs within the subexponential decoding 
function. In this case, the bottleneck is a loop with varying number of iterations. It can be replaced with a switch-case 
set without break at the end. This technique is known as the Duff’s device method. It allows replacing several loop iter- 
ations through sequentially execution of the instructions without the need for conditional transitions. The bit shift 
amount is a constant that does not need to be read from the register - loop counter. 


Code Listing 1: The original literal decoding function 


static int vp9 read _literal(vp9 reader “br, int bits) 
{ 
int z = 0, bit; 
for (bit = bits — 1; bit >= 0; bit —) 
z |= vp9_read_bit(br) << bit; 
retum Z; 


j 


Code Listing 2: Modified literal decoding function 


static int vp9 read _literal(vp9 reader “br, int bits) { 
register int z = 0; 
switch(bits -1){ 

case 6: z |= vp9_ read(br, 128) << 6; 

case 5: z |= vp9_read(br, 128) << 5; 

case 4: z |= vp9_read(br, 128) << 4; 





Arzumanyan R. V. Arithmetic coder optimization for compressing images obtained through remote probing of water 


bodies Ap3ymanan P,. B. Onmumu3auua apudmemuueckozo Kodepa O/1A CHCAMUA UZO0Opastcenuu, NOIYUeEHHbIX NDU OUcmManuUuonnom 
———————— Eee a SE 


case 3: z |= vp9_read(br, 128) << 3; 
case 2: z |= vp9_read(br, 128) << 2; 
case 1: z |= vp9_read(br, 128) << 1; 
case 0: z |= vp9_read(br, 128); 


break; 
j 


retum Z; 


j 





Another bottleneck is the calculation of the number of literal bits in the while loop [13]. This is a worse solution 
because the number of loop iterations is unpredictable. Instead, a fast bit-counting algorithm was used [14, 15], which 
performs the calculation for a fixed number of steps without conditional transitions. 


Code Listing 3: Fast counting of the number of bits in a literal 


Unsig ned intv; // 32-bit argument 
Register unsig ned intr; // variable for the number of bits 
register unsigned int shift; 

tr =(v > OXFFFF) << 4; 

V>rr= PF; 

shift = (v > OxFF) << 3; 

v >>= shift; 

r = shift; 

shift = (v > OxF) << 2; 

v >>= shift; 

r = shift; 

shift = (v > 0x3) << 1; 

v >>= shift; 

r = shift; 

rl (v >> 1); 





Numerous starts of the reference and modified codecs were made to measure the runtime. In this case, their av- 
erage run-time was compared. It has been found that the performance of the modified codec is 5.21% higher. The in- 
crease in overall performance for arithmetic decoding was 7.33%. 

Conclusions. The operation of the arithmetic coder as a component of the video codec has been optimized 
using the example of the Google VP9 standard. To solve this problem, a modification of the method for analyzing the 
average-case complexity of algorithm has been proposed. The approach is based on the partitioning of the set of inputs 
into equivalence complexity classes. The considered method enables to predict the average-case complexity of the 
algorithm when the number of steps of the algorithm and the time of its execution depend on difficult-to-measure 
parameters, which is typical of the context-adaptive arithmetic coding. The proposed method has been applied to 
optimize the speed of the arithmetic binary coder (using the example of the Google VP9 codec) for the image 
compression problems obtained under remote sensing of water bodies. The research results make it possible to apply 
advanced methods of compressing photo and video data obtained through the aerial photography of water bodies. Thus, 
it is possible to increase an amount of the accumulated data, to improve the visual quality and resolution of the footage 
by 25-34% (according to the SSIM visual quality metric) and to increase the speed of the arithmetic coder by 7%. 
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Introduction. The solution to the problem on the centroidal 
motion control synthesis (guidance problem) of an unmanned 
aerial vehicle (UAV) with long-range capabilities in the boost 
phase is considered. Control condition requires optimum fuel 
consumption. The principle of dynamic programming consid- 
ering the restrictions to the vector modulus of the thrust output 
is used to solve the problem. The implementation of terminal 
guidance requires the formation of control as a function of the 
object state at the end of the ascent phase. The attainment of 
these boundary conditions determines the further transition to 
the ballistic flight phase. 

Materials and Methods. Bellman’s principle of dynamic pro- 
gramming is the most reasonable from the point of view of the 
implementability of the computationally efficient on-board 
algorithms and the solution to the problems in the form of 
synthesis. With natural scarcity of thrust and energy resources 
on board, this principle enables to obtain solutions free from 
the switching functions. In this case, the optimal control is a 
smooth function (without derivative discontinuity) of the cur- 
rent and final parameters of the UAV. 

Research Results. A new algorithmic method for the synthesis 
of terminal motion control is developed. Its difference is that 
the UAV movement control in the ascent phase is formed by 
the function of the motion actual and terminal parameters. 
This ensures movement along an energetically optimal 
trajectory into the given region of space. The problem solution 


results enable to build closed terminal guidance algorithms for 


” The research is done within the frame of the independent R&D. 


Beedenue. CratTba MOcBAIIeHa pellleHHIO 3aya4H CHHTe3a 
ylipaBJIeHHa JIBYXKeHHeM IleHTpa Macc (3aqa4a HaBeyleHHs) 
OecHHIOTHOrO JIeTaTesbHOrO amnapata (bJIA) c OonpmoL1 
JaJIbHOCTbIO MOJIeTa Ha pa3srOHHOM yyacTKe. YcuIOBHe yupaB- 
JICHHA: ONTHMAJIbHbIM pacxoy TomIMBa. Jia pelieHua 3a7aqn 
HCHOUIb3yYeTCA TIPHHUMI THHAaMMYecKOrO MporpamMMupoBaHuA 
C y4eTOM OrpaHHyeHHt Ha MOJlyJIb BEKTOpa TAM JBUTaTessA. 
Peamn3allua TEPMHHAJIbHOrO HaBeyeHuA TpeOyeT CPOpMupoBa- 
HUA yiIpaBleHud Kak (yHKUMM COCTOAHUA OObEKTAa B KOHIIC 
yuacTka BbIBeyeHua. JlocTuxKeHHe 3THX TpaHH4HbIX ycuOBUi 
olpeyeiaeT TaibHeiiiMu mepexoy kK OammMcTHYecKON da3ze 
moyleTa. 

Mamepuaazoi u memoooi. Wpunuun uHamuyeckoro mporpam- 
MupoBaHna beswimMaHa ABJIAeTCA HaMOOIIee pallMOHaJIbHbIM C 
TOUKH 3PeCHHA peallMu3yeMOcTH IPeKTHBHBIX B BbIYHCIIMTEIb- 
HOM OTHOMICHHH OOPTOBHIX aJITOPHTMOB UH pellleHHA 3aadu B 
(bopMe cuHtTe3a. IIpu ecTrecTBeHHOM OrpaHHyeHHOCTH BeJIM- 
4HHbI TATH HW 9HEPreTH4YeCKMX PecypcoB Ha OoOpTy MaHHbIi 
TIpHHUMM MO3BOJAeT MOJYYINTb pellleHua, He coOseprxkaHe 
(@YHKUMU MepekHOueHHA. ONTHMasIbHOe yiIpaBlIeHve B 3TOM 
ciyyae ABJIAeTCA TaKOU PyHKUMet (Oe3 pa3spbiBa IpOu3BO- 
HOM) TeKYUIMX UW KOHCUHBIX MapamMeTpos BJIA. 

uccIeOOBAHUA. 

clocoo 


Pezylbmambl Pa3pa0oTaH HOBBIM 


asIrOopHTMuYecKHi CHHTe3a TePpMMHasIbHOrO 
yipaBsieHuad BMyKeHuA. Ero oTIM4Me B TOM, 4TO yiIpaBsIeHHe 
mpwxkeHnueM BbJIA Ha pa3roHHOM y4acTKe TpaeKTopHu 
(bopmupyetca PyHKIMeH TeKyUIMX MH KOHCYHbIX MapaMeTpoB 
mBwxKeHua. Takum o0pa30M oOecneuMBaeTCA JIBYHXKeHHe TO 
9HepreTHyeCKH OMTHMaJIbHOM TpaeKTOpuH B  3aaHHyto 
oOmacTb wWpoctpaHctBa. Pe3yiIbTaTbl pelleHua 3aa4u 
TIO3BOJIAIOT CTPOUTb 3AMKHYTbIC aJITOPUTMbI TePMHHaJIbHOrO 
HaBeyleHHA JIA pa3rOHHOrO yyacTKa TpaexTopuu BJIA c 
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the boost phase of the UAV trajectory with long-range 
capabilities. Such algorithms have good convergence and 
injection accuracy due to the prediction of parameters during 
the flight at a shorter time interval. 

Discussion and Conclusions. The most preferred is_ the 
principle of dynamic programming. It should be used when 


solving the problem on the centroidal motion control synthesis 


XOPOMIeH CXOJMMOCTbIO MH TOUHOCTHIO BbIBeCJICHHA 3a CUCT 
IIPOrHO3HpOBaHHA TapaMeTpoB B porecce oseTa Ha 
COKpalljarollleMca HHTepBasie BPeMeHH. 

O6écyacdenue u 3akmoueHua. Handomee lpexqnNOuTUTeIbHbIM 
IIpecTaBsAeTca IIpHHIMH MHAMUYeCKOrO 
porpamMmMupospanna. MMeHHo ero cileyeT MCHOJIb30BaTb TIPU 
pelleHHu 3aja4dn CHHTe3a ONTMMAIIbHOrO MO pacxoyy TOMIMBa 


ylipaBJIeHHa JIBHXKeHHeM IleHTpa Macc (3aqa4a HaBejleHHs) 


(guidance problem) of the UAV with long-range capabilities BJIA c OoMbUIO WaIbHOCTbIO HoseTa Ha pa3srOHHOM y4acTKe. 


in the boost phase. 
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Introduction. Currently, the capabilities and scope of applicability of the unmanned aerial vehicles (UAV) 
have increased significantly. This is primarily due to the UAV flying range extension. Vehicles with a ballistic phase of 
flight should be provided with the control in various areas, including the UAV positioning. For the control design, it is 
reasonable to apply the principle of terminal guidance. The control should be developed as a function of the terminal 
motion variables, and not as a temporal function. 

A great many publications are devoted to this problem solution, but the task described above holds relevance. 
In particular, the development of highly computationally efficient algorithmic methods of terminal guidance 1s of inter- 
est. At the same time, the features universal for various types of launch vehicles should be considered. They have adap- 
tive features and in a certain sense meet the optimality requirements. 

Considering the guidance task, such a control, which uses the minimum amount of fuel, is valuable. The syn- 
thesis of optimal control is based on the use of the Bellman dynamic programming method [1]. 

Materials and Methods. The solution to the problem of the synthesis of optimal control of the UAV motion in 
the boost phase of power-on flight has been studied in many papers. However, the optimal control solution obtained in 
most cases is reduced to the implementation of a time or parametric program. We will solve the problem of finding the 
optimal control in the UAV ascent phase in the following formulation. The UAV motion parameters are known as Xo, 
Yo, Zy) coordinates of the current path point, obtained through solving the problem of navigation. X,, Y,, Z, parameters of 
the final point satisfy the boundary state that fixes the transition to the ballistic phase of flight determined by the hyper- 
surface in the phase space. The condition of S; (X,., Y,, Z,.) = 0 1s satisfied by a whole set of finite parameters. It is re- 
quired to synthesize optimal control in the task of UAV guidance, which ensures its transfer from the initial state to the 
hypersurface of the final conditions. 

The optimality criterion is the amount of fuel consumed in the boost phase: 


m(t) = ri(t)dt (1) 


where m(t) 1s fuel mass flow rate. 


The following system of differential equations is adopted as a mathematical model of the centroidal motion: 
R®=Ve), (2) 
V(tV=WOEO+ 87), 
where R(f) is radius-vector, Z(7) is terrestrial attraction vector, W(t) is module of control acceleration vector. 
The module of the control acceleration vector (W(t) ) is a specified temporal function, and it is determined by the 


UAV engine performance. Unknown is the unit vector of control acceleration — E(t). When solving the optimization 


problem, it will determine the required properties of the UAV guided motion in the boost phase. There is no need in the 
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thrust vector control (in a sense, this is apparent acceleration amount — W(t) ). A more rational approach is to maximize 


this value for the engine booster. This ensures the application of the solutions obtained for the case of using solid fuel 
engines. The UAV engine performance is stable enough, and at m(t) = const constant fuel-flow rate, the optimal con- 


trol is determined by the combustion duration. In this case, the functional (1) will be a function of the upper limit of 
integration. Hence, the problem of minimizing the fuel amount turns into an equivalent task of minimizing the flight 
time, and speed-of-response will be the optimality criterion. The task of the synthesis is to find the orientation of the 


thrust vector of the UAV engine, which is determined by the direction cosines of the thrust vector ( P(t) ) as a function 


of the current parameters and the final state. 

We admit two assumptions. The first is as follows. Since a large part of the boost phase lies outside the dense 
atmosphere, we will not consider the angle rate of evolution of the thrust vector in space (longitudinal axis of the UAV) 
as constraint. The second assumption is the following. Assume that the UAV motion is passing in a predetermined 
plane. 

The UAV final state for the transition to the ballistic phase of flight is fixed by satisfying the following bound- 
ary condition [2]: 

Si = (Vy, =x )[V (y, -y,)-V,, ie —%, )|- 
XX TV, (3) 


U1 y k k 


Here, x,,),,X,,¥, are, respectively, the coordinates of the initial point and the point of the starting of the UAV opera- 


tion on the final path segment. The current value of the boundary condition (S,,(¢,) ) is a measure of non-compliance 


with the condition (3). A mathematical notation for this condition corresponds to the hypersurface, which is a smooth 
function of phase coordinates and describes the entire family of possible UAV ascent trajectories [2]. 

The boundary conditions are specified for the central field of attraction. For this case, the equations of the UAV 
motion will have the following form: 


x, =%as 
Tt 
{= XP COSC, 
"x =x, 4) 
: , 
Tt 
Xx, =-—x,+—Pcosa,. 
3 3 2 
r m 


Here, =, = f-M, 1s constant of the central field of the earth’s attraction equal to the product of the gravitational con- 


stant (f) and the mass of the Earth (M,); r=(x; +x, )'", x, =x, x,=V,,x,=y,x,=V, ; P is thrust vector value; X4 is aero- 


x9 "3 
dynamic drag force; Y, is aerodynamic lift. 
Since the UAV is fitted up with a steerable thruster, its direction cosines will determine its pointing. In this case, 
the handling constraint will be determined by the ratio: 


|PO|={[PO cosa, (OF +[PO cose, One <|P°(o). (5) 


The handling constraint (5) is “hypersphere restriction”. This implies the solution in which optimal control is 
not piecewise constant, without switching. In this case, the time optimal ascent trajectory in the phase space has no “an- 
gles” and no discontinuities of the derivative. The time optimal control will be a nonlinear, continuous function of the 
boundary conditions (3) [3]. 

The synthesis task is formulated as follows. The control object from an arbitrary current state, taken as the ini- 
tial one and determined by the current value (Sj;,) at the time (f), 1s transferred to the hypersurface of the boundary con- 
dition S|, =0 at the time (¢,). This takes into account the control acceleration amount (engine thrust module). It is re- 
quired to find the optimal control in the form of synthesis, which provides such a transfer in the shortest possible time. 

The plant state at the final instance satisfies the boundary condition S,[X(t,)]=0 and determines the moment 


of transition to the ballistic phase of flight. 
In accordance with R. Bellman’s dynamic programming principle [4], the necessary and sufficient condition 
for optimality of the formulated problem will be the ratio: 


nip] 1.0) =-], (6) 


ueU i=l 


Considering the mathematical model of the object specified by the system (4), the dynamic programming equa- 
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tion [4] with the optimality criterion, is determined by the relation: 


0 0 0 0 
np a [Fe - tp |S +x ca -(3s-22)2 |-+. (7) 


a 4 
Oe Ox, r m Ox, Ox, r m Ox, 














where P, = Pcosa,,P, = Pcosa,. 


Due to the methodology of Bellman’s dynamic programming, the minimization can be carried out through the 
application of the Schwartz inequality [3] to the relation (5). This will significantly simplify the solution to the optim1- 
zation problem. Then, the expression (7) with regard to inequality (5) will be determined by the relation: 














2 2 1/2 
0 0 0 0 
gi oe | | (8) 
wU | m ~ Ox, Mm” OX, m OX, m OX, 
In this case, the Hamilton-Jacobi equation will be as follows: 
OT’ 1%, OT” 








Tt 
0 0 
x ee, ye (LX 
Ox, I Ox, 1 


: ae ; 9 
1 oT 1 oT ) 
x} | — +| — =I, 
m Ox, m OX, 


The equation (9) is solved considering the given boundary condition: 























(X15 Xp> Xo %,) ES - (10) 
The calculation of the partial derivatives of S for phase variables (x,,x,,x,,x,) gives the following dependenc- 
es: 
os” 
— (-2x,x,x, nee Wee 90 ee ee 
1 
os* ; 
Fs = (2x; Ky 2X Ae 2K Xi ON Xa, ) 
<a (11) 
os ; ; 
Ay = wo KH 2L Me Ht HW, ), 
3 
os* ; 
rs = (2x; Ki ALN LK KG PAH, PX NX ). 
4 
Through substituting the relations (11) into (9) and performing simple transformations, we obtain the depend- 
ence: 
P 0 
{= @, SR AZ N i, oe 
r m : 
x{ (2x x, LU NN LK Ns NN +X,%5%,, ) + (12) 
or: 
+(2x7%, SZ Oh a, ) lian — —] 
1 
with the boundary conditions 
T°’ (S,) =0, at (x,, x,,%,,x,) €5,,. (13) 
The expression in curly brackets in (12), denoted by S’(x,), can be written in a compact form: 
oT” 
S’(x.)- =-—]. 14 
(x; ) as (14) 


The equation (14) with boundary conditions (13) can be solved in various ways, for example, by the method of 
characteristics [5]. However, it is more rational to use the following method. 
The expression (11) defines the optimal control structure [6]: 


ae (x) = —PRe”" (x), 


oe - (15) 
We (x) — WE" (x), 
where E(x) is unit thrust vector. 
The ratio (15) in the scalar form: 
P” = Pcosa,, 
P™ = Pcosa,, ~ 
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ar® 
Ox, 
COSQ, = rT, Ung 2? 
OT” oT” 
|| Sy 
Ox, Ox, 
ax, 
COSQ, = 


2 2 1/2 * 
OT” OT” 
—— + —__ 
Ox, Ox, 


The equation (9) and relations (17) can be put in a more convenient form for the following reasons. According 
to the formation of a set of perturbed trajectories, through varying the control, it is possible to construct a hypersurface 
with equal ascent time, that is, an isochronous surface, in the phase space. Indeed, for each phase trajectory point, we 
will calculate the optimality criterion value: J(t,), J(t),...,/(t). 

Hence, we obtain a set of trajectories for each fe [t,t]. Owing to the continuity of the dependence of x(¢) and J 
on the variable control, a set of trajectories forms a surface in the phase space (X). This boundary surface formed by the 
set of vectors x[t;,/(t;)], 1s convex and smooth. For transition conditions to the ballistic flight phase, the isochronous 
surface has a tangency point with the hypersurface of the boundary conditions. 

Under the qualitative implementation of the optimal control, the distance in the phase space between the hyper- 
surface of the boundary conditions S;=0 and the isochronous surface T (x, R,, t) = 0 will decrease. During some time, 
the two surfaces will have a common point (Fig. 1). 


Sik (x, Ry, t) = 0 
Xi 


SA 


T (x, Ry, t) =0 


AG 


Fig. 1. Geometric interpretation of hypersurfaces in phase space under guidance in boost phase 


At the point S; and 7°, they have a common tangency, that is, their gradients coincide in the phase space. In this 
way: 

- existence of the optimal trajectory x,,,(¢) and optimal control u,,, (¢) is noted; 

- termination of the boost phase of the trajectory (¢ = ¢,) and the transition to the ballistic flight phase are deter- 
mined. 

At the tangency point, the isochronous surface and the hypersurface of the boundary conditions have a com- 
mon tangency and a normal. Mathematically, the condition for the existence of a common normal is determined by the 


expression [2]: 
0 
or") _ 4{@.\. a 
Ox, Ox, 


where A is constant determined from the analysis of the convexity of both hypersurfaces. 
From the relations (11) and (18), we transform (17). We obtain the expression for the direction cosines of the 
thrust vector as a function of the current and final motion parameters: 
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4- (—2x,x,%, eee 2x; X, oe ) (19) 


2 2 9 
2 2 2 1/2 
(2x, x3, Nay Nag 20 —x{X,, ) + (2x,x LN Ky Nyt LN XN Py Ny Nyy +X,%;%,, } | 


2 
(2x; Ky 2X NX, PLN HX NHK, FX NX, ) 


B- (20) 


2 2 1/2 ? 
2 2 2 
| (2x, x,3, +X,X,X,, +2X,X, —X4%,, ) + (2x,%, LN kh 2K hin NN, +X,%5%,, } | 


where 4=cosa,, B=cosa,. 
From the relations (19), (20), it is easy to obtain a parameter natural for this UAV type to determine the orien- 


tation vector of the control acceleration (in a certain sense, of the thrust vector), the pitching angle: 
COS, _ Coes = 20K N, 2H X,, PX AN, FH MX, ) (21) 


O(S) = arctg 
(-2x,x,%, to tO a, ) 





O 
COS QL, 

In our case, the UAV motion takes place in the given plane, and the yaw angle is  (S') = 0. Using the proposed 
methodology for the synthesis of terminal optimal control, an algorithm for calculating the pitching angle was devel- 
oped [2], and the computational simulation of the UAV flight was carried out [7]. 

Research Results. The numerical studies have been conducted using the software that implements the pro- 
posed method. We are talking about the algorithmic software for terminal guidance of ballistic aircraft based on the 
solution to boundary problems of ballistics. The corresponding computer program was registered in 2013. 

When modeling, a hypothetical accelerating tool was used with the characteristics and initial conditions for the 


UAV launch ascent, given in [8]. 
The initial conditions of the simulation and the calculation results of the early ascent are given in Table 1. 


Table 1 
Initial conditions for UAV launch ascent 


xin 
5.6323 —0.0000 322.6757 90.00 0.00 0.000 
5.6323 11.1626 322.6756 90.00 1.98 —0).048 
21.4565 69.4705 322.6727 69.34 12.13 —7.881 
155.6274 259.0827 322.6335 60.00 35.9] 0.011 


368.3725 464.9493 322.544] 50.02 43.59 —2.076 
662.2567 679.2317 322.4053 43.44 42.79 —2.582 
1044.7345 915.1428 322.2176 38.74 40.11 —2.677 
1429.1498 1131.2610 322.0334 36.01 37.94 —2.501 





Fig. 2 shows the pitching-angle profile in the boost flight phase. 
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Dependence of pitching and attack angles on apparent speed rate 





0 1 000 2 000 3000 4000 3000 6000 F000 


— pitching angle — attack angle 


Dependence of UAV entrance angle on apparent speed rate 





0 1 000 2 000 3000 4000 3 000 6 000 F000 


— UAV entrance angle 
Fig. 2. Variation of UAV pitching and entrance angles 


Fuller information on the simulation results given in [8] is as follows. 

1. The resulting control is suboptimal due to the use of boundary conditions in an analytical form. To improve 
the accuracy of the boundary conditions in the navigation algorithm, it 1s necessary to introduce an up-date algorithm, 
which is based on attracting more accurate models of the Earth’s gravitational field. Hence, the launch ascent accuracy 
will be improved. 

2. The effect of random disturbances is compensated by the adaptive properties of terminal guidance, as well 
as by predicting the motion parameters and the formation of control at a decreasing time interval each time until the end 
of the launch ascent process. 

3. The obtained estimates under modeling allowed us to rationally select the discreteness of the navigation 
and guidance algorithms and thereby limit the range of requirements for the FMC during its implementation. 

4. The implementability of guidance algorithms based on the developed methodology on modern onboard 
computers creates no difficulties. The required response rate is (1-1.5)x10° k.o./s. 

Discussion and Conclusions. Thus, the principle of dynamic programming seems to be most preferable. It 
should be used when solving the problem on the centroidal motion control synthesis (guidance problem) of the UAV 
with long-range capabilities in the boost phase. The well-known remark on the applicability of dynamic programming, 
the so-called “curse of dimensionality”, is inappropriate in the task of developing control as a function of the final state 
[5]. In addition, the use of boundary conditions in an analytical form simplifies essentially the formation of suboptimal 
control and enables to change strategically the flight missions. This scales up the applicability of this algorithm for 
UAV of various purposes. 
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